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RESPONSE  OF  STRUCTURAL  SYSTEMS 
TO  NONSTATIONARY  RANDOM  EXCITATION* 


ABSTRACT 


In  this  work,  the  response  of  lumped  parameter,  second  order  systems  to 
nonstationary  random  excitations  is  examined.  Included  is  a  brief  intro¬ 
duction  to  the  probabilistic  theory  of  structural  dynamics  and  various 
basic  concepts  required  for  subsequent  work. 

More  specifically,  the  second  order  central  moment  (covariance)  response 
of  structural  systems  to  random  excitations  is  studied.  In  the  course  of 
the  analysis,  an  approximate  method  for  the  calculation  of  system  re¬ 
sponse  to  a  class  of  nonstationary  excitation  processes  is  constructed. 
This  class  of  excitations  we  have  called  "slowly  varying"  nonstationary 
random  processes.  By  this  is  meant  that  the  nonstationary  variation  of 
the  correlation  functions  of  the  process  is  small  compared  with  the  time 
variation  of  the  impulse  response  functions  of  the  system  considered. 

It  is  shown  how  this  approximation  technique  may  be  applied  to  the 
estimation  of  inertial  loads  in  the  structural  members  of  a  payload  during 
the  launch  phase  of  flight.  Employing  previous  rocket  engine  test  data, 
the  excitations  to  the  payload  are  idealized  as  a  "slowly  varying”  non¬ 
stationary  random  excitation.  An  approximation  procedure  is  then  de¬ 
veloped  for  the  calculation  of  the  second-order  central  moments  of  the 
payload  response. 
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CHAPTER  I 


INTRODUCTION 

1  - 1  General  Nature  of  the  Problem 

This  work  is  concerned  with  the  estimation  of  inertial 
loads  in  a  payload  structure  during  the  launch  phase  of  flight. 
It  has  been  motivated  by  research  on  this  problem  with  regard 
to  the  Lincoln  Experimental  Satellite  (LES)  series  for  which 
the  Titan  III  launch  vehicles  have  been  used.  Thus,  it  is 
toward  the  LES  program  that  this  thesis  is  generally  directed. 

For  LES-5,  structural  qualification  was  based  on  a  dynamic 
analysis  of  the  payload  system  which  consisted  of  the  dispenser 
truss  and  all  satellites.  To  accomplish  this  analysis  the 
satellite  contractors  supplied  Martin-Marietta  Corporation  with 
the  relevant  characteristics  of  the  satellites.  Martin- 
Marietta  then  determined  the  system  response  to  excitations 
obtained  from  3  of  27  booster  engine  test  firings.  The  analysis 
was  carried  out  only  for  the  core  engine  cut-off  phase  of 
flight  since  it  is  during  this  period  (approximately  2  seconds) 
that  the  payload  is  subjected  to  the  most  severe  dynamic  load¬ 
ing  . 

It  was  subsequently  assumed  that: 

A-  The  dynamic  excitations  are  such  that  the  maximum 
response  of  the  payload  system^  ^  occurs  in  the  low  frequency 
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range  (less  than  50  Hertz). 

B.  The  six  components  of  rigid  body  acceleration  of  the 
booster/payload  interface  are  virtually  independent  of  the  pay- 
load  impedance  for  payloads  in  a  given  weight  category. Thus, 
in  estimating  the  low  frequency  response  of  a  satellite,  the 
acceleration  time  histories  at  the  booster/payload  interface 
were  assumed  to  be  the  basic  dynamic  inputs.  We  shall  adhere 
to  this  assumption  in  the  following. 

In  this  work,  we  shall  be  dealing  with  the  low  frequency 
response  of  a  payload  to  excitations  obtained  from  the  accelera¬ 
tion  time  histories  mentioned  above. 

1 . 2  Interface  Accelerations  Considered  as  Random  Processes 

Typical  plots  of  the  booster/interface  accelerations  are 
shown  on  page  12.  The  non-reproducibility  of  these  data  are 
to  be  noted.  Indeed,  from  the  structural  designers  point  of 
view,  these  time  histories  can  only  be  described  as  random. 

The  apparently  random  nature  of  the  data  suggest  the  following 
approach  to  the  estimation  of  structural  member  loads: 

A.  The  27  time  histories  of  the  six  acceleration  com¬ 
ponents  are  considered  to  be  ensemble  members  of  a  nonstationary 
random  process. 

B-  Statistical  parameters  of  the  stochastic  response  of 
the  payload  are  to  be  calculated  with  the  interface  accelera¬ 
tions  considered  as  the  stochastic  inputs. 
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ACCELERAT I  ON C  G 


Fig*  1  Tne  coordinate  system 
for  the  interface  accelerations 
is  shown  at  left;  (a).  A  compos¬ 
ite  plot  of  all  the  z  accelera¬ 
tion  time  histories  from  the  27 
booster  engine  test  firings  are 
shown  below;  (b). 


( b ) 
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The  most  important  response  parameters  to  be  obtained  are 

the  means  and  mean  squares  of  the  loads  on  the  payload  structural 

members.  For  if  x  is  one  of  the  stresses  in  some  member  of 

2 

the  structure,  the  expectation  value  of  x,  and  Y^  the  expecta- 

2  (2) 

tion  of  x  then  in  the  customary  notation  : 


p{  lx-^xl  1  kcJ}  1 


(1.1) 

(1.2) 


where  P|_A]  is  the  probability  of  event  A,  and  a  is  called  the 


x 

variance  or  second  central  moment  of  x.  We  thus  see  that  a 

knowledge  of  P  and  Y  allows  us  to  determine  an  upper  bound 
x  x 

to  the  probability  that  the  member  stresses  will  exceed  the 
values  for  which  failure  occurs.  Since  the  evaluation  of  the 
means  is  relatively  simple,  we  shall  be  mainly  concerned  with 
the  determination  of  the  variance  of  the  member  loads.  To 
obtain  these  it  is  necessary  to  determine  the  complete  set  of 


★ 


second-order  central  moments  of  the  payload  displacements. 

In  summary,  the  approach  to  be  presented  involves  the 

construction  of  a  "stochastic  equivalent"  to  the  interface 

accelerations  and  the  subsequent  evaluation  of  the  second 

*  See  reference  (2),  p.  24  or  Chapter  2  for  definition  of 
second  central  moments.  See  also  Section  2.8  for  a  justifica¬ 
tion  of  this  assertion. 
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moment  response  of  a  satellite.  In  general,  the  object  of 
this  thesis  is  the  development  of  a  plan  of  analysis  whereby 
the  approach  may  be  conveniently  applied.  The  two  major  aspects 
of  the  problem  are  discussed  in  more  detail  below. 

★ 

1 . 3  Statistics  of  the  Interface  Accelerations 

In  the  analysis,  three  basic  assumptions  have  been  made 
concerning  the  general  nature  of  the  interface  accelerations. 

A.  The  interface  accelerations  are  the  response  of  a 
lumped  parameter,  second  order,  linear  system  (i.e.  the  Titan 
III  launch  vehicles)  to  a  nonstationary  vector  shot  noise 
excitation . 

B.  The  intensity  functions  of  the  components  of  the 
vector  shot  noise  are  the  same  to  within  multiplicative  con¬ 
stants  . 

C.  The  excitations  to  different  normal  modes  of  the  booster 
are  statistically  uncorrelated. 

In  chapter  VI  these  assumptions  are  discussed.  Furthermore 
it  is  shown  that  under  the  above  assumptions,  the  generalized 
spectral  density  matrix  of  the  six  booster-payload  interface 
accelerations  can  be  determined  from  the  general  properties  of 
the  time  histories  and  appropriate  averages  of  the  Fourier 
transforms  of  these  time  histories. 

*  For  definitions  of  statistical  terms,  see  references  (3), 

(2)  or  (4). 
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1.4  Evaluation  of  Second  Moment  Response 


'The  techniques  for  the  calculation  of  the  response  of  a 
system  to  both  stationary  and  non-stationary  random  excitation 
are  well  established.  However,  the  evaluation  of  response 
to  nonstationary  processes  is  considerably  more  taxing  than  in 
the  stationary  case0  Indeed,  for  the  simplest  systems,  the 

analytic  expression  for  the  second  moment  response  may  be  quite 

•  ^  <5'6> 
unwieldy . 

From  Section  1.3,  the  accelerations  are  assumed  to  repre¬ 
sent  the  output  of  a  nonstationary  random  vector  process. 
Furthermore,  they  are  the  inputs  to  a  complicated  structural 
system.  Therefore,  the  need  for  reliable  approximation  pro¬ 
cedures  is  evident. 

One  such  procedure,  developed  by  Bucciarelli  and  Kuo, 
deals  with  excitations  of  the  form 

x(t)  =  g(t)n(t)  (1.3) 

where  n(t)  is  a  stationary  random  process  and  g(t)  is  a  deter¬ 
ministic  function  of  time.  When  g(t)  satisfies  certain  restric- 

( 8 ) 

tions,  Reimann's  lemma'  may  be  used  to  eliminate  the  high 

frequency  terms  in  the  expression  for  the  second  moment  response. 

This  analysis  has  been  applied  successfully  to  single  degree 

(9) 

of  freedom  systems/ 

The  present  work  is  largely  concerned  with  the  improvement 
and  generalization  of  the  approach  mentioned  above.  A  procedure 
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is  developed  for  the  approximate  calculation  of  the  second 
moment  response  of  a  payload  system  to  excitations  of  the  form 
presented  in  Section  1.3.  An  analysis  of  the  statistical 
properties  of  the  interface  accelerations  is  also  included  but 
no  attempt  will  be  made  at  a  detailed  calculation  of  the  confi¬ 
dence  intervals  of  the  estimated  statistical  parameters. 

1 . 5  Summary  of  Contents 

In  Chapter  II  we  review  the  aspects  of  the  probabilistic 
theory  of  structural  dynamics  relevant  to  this  study.  Since 
the  proposed  approximation  technique  involves  extensive  use  of 
the  Fourier  integral,  a  discussion  of  its  elementary  properties 
is  included  in  Chapter  III.  At  the  conclusion  of  Chapter  III 
we  present  two  related  theorems  which  form  the  basis  of  the 
approximation  procedure.  The  general  response  formulation  of 
Chapter  IV  leads  directly  to  the  application  of  the  theorems 
developed  in  Chapter  III  to  the  approximate  evaluation  of  the 
second  moment  response  in  Chapter  V.  Finally,  in  Chapter  VI 
we  discuss  the  statistical  analysis  of  the  interface  accelera¬ 
tions  and  outline  a  general  procedure  for  the  calculation  of 
the  variances  of  the  member  loads  in  an  actual  payload  structure. 
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CHAPTER  II 


PROBABILISTIC  THEORY  OF  STRUCTURAL  DYNAMICS 

2 .1  Introduction 

Here  we  review  certain  aspects  of  the  theory  of  structural 
dynamics  and  of  the  theory  of  random  processes  needed  to  formu¬ 
late  the  response  of  a  structural  system  to  random  excitations. 
In  Sections  2.2,  2.3  and  2.4  are  presented  the  deterministic 

equations  of  response.  Sections  2.5  through  2.8  include  a 
discussion  of  certain  results  in  random  theory.  The  chapter 
concludes  with  the  formulation  of  the  second  moment  displace¬ 
ment  response  of  a  structure  to  random  excitations  and  a  con¬ 
sideration  of  two  special  classes  of  random  processes. 

2 . 2  The  Idealization  of  Elastic  Structures 

a.  The  system  has  a  finite  number,  N,  of  degrees  of 
freedom.  Thus  the  state  of  the  system  can  be  specified  by  a 

vector,  [x] ,  whose  components  are  the  generalized  coordinates 

★ 

of  the  system. 

b.  The  inertial  properties  can  be  represented  by  a 

symmetric  matrix  of  constants,  [m]  ,  called  the  mass  matrix . 

_  (Nx  N ) 

*  For  an  explanation  of  the  matrix  notation  see  reference 

(10)  . 
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The  total  kinetic  energy  of  the  system  during  its  motion  is 

-  *  [m]  U1  . 

2 

c.  Each  component  of  the  structure  is  assumed  to  be 

linearly  elastic.  Hence,  the  total  potential  energy  due  to 

elastic  deformations  is  —  L xj  [k!  [x]  where  [k]  is  a  symmetric 

2 

matrix  of  constants,  called  the  stiffness  matrix . 

d.  Because  only  small  motions  are  considered,  the 

potential  energy  due  to  the  action  of  externally  applied  forces 

(11) 

can  be  written  as  -  (_FJ  {x}  .  The  components  of  {f}  are  the 

components  of  the  generalized  forces  corresponding  to  each 
generalized  coordinate . 

e.  Forces  which  retard  the  motion  (damping  forces)  are 

assumed  to  be  linearly  related  to  the  generalized  velocities. 

Hence  the  column  vector  of  damping  forces  equals  -[c]  {x}  where 

[Cj  is  a  symmetric  matrix  of  constants. 

Under  the  above  assumptions,  the  equations  of  motion  in 

(12) 

terms  of  the  generalized  coordinates  are 

[m3  {x}  +  [c]  [x]  +  [K]  lx]  =  [F]  (2-1) 

For  a  linear  elastic  structure,  we  have 

IS}  =  [J]  {x}  (2.2) 

where  the  elements  of  {s}  are  the  components  of  the  internal 
forces  within  each  member  of  the  structure.  [jj  is  a  matrix  of 
constants  depending  on  the  elastic  properties  of  the  structural 
members . 
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Once  equations  (2.1)  have  been  written  for  a  given 
structure,  they  may  be  greatly  simplified  according  to  the 
specific  nature  of  the  problem.  We  present  the  following 
example  which  is  of  utmost  importance  to  the  present  study. 


Example  2._1  Consider  the  case  in  which  the  structural  system  is 
not  subjected  to  external  forces  but  L(6  £  L  <  N)  of  its 
generalized  coordinates  are  prescribed  functions  of  time. 

We  can  partition  the  matrices  of  Eq .  (2.1)  to  obtain 


Mrr_j_^rf 
Mfr  1  MffJ 


R 


C  C  * 
__rrj  _r_f 

c  fc 

.  fr  I  ff. 


R 


K_rrJ  “rf 

K  !  K 
.  fr  !  ffJ 


JR 


(2.3) 


where  {/  }  are  those  generalized  displacements  which  are  given 

R 

functions  of  time,  and  [F  }  are  the  reaction  forces  at  the 

R 

locations  where  displacements  are  prescribed. 

In  expanded  form,  Eq .  (2.3)  can  be  written 


[Mrr-  UR)  +  r.HrfMif}  +  [crr](f.R)  ♦  tcrfK»f] 


+  [K  1U  }  +  [K  j  ( A  }  =  (F  } 
rr  r  rf  f  r 


(2.4) 


rV!V  +  +  [Cfr1(;R}  +  [C££K!£1 


+  -Kfr:|;R]  +  -  101 


<2. 5) 
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Now,  we  define  the  relative  displacement  vector,  [u] ; 


Uf]  =  -[Kff]_1[Kfr]  {Ar}  +  {u} 


(2.6) 


Substitution  of  Eq .  (2.6)  into  Eq .  (2.5)  yields 


[M  -Mu]  +  [CffKu]  +  [K  ](u]  =  [P]{A  }  +  [R]{A  } 
ft  rt  rr  r  R 


(2.7) 


where 


(2.8) 


(2.9) 


Substitution  of  Eq .  (2.6)  into  Eq .  (2.4)  gives  expressions 


for  the  reaction  forces  m  terms  of  (u} ,  (A  }  and  their  deriva- 

R 


tives.  Thus,  in  this  case,  (2.7)  is  the  matrix  equation  of 
motion . 

In  the  case  for  which  the  damping  forces  are  small,  the 


(12) 


last  term  on  the  right-hand  side  of  Eq .  (2.7)  may  be  neglected. 


Hence,  for  small  damping,  the  approximate  equations  of  motion 
are 


[Mff]  fu}  +  [Cffj(u}  +  [Kff]  [u]  =  [PlUR] 


(2.10) 


(N-L) x (N-L) 

where  [p]  is  again  given  by  Eq .  (2.8). 

2 . 3  Normal  Equations  of  Motion 

For  the  moment,  we  shall  consider  Eq .  (2.1)  without  the 

damping  term  and  with  the  applied  forces  set  equal  to  zero. 


[mMx]  +  f  K]  lx}  =  (0} 


(2.11) 
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That  is,  we  shall  consider  the  free  vibration  of  an  undamped 
structure . 

Since  (2.11)  is  linear  with  constant  coefficients,  we  may 
try  a  solution  of  the  form 


[x]  =  {a]  exp(ixt)  (2.12) 

uhere  [a]  is  a  constant  vector. 

Substitution  of  (2.12)  into  (2.11)  yields 

-[m]x  [a]  +  [ K] (a }  =  {0}  (2.13) 

after  the  complex  exponential  factor  has  been  divided  out. 

The  condition  to  be  satisfied  for  (2.13)  to  yield  a  non-trivial 


solution  for  the  components  of  {a}  is 

I [k]  -  x2[m]|  =  o 


(2.14) 


This  is  an  algebraic  equation  of  order  N  in  x  ,  called 


the  characteristic  equation . 


(12) 


Let  us  denote  the  solutions 


2  2 

by  x  (i  =  1,2,...,N).  To  each  separate  root  of  (2.14),  jj  , 
i  i 

there  corresponds  a  solution  co  (2.13)  which  we  denote  by  {a.}. 

i 

th 

Now,  let  us  define  the  matrix,  Lcpj  ,  whose  i  column  is 

( NxN) 

composed  of  the  elements  of  {a.}  and  make  the  transformation 

(xl  =  [cpHU  (2.15) 

It  can  be  shown  that  the  application  of  this  transformation 

to  Eq .  (2.13)  with  the  term  [f]  produces  uncoupled  equations 

r  3  (H/12) 

of  motion  in  terms  of  the  components  of 

=  [0] 

(equation  continued) 
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where 


C~m-J  =  [cp1T[M3[cp] 

C~~  K — ^  =  [cp]T[K][cp]  = 

i 

(Ql  =  [cpl  T{  f} 


(2.16) 


To  extend  (2.16)  to  the  case  wherein  damping  is  present, 

we  assume  that  [c]  is  proportional  to  either  [k]  or  [m]. 

Under  this  assumption  the  transformation  (2.15)  also  diagonalizes 

[c].  It  can  then  be  shown  that  the  full  equations  of  motion 
(12) 


are 


Vi)  +  +  \r*2~^{z}  =  ["  m-j  1  {q} 

11  1  (2.17) 

where  [’"m-J  =  [cp]  [M][cp] 

{Ql  -  [cp]T[F],  2["!.x.^  =  l^mO_1[cp]T[c][cF] 

1  1 


The  s \  are  called  the  natural  frequencies  of  the  structural 
system  each  of  which  is  associated  with  a  particular  mode  of 


free  vibration. 


(12) 


The  Z  are  called  the  damping  ratios, 
i 


Eq .  (2.17)  are  the  normal  equations  of  motion  and  the  are 

(11) 

the  normal  coordinates. 


We  shall  assume  hereafter  that  % .  (i=l ,  2 ,  . . . , N)  <<  1.0  since 


l 

this  is  true  for  most  engineering  structures  and  particularly 
for  payload  structures.  When  this  holds,  the  free  motion  of 
the  system  takes  the  form  of  lightly  damped  oscillations.^^ 
Finally,  from  (2.15),  it  can  be  seen  that  the  member 
stresses,  [s]  are  related  to  the  normal  coordinates  by 

fS]  =  (2.18) 
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Example  2.2:  Derive  the  normal  equations  corresponding  to 

Example  2.1. 

We  write  the  equations  of  undamped,  free  vibration: 

[M  j{u)  +  [Kff]  [u]  =  0 
ff  ff 

From  (2.13)  and  (2.14)  the  natural  frequencies  and  the  matrix 
[pj  are  obtained  from 

([Kff]  -  x2[Mff]){a}  =  0 

and  | [K  ]  -  x2[Mf f J |  =  0 

Assuming  that  [c  ]  is  proportional  to  [m^]  or  we 

can  write 

[v]  +  2[^.x.<l{h}  +  [''x  ; — ]  { h}  =  lT~m~-!]-1[cp  JT[  P]  { A  } 

1  1  1  r 

where  {u}  =  [cp]{r,}  (2.19) 

2[~~-;.x.^3  =  1[cpjT[c  ][cp] 

ii  f  f 

[PJ  =  [Mf f 1 [ Kf f [Kfrl  -  [Mfr] 

Since  equations  (2.17)  and  (2.19)  will  form  the  basis  of 
much  later  work,  we  consider  their  solution  in  the  next  section. 


2  -  4  Response  Formulation  to  Deterministic  Forces 

t  h 

Consider  the  K  normal  equation  from  (2.17); 


k 


2:  x  f  +  x  % 

k  k  k  k  k 


Vu 

m 

k 


and  suppose  that  Q  (t)  is  a  unit  delta  function; 

k 

Q  (t)  =  6  ( t  —  t  ) 
k 


(2.20) 


(2.21) 
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We  shall  denote  the  resulting  5  (t)  by  h  (t,T).  h  (t,T) 

k  k  k 

th  (12) 

is  the  impulse  response  function  of  the  k  mode. 


The  solution  to  Eq.  (2.20)  is  then  easily  obtained: 


(2) 


h  (t-T)  = 

k 


“kd^ 

0 


e  K  K  sin  uu  (t-T)  t-T  >  0 


kd 


t-T  <  0 


(2.22) 


where  x 


kd 


=  x  ./l-C2  . 


Now,  any  arbitrary  Q  (t)  may  be  constructed  from  a  series 

k 

of  impulses : 

CD 

Q  (t)  =  f  Q  (T)6(t-T)dT  (2.23) 

k  k 

assuming  that  Q  (t)  =  0  for  t  <  0. 
k 

Hence,  when  the  structure  starts  its  motion  from  rest  at 
t  =  0,  the  response  of  the  k^  normal  displacement  is  ^  ^ 


(t)  =  J  h  (t-T)O  (T)dT 

°  0  k  K 


(2.24) 


with  h  (t-")  given  by  Eq .  (2.22). 

k 

Since  this  equation  holds  for  all  the  normal  coordinates 
of  the  system,  it  represents  a  complete  response  formulation 
to  deterministic  excitations  which  commence  at  t=0. 


Example  2.3:  Obtain  the  response  of  the  system  considered  in 

Example  2.1,  assuming  that  [X  }  =  0  for  t  <  0. 

r 

We  can  write  (2.24)  in  the  matrix  form: 
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{%}  =  J  C^h  (t-T)^3[Q(T)  dT 
0  K 

Now,  [§}  is  replaced  by  {r}  and  [Q(t)}  is  replaced  by 
-1  T 

["mO  [cp]  [P]{A  (t)}.  He  nee,  we  have 

r 

t 

[ri]  =  j  [-h  (t-T)<)[T]{A  (T)}dT  (2.25) 

0  k  r 

where 

[T]  =  ]  -  [M  ])  (2.26) 

ff  ff  fr  fr 

and  the  h  (t-T)  (k=l ,  2,  .  • . ,  N-L)  are  given  by  Eq.  (2.22). 
k 

In  the  next  several  sections  we  consider  the  response  of 
a  structure  to  random  excitations. 

2 . 5  Brief  Discussion  of  Stochastic  Processes 

Probability  theory  will  not  be  discussed  here  but  we 

shall  give  a  heuristic  description  of  a  random  process. 

Suppose  we  have  a  family  of  random  variables,  X^(t^)_ 

(i=l , 2, . . . , n)  which  are  distinguished  by  the  value  of  some 

parameter,  t,  within  a  range  of  values,  (t  , t  ) .  That  is, 

0  n 

there  is  a  one-to-one  correspondence  between  each  random 

variable,  X.,  and  a  distinct  value  of  a  parameter  t,  t=t  .  If 
i  i 

we  conduct  successive  experiments  whose  outcomes  correspond 

to  the  vector  {x  ,X  ,X  ,...,X  }  and  plot  each  of  the  components 
1  2  3  n 

of  the  vector  against  the  corresponding  values  of  t ,  we  obtain 
for  each  experiment  a  function  defined  only  at  the  points  t=t 

i 

(i=l,...,n).  Figure  2  shows  the  outcomes  of  two  such  experiments. 
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When  the  number,  n,  of  random 
variables,  and  of  correspond¬ 
ing  values  of  the  parameter, 
t,  increases  without  limit 
within  the  t  interval  under 
consideration,  we  have  a  con¬ 


tinuously  parametered  stochastic 


process  in  the  parameter  t, 

(4) 

which  we  denote  by  X(t). 


For  our  work,  the  parameter  t 
Fig.  2  will  henceforth  be  associated 

with  time. 

To  completely  describe  a  random  process,  we  must  specify 
something  like 

f  (x  (t  )  ,  X  (t  )  ,  - -  X  (t  )  ) 

x  1  1  2  2  n  n 

where  f  (...)  is  the  n*'*1  order  joint  probability  density  of 

X  *  (2) 
the  random  variables  X_  (t_  X  (t  )  .  Since  this  is 

11  n  n 

usually  impractical,  we  must  be  content  with  more  limited  infor¬ 
mation  such  as  a  determination  of 

f  ( x  ( t ) ) ,  and  f  (x ( t  ),  x(t  ))  for  all  t,  t  ,  and  t  . 
x  x  1  2  12 

or  of  other  first  and  second  order  statistics. 


*  We  denote  a  random  variable  by  a  capital  letter  whereas 
we  denote  the  specific  values  it  may  assume  by  the  small  letter. 
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In  the  next  sections,  the  ensemble  averaging  operation 
is  defined,  and  its  elementary  properties  are  considered. 


2 . 6  The  Ensemble  Average 

The  ensemble  average  of  a  function  of  several  random 

variables  is  defined  as  the  integral  of  the  function  multiplied 

by  the  joint  probability  density  of  the  random  variables  taken 

over  the  entire  state  space.  For  example,  if  G(X  ,X  ,  .  .  .X  ) 

12  n 

is  a  (Borel)  function  of  the  n  random  variables  X^ ,  X^,  and 
X^  then  its  ensemble  average,  denoted  by  E[ G (X^ , X^ / • • • X^ ) ]  is 
(Ref.  2,  p.  23,  Eq.  (2-58)) 


CD  00  CD 

E[g]  =  J  J  .  .  .  J  G(x  ,  x  ,  . 

X_  =-0D  X  X  =  -»  1 

1  2  n 


.x  )F(x  ,x  ,...x  )dx  ...dx 
n  1  2  n  1  n 


(2.27) 


provided  that  the  integral  remains  finite  if  G  is  replaced  by 
| G | .  Here,  f (x^ , , 0 . . , x^ )  is  the  joint  probability  density 

function  of  the  random  variables  X  ,  X  ,...,X  . 

12  n 

Using  the  ensemble  averaging  operation,  we  may  define 
certain  useful  statistics  of  random  processes.  It  is  to  be 
noted  that  for  a  specified  value  of  t,  the  stochastic  process 
X(t)  is  a  random  variable  and  for  a  different  value  of  t  the 
process  is,  in  general,  a  different  random  variable.  Thus,  for 
the  stochastic  process,  X(t),  we  may  define  the  following 
statistical  averages  ^ ^ 
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a.  The  mean  of  X(t)  is  E{x(t)}  and  is  usually  denoted 


by  u  (t) 
x 


The  mean  square  of  X  is  E{x  (t)} 


2 

c.  The  variance  of  X(t)  is  E{x(t)  -  u  (t))  j  =  a 

-  x 

d.  For  n  processes,  X  (t)  .  ..  X  (t),  the  average 

1  n 

K. .  =  E[(X. (t)  -  U  (t)) (X. (t)  -  U  .  (t))] 

i  xi  3 

is  the  second  order  central  moment  of  X. (t)  and  X  (t) .  When 

1  j 

U  .  (t)  -  u  .  (t)  =  0  it  is  simply  the  second  order  moment  of 
xi  x  j 

X^(t)  and  X  (t)  •  In  any  case,  we  shall  most  frequently  use 
the  latter  designation. 

e.  The  autocorrelation  of  X(t)  is 

R  (t  ,t  )  s  E[ X (t  )X(t  )] 
x  1  2  12 

f.  The  autocovariance  of  X(t)  is 

c  (t  ,t  )  =  e[  (X  (t  )  -  u  (t  ))(X(t  )  -  u  (t  ))] 
xl2  1  xl  2  x2 

■  Rx(tl't2)  -  ux,tl)ux,t2) 

g.  The  cross -correlation  between  two  random  processes 

X  ( t )  and  Y (t )  is 

R  (t  ,t  )  =  E [X (t  )  Y (t  )  } 
xy  1  2  12 

h.  The  cross- covariance  between  X(t)  and  Y(t)  is  defined 


as 


xy 


(t,.t  )  =  E[ (X (t  )  -  U  (t  ))(Y(t  )  -  u  (t  ))] 

1  2  1  x  1  2  y  2 

=  R  (t  ,t  )  -  u  (t  )u  (t  ) 
xy  1  2  x  1  2  2 


A  most  important  result  associated  with  definitions  a.  and 

(2,4) 

c.  is  the  Tchebycheff  inequality: 
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(2.28) 


P(|x(t)  -  a  (t)  I  >  ka  (t) }  < 
x  x 


where  p{a]  is  the  probability  of  event  A.  The  existence  of 

this  inequality  makes  the  evaluation  of  the  mean  and  variance 

of  a  random  process  of  the  utmost  importance. 

One  more  important  ensemble  average  will  now  be  intro- 

*  (13) 

duced.  If  we  take  the  Fourier  transform  of  each  of  the 

ensemble  members  of  a  stochastic  process  X(t),  we  obtain  under 
suitable  conditions  another  random  process/  X(j.),  defined  in 

the  frequency  domain.  The  generalized  spectral  density  of 

** 

X (t )  can  be  defined  as 

(xi  ,  X  )  =  E[x(x  )X* **(x  )] 
x  2  1  Z 

-  E[X(x1)]E[X*(x2)]  (2.29) 

This  quantity  is  related  to  the  autocovariance  function 
by  the  formulae 

OD 

ffi  (x  ,x  )  =  "  C  (t  , t  )exp(-i (x  t  -X  t  ) )dt  dt  (2.30) 

x  1  2  ‘  Xs  xx  ■*-  2  1122  12 


- -  1  a  (x  ,x  )exp(i(x  t  -X  t  ))dx  dx  (2.31) 

(2n )  —  X  1  2  X  1  2  *  12 


Similarly,  for  two  random  processes,  X(t)  and  Y(t),  we 
have  the  generalized  cross-spectral  density: 


*  See  also  Chapter  III  for  the  definition. 

**  For  a  full  discussion  of  this  material  see  Ref.  (2), 

section  3.8. 
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(2.32) 


®xy  (‘"‘l'  X2)  =  e^x^1)Y*(x2)]  -  E[X(o)i)]E[Y*(iu2)] 
which  is  related  to  the  cross-covariance  function  by 

CD 

®xy(Va2)  =  I  J  CXy(tl't2)eXP('i(Xltl~X2t2))dt1dt2  (2*33) 

oo 

Cxy(tl,t2)  =  "2-1  2  H  ®xy(u;i'Ji2)eXP(i(UJltl"‘iJ2t2))d,',,ldX’2  (2*34) 

2.7  Theorems  on  Linear,  Time  Invariant  Systems 

Now  we  shall  consider  two  theorems  concerning  linear 
time  invariant  systems,  excited  by  stochastic  time  functions. 

Let  the  input  process,  X(t),  be  related  to  the  response  of  the 
system,  Y(t)  by 

Y (t )  =  L[X(t)]  (2.35) 

where  L  is  an  operator  in  t  which  may  be  algebraic,  differential, 
integral,  or  any  combination  of  these.  The  system  is  linear 
when 

L[aiXi(t)  +  a2X2(t)]  =  a1L[X1(t)]  +  a^x^t)]  (2.36) 

where  X^(t)  and  Y2(t)  are  two  different  input  processes  and 

a  and  a  are  constants.  The  system  is  time -invariant  when 
1  2 

Y (t  +  €)  =  L[ X ( t +  e  )  1  ( 2  o  37  ) 

That  is,  when  the  input  is  translated  in  time  by  e,  the  output 

is  translated  in  the  same  way. 

For  a  system  satisfying  Eqs.  (2.36)  and  (2.37)  the  follow- 

(4) 

ing  result,  within  certain  restrictions,  holds  true: 
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E[  Y  (t )  ]  =  L[E[X(t)]] 


(2.38) 


i.e.,  the  operators  E[  ]  and  L[  ]  commute. 
As  a  result  of  (2.38): 

Cxv(tl't2)  =  \  [C 
XY  ±  z  t  ^  xx  x  2 

c  (t  ,t  )  =  l  rc  (t  ,t  )] 
yy  1  2  t^  xy  1  2 


and  so,  C  (t  ,  t  )  =  L  Tc  (t  ,t  )] 
yy  i  2  tx  xx  i  2 


(2.39) 

(2.40) 

(2.41) 


where  L  and  L  are  the  operator  L  expressed  in  terms  of 

tl  t2 

t  and  t  respectively. 

12 

1hus#  from  (2.41)  the  autocovariance  of  the  output  can 
be  expressed  in  terms  of  the  autocovariance  of  the  input. 

2 . 8  Autocovariance  Response  of  the  Normal  Coordinates 

The  form  of  Equation  (2.24)  shows  that  the  structural 

system,  as  we  have  idealized  it,  satisfies  (2.36)  and  (2.37) 

and  is  therefore  a  linear,  time  invariant  system.  If 

is  a  stochastic  process  then  it  follows  that  ?^(t)  is  stochastic. 

The  response  equations  for  -  (t  )  and  %  (t  )  are 

k  1  i  2 


W  =  'o  VVTi,Qk<Ti)dTi 


■>*"  2 

:.(t  )  =  j  h  (t  -T  )Q.(t  )dT 

3  2  °0  3  2  2  3  2  2 


(2.42) 


(2.43) 


Taking  ensemble  averages  of  both  sides  of  the  above 
equations  we  obtain: 
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(2.44) 


W  =  J0  VYTi)E[VTi)ldTi 


U  .  (t  )  =  r  h.(t  -T  )E[Q.(T  )]dT 

3  2  o  0  2  2  32  2 


(2.45) 


Multiplying  c  (t  )  by  c  (t  )  and  taking  the  ensemble 
'k  1  * j  2 

average  and  using  (2.38),  (2.42)  and  (2.43): 


tl  r-fc2 


E[?  (t  )?.(t  )]  =  R  .  (t  ,t  )  =  J*  r  h  (t  -t  ) 

k  1  3  2  k3  1  2  J  o  o  k  1  1 


-h  (t  —  )E[Q  (-  )Q.  (t  )]dT  d-  (2.46) 

322  k  1  3  2  12 


Multiplying  (t  )  by  u  (t  )  and  employing  (2.44)  and 
k  1  j  2 


(2.45)  : 


^1*2 

U,  (t  )u  .  (t  )  -  I  !,  h  (t  -T  )h  (t  -T) 


k  1  3  2 


0  0 


k  1  1  3  2  2 


•  E[Q  (t  ) jETq  (t  ) ]dT  d’ 

k  1  3  2  12 


(2.47) 


Now,  if  we  subtract  Eq.  (2.47)  from  (2.46)  and  employ 
definitions  g.  and  h.  of  section  2.6  we  see  that 

t  -I  t 

_  r  1  -  2 


C  (t  , t  )  =  *  h  (t  --  )h  (t  — t  ) c  (t  .t  )dr  d- 

k  j  1  2  ^  J  ‘  ' 


0  0  k  1  1  j  2  2  Qkj  1  2  12 


(2.48) 

th 

where  C,  .  (t  ,t  )  is  the  cross  covariance  of  the  k  and  j 
K.  J  1  Z, 

normal  coordinates  and  C  ,r  )  is  the  cross-covariance  of 

th  th  °kj  1  2 

the  k  and  3  components  of  the  excitation. 

Furthermore  from  definition  d.  section  2.6  it  is  apparent 


that 
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-  Kij<t) 


(2.49) 


[Ck 


j<trt2)1 


W* 


Therefore,  from  Eq.  (2.48): 


\i(t)  =  j  j 
KD  0  0 


=  r  '  h  (t-T  )h  (t-T  )c  (T  ,t  )di  dT 
k  1  3  2  1  2  12 


(2.50) 


and 


Now,  from  (2.18)  and  (2.38): 

E[{sn  =  tJKcpMun 

E[(S)LSJ]  =  [jj[cp]E[  U3  L5J  ][cp]T[  j]T 


so  that  from  definition  d.,  section  2.6: 

[It  J  =  [J][cp][tc  (t)][c?]T[j]T 
skj  k3 


(2.51) 


where  [  K  1  is  the  matrix  of  second  order  moments  of  the  com¬ 


ponents  of  the  stress  vector.  The  diagonal  terms  of  [K  (t)] 

Skj 

are  the  variances  of  the  member  loads.  These  can  be  determined, 

according  to  Eq.  (2.51),  only  when  [K  (t)j  is  first  known. 

kD 

Thus,  in  later  sections  we  shall  employ  Eq.  (2.50)  repeatedly. 

In  the  final  section  of  this  chapter  two  important  classes 
of  stochastic  processes  which  will  be  useful  in  later  discussions, 
are  introduced. 


2  o  9  Two  Special  Classes  of  Stochastic  Input  Processes 

Here  we  consider  weakly  stationary  and  shot  noise  processes. 

A  stochastic  process  is  said  to  be  weakly  stationary  when 

its  mean  is  a  constant  and  its  autocorrelation  can  be  expressed 

( 2 ) 

as  a  function  of  i*e*'  x(t)  i-s  weakly  stationary  when' 
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(2.52) 


U  (t)  =  constant 
x 


and 


R  (t  ,  t  )  =  R  (t  -t  ) 
x  1  2  x  2  1 


(2.53) 


where  R  (t  -t  )  is  an  even  function, 
x  2  1 

It  can  be  shown  that  if  the  input  to  a  linear  time  invariant 

(3,4) 

system  is  weakly  stationary,  so  is  the  response  of  the  system. 

Furthermore,  it  can  be  shown  that  the  generalized  spectral 

(4) 


density  of  a  stationary  process  is  of  the  form 

CD  (x  ,  x  )  =  S  (u)  )6  (uu_— uu_  )2rr 
x  1  2  x  1  2  1 


(2.54) 


-IX,  T 


and 


S  (x  )  =  J  R  (t  )  e  1 

X  1  X 


dr 


(2.55) 


where  t  =  | t | 

Finally,  for  a  stationary  excitation,  the  second  order 
moment  between  the  and  normal  displacements  is  given 


by 


(3) 


1  * 

<  .  =  ~  H  U)H  (x)S  (t)dx 

^3  2  la>  k  j  x 


(2.56) 


where  H  (x)  and  H  (x)  are  the  Fourier  transforms  of  h  (t)  and 
k  j  k 

h  (t)  respectively  and  are  called  the  complex  frequency  response 
3 

_  th  ,  . th  , 

functions  of  the  k  and  j  modes. 

Finally,  we  say  that  a  random  process,  S(t),  is  a  short 

^  (2) 

noise  when  its  mean  and  covariance  functions  are  given  by 

U  (t)  =  0  (2.57) 

s 

t2-58) 
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where  l(t)  is  called  the  intensity  function  of  the  shot  noise. 

When  each  of  the  components  of  a  vector  { F]  is  a  shot 

(2) 

noise  process,  we  call  [Fj  a  vector  shot  noise.  In  this 

case : 

[C  (t  ,t  n  =  [I.  .  (t  )]6  (t  -t  )  (2.59) 

sij  1  2  ID  1  2  1 

This  concludes  the  discussion  of  the  probabilistic  theory 
of  structural  dynamics0  It  is  hoped  that  this  brief  intro¬ 
duction  will  suffice  in  providing  the  foundations  of  our  later 
work . 
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CHAPTER  III 


ELEMENTARY  DISCUSSION  OF  THE  FOURIER  INTEGRAL 

3.1  Introduction 

In  this  chapter,  we  define  the  Fourier  integral  and 
describe  its  elementary  properties.  The  applications  of  the 
Fourier  integral  to  analysis  of  system  response  to  both 
deterministic  and  random  excitations  is  also  developed.  The 
chapter  concludes  with  the  derivation  of  an  approximate  inver¬ 
sion  formula  which  is  shown  to  be  closely  related  to  Reimann's 
lemma • 

★ 

3.2  Definitions  and  Simple  Theorems 

We  shall  define  the  Fourier  integral,  F(x),  of  a  function 
of  time,  f (t ) ,  as 

F(jl)  =  J  f(t)  e  1X  dt  (3.1) 

—  00 

where  the  parameter  oj  is  called  the  frequency  and  is  always 
a  real  quantity,.  Note  that  from  (3.1),  F(x)  is,  in  general,  a 
complex  quantity.  F(x)  is  also  called  the  Fourier  transform 
of  f(t).  Henceforth,  we  shall  adhere  to  this  designation. 

It  is  known  that  f(t)  may  be  represented  in  terms  of  its 
Fourier  transform  by  the  following  equation,  known  as  the 
inversion  formula : 

*  The  material  in  this  section  is  taken  from  Ref.  (13). 
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f  (t) 


_1 

2tt 


(3.2) 


I  F<*> 


iojt  , 
e  dua 


The  above  representation  is  valid  only  at  continuity  points 
of  f(t).  However,  the  inversion  formula  (3.2)  also  holds  at 
discontinuity  points  if  it  is  assumed  that 

f(a)  =  i  (f(a+)  +  f (a-) )  (3.3) 

2 

where  t=a  is  a  point  of  discontinuity.  f(a  )  and  f (a  )  denote 
the  value  of  the  function  as  t  approaches  a  from  the  right 
and  from  the  left  respectively. 

Sufficient  conditions  for  the  validity  of  (3.2)  are  given 
in  Ref.  (13),  p.9.  The  functions  we  shall  deal  with  will  be 
assumed  to  satisfy  these  conditions. 

When  a  function  of  time,  f(t),  and  a  function,  F(r),  are 
related  as  in  Eqs.  (3.1)  and  (3.2),  we  shall  indicate  the 
relationship  by 

f(t)  ^  F (uu) 

We  shall  now  consider  the  following  simple  theorems  : 


Theorem  1 :  Linearity 

It  is  apparent  from  the  definitions  that  if 
f  (t)  ♦-»  F1(x)  and  f  (t)  *-*  F2  (x) 

then 

(a1fi(t)  +  a^f2  (t)  )e  litdt  =  a;LFi(x)  +  a2F2(x) 
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or 


a  f  (t)  +  a  f  (t)  ♦-»  a  F  (n)  +  a  F  (x) 
11  22  11  22 


(3.4) 


The  above  can  be  extended  to  finite  sums 

a  f  (t)  +  ...  +  a  f  (t)  a  F  (x)  +  ...  +  a  F  (x) 

11  nn  11  nn 

but  the  extension  to  infinite  sums  is  not  always  valid. 

Thus,  as  we  see  from  (3.4),  the  Fourier  integral  operator 
is  linear. 


Theorem  2:  Time  Shifting 

Let  us  replace  t  in  Eq. 
prescribed  (real)  constant. 


(3.2)  by  t-t  where  t  is  a 
0  0 

Then  d  (t-t  )  =  dt  and 
0 


co  ix(t-t  )  ®  -ixt  .•  4- 

f(t-t  )  =  —  f  F  ( x )  e  0  dx  =  —  r  [  F  (x  )  e  °]e  dx 

o  2r  2TT  i. 

or 

~ixt0 

f(t-t  )  F  ( x )  e  (3.5) 

0 

where  F(x)  is  the  transform  of  f(t). 


Theorem  3 :  Time  Differentiation 

Differentiating  both  sides  of  Eq .  (3.2)  with  respect  to 

time,  we  obtain: 


or, 


df  _  1 

dt  2-  1 

df  (t) 
dt 

Doing  this 


ixt 

[F(x)ixje  dx 

++  ix  F  (x ) 
n  times,  we  have 
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(3.6) 


Theorem  4;  Symmetry 


If  we  rewrite  (3.2)  in  the  form 


oo 


2n  f(-t)  =  J  F(x)e  XXtdx 


and  interchange  t  and  x  we  obtain  a  statement  of  the  sym- 
metrical  character  of  Eqs.  (3ol)  and  (302): 


(3.7) 


F  (t )  2-  f(-x) 


Thus,  theorems  obtained  in  the  time  domain  may  be  easily 
extended  to  the  frequency  domain. 

We  shall  now  state  without  proof,  two  theorems  funda¬ 
mental  to  the  application  of  the  Fourier  transform  to  the 
analysis  of  linear,  time  invariant  systems. 

Theorem  5 :  Time  Convolution 

Let  g(t)  and  f(t)  be  two  time  functions  for  which  the 

inversion  formula  is  valid.  The  integral  j  f (t-7 ) g (t ) dT  = 

00 

f  (')g  (t-r)dt  is  the  convolution  of  f  ( t )  and  g(t)  in  the  time 


domain.  We  shall  use  the  notation 


00 


l) 

—  CO 


f(t—  )g(-)d-  =  f  ( t )  *g  ( t ) 


(3.8) 


The  Fourier  transform  of  f(t)*g(t)  is  given  by 


f(t)*g(t)  ~  F(x)G(x) 


(3.9) 
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where  F(x)  and  G(x)  are  the  Fourier  transforms  of  f(t)  and 
g  (t )  respectively. 


Theorem  6:  Frequency  Convolution 
From  Eqs.  (j.7)  and  (3*9) 

CD 

F(t)g(t)  ♦-*  — ^  f  F(x-y)G(y)dy 

2tt  l» 

or  f(t)g(t)  ♦— ♦  — -  F(oj)*G(x)  (3.10) 

2tt 

where  now,  the  convolution  is  taken  in  the  frequency  domain. 

For  much  of  our  work,  we  shall  need  the  de lta  function , 

6  (t-t^)  .  This  is  not  a  function  in  the  usual  sense,  since  it 
is  not  defined  by  its  values  for  given  values  of  the  argument, 
but  by  the  following  integral  properties 


^  (t-tQ)p(t)dt  =  cp(tQ) 


(3.11) 


where  ^(t)  is  an  arbitrary  function,  continuous  at  the  given 


point  to.  Similarly: 


d  Mt-t  )  n 

- —  cp  (t)dt  =  (-1)  - 


(3.12) 


dt  dt 

The  Fourier  transform  of  the  delta  function  is  readily 
obtained  from  (3.2)  and  (3C11): 


00  -ixt 


F(x)  =  „  e 


P>  (t ) dt  =  1 


or 


6  (t) 


(3.13) 
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Applying  "Theorem  2  to  (3.13)  we  see  that 


-iuut 

6(t-tQ)  *-♦  e 

From  this  result  and  Eq.  (3.7)  we  have 
e  2^  6  (a;-a0) 


(3.14) 


(3.15) 


Since  the  delta  function  appears  often  in  later  discussions 
we  shall  make  frequent  use  of  the  above  results. 

In  the  next  section,  the  Fourier  transform  is  applied  to 
the  determination  of  the  response  of  a  linear  system  to  deter¬ 
ministic  excitations. 


3 . 3  System  Response  to  Deterministic  Excitations 

Consider  the  normal  coordinate  of  a  structural  system. 

Its  response  to  a  deterministic  excitation,  Q  (t),  is  given  by 

K 

Eq.  (2.24): 


„t 

?  (t)  =  ,  h  (t-T)Q  (T)dT  (3.16) 

k  q  k  k 


Since  O  (t )  is  zero  for  t  <  0  and  (t— t )  is  zero  for 
>  t,  the  limits  of  integration  in  (3.16)  can  be  replaced  by 


-or  and  oc: 

?k(t)  =  f  \(t"T)Q  (T)dT 

—  CO 


(3.17) 


Taking  the  Fourier  transform  of  both  sides  of  (3.17)  and 
applying  Eq .  (3.9)  gives 
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(3.18) 


JL,  ( Ju)  =  H  (x)0  (ju) 
k  k  k 

A - /  * 

where  ^  (x),  H  (x)  and  Q  (x)  are  the  transforms  of  ?  ,  h  (t) 

"~7  k  k  k  k 

k 

and  Q  respectively.  We  can  express  (3.18)  more  conveniently 
k 

in  matrix  notation  as  follows 


(uu) 


[  Hk  3  {  ^(X)  } 


(3.19) 


For  the  impulse  response  function  given  by  (2C22)  the 

.  (2) 

Fourier  transform  is 


H  (x)  =  ST  ~2 — T - 

^  k  x  -x  +2iC  x  oj 

k  k  k 


(3.20) 


Eq .  (3.18)  shows  that  once  the  transform  of  the  excitation 

is  obtained,  the  transform  of  the  response  can  be  readily 
calculated  and  application  of  the  inversion  formula  (3.2)  gives 
the  corresponding  time  function. 

3 . 4  Formulation  of  Response  To  Random  Excitations 

Recall  the  definition  of  the  generalized  spectral  density 
of  a  random  process,  X(t): 

®x(jbl,X2)  =  E[X('X1)X*(X2)]  -  E[X(x1)]E[X*(x  )J  (3.21) 


This  can  be  immediately  extended  to  a  vector  process. 


{ X± ( t ) 1  ; 
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[®  (a;^  - oj  )  1  =  E  [  { X .  (w  )}  |.X?(x  )j] 

i]  1  2  l  1  l  2 


-  E[{Xi(u)  )]]E[LX*(U)^  1 


(3.22) 


where  [ffi. .  (x  ,x  )]  is  the  matrix  of  generalized  spectral  densities 
ID  1  2 

of  the  components  of  the  vector  {X. (t)}. 

[©  (x  /X  )]  is  related  to  the  matrix  of  covariance 
ij  1  2 

(2) 

functions,  [  Ch  (t^ ,  t  ^ )  3 


[©.  (x  ,x  )]  =  ('I'  [C,  (t  ,t  )]exp(-i(u)  t  -  x  t  ))dt  dt 

ID  12  1®  ID  1  2  1122  1  ; 


(3.23) 


tc.  (t  ,t  )j  =  - 

x3  1  2  (2tt) 


-  j  J[® .  . (x  ,  iJexp(i  (x  t  +  x  t  )  )dx  dx^ 
2_®13l2  11  22  12 


(3.24) 

With  the  above  definitions  and  the  results  of  section 

3.3,  we  can  relate  the  spectral  density  of  [5  (t)}  to  the 

i 

spectral  density  of  {Q  (t)}. 

k 

Replacing  x  in  Eq.  (3.18)  by  x  ,  we  have 


)]  =  1>H  (x  {6  (X.)5 

1  k  1  K  1 


(3.25) 


Now,  replacing  x  by  taking  the  transpose  and  then  the 


conjugate  of  Eq .  (3.18),  we  have 


(x2) 


'K  (,Vj  ^>2^ 


(3.26) 


From  (3.19)  and  (3C20),  the  quantitites 


E[  (S'  (x.  )  }  Ly  (x„)  j  j  and  E[  {£.  (x  )  }  jE[  I  (x  )  j  ]  are  given 


'*1 


by 
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E[  (w  )H3^*(w  3  =  (UJ  )^EttQ  (w  )}  LQ*(u)  )J  3 

v  ^  i>.  2  k  1  k  1  k  2 


•C^H  (U)  )>0 
k  2 


(3-27) 


e[  [3T  (uj  ))3  e[  (_3T  (id  )j  3 

k *  1  k  2 


(>H  (x  )<Ie[(Q  (x  )}jE[LQ  (x  )j3I>H*(w  )•<! 
k  1  kl  k  2  k  2 


(3.28) 


Subtracting  (3.28)  from  (3.27)  and  employing  (3.22),  we 
obtain 

[ffir(x1,x2)3  =  ['^Hk(x1K][ffiQ(x1,X2)3[^H*(X2)'>.3  (3.29) 

A 

where  [ffi  (x^fA^Jj  is  the  matrix  of  generalized  spectral 
densities  of  the  vector  process  {Q  (t)}« 

K 

It  can  be  seen  that  the  important  result  (3.29)  is  the 

counterpart  of  Eq.  (2.48)  in  the  frequency  domain.  When 

[ffi  (x  ,x  )]  is  known,  [3L  (x  ,  x  )j  may  be  readily  obtained  and 
Q  1  2  -12 

the  correlation  matrix  of  the  response  calculated  by  a  double 
application  of  the  inversion  integral,  Eq.  (3.2) 


3.5  Evaluation  of  the  Inverse  Transform 

It  is  clear  from  the  preceeding  sections  that  efficient 
application  of  the  Fourier  transform  methods  to  both  deter¬ 
ministic  and  random  response  of  a  system  depends  upon  the 
application  of  the  inversion  formula,  Eq .  (3.2).  Therefore, 

we  shall  consider  the  calculation  of  the  integral 

CD 

i  «  ixt 

—t-  F(i)e  dx  when  F(x)  is  a  given  function  of  «u. 
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The  integral  may  be  conveniently  calculated  if  we  allow 
x  to  assume  values  in  the  complex  plane.  That  is,  we  replace 
x  by  the  complex  number  z  =  x  +  iy  and  consider  the  integral 


1  «  lZt 

I  =  —  0  F (z ) e  dz  (3.30) 

2n  "  c 

where  c  is  a  conveniently  chosen  closed  curve  in  the  complex 

plane.  In  this  discussion  we  assume  that  F(z)  is  analytic 

everywhere  in  the  finite  z  plane  except  at  certain  isolated 

points  a  ,  a  ,...,  a  where  it  has  poles  of  various  orders. 
12  n 

When  this  is  true  we  may  use  the  calculus  of  residues  to 
evaluate  the  inversion  integral.  To  do  this,  however,  we  must 

introduce  the  following  essential  theorems,  stated  without 

★ 

proof . 

Theorem  1 : 


Fig.  3(a) 

if  C  is  in  the  first  and/or 
R 

known  as  Jordan's  lemma. 

*  See  reference  (14),  p. 


Suppose  that,  on  a  circular  arc  C 

R 

with  radius  R  and  center  at  the 

origin,  F(z)  0  uniformly  as  R  oc. 

Then : 

„  imz 

lim  e  F(z)dz  =  0  (m>0)  (3.31) 

second  quadrants.  This  result  is 

556. 


36 


Theorem  2 : 


Now  let 
and/or 


C  extend  over  the  third 

R 

fourth  quadrants.  Then 


*  -  imz 

lim  J  e  F(z)dz  =  0  (m>0)  (3.32) 


R"“  CR 

Fig.  3(b) 

With  these  theorems  a  procedure  for  evaluating  the 
inversion  integral  can  be  outlined.  We  shall  consider  separately 
the  two  cases  t  >  0  and  t  <  0. 

t  >  0:  Let  the  path  of  integration  in  Eq.  (3.30)  be  taken 

in  the  counterclockwise  direction 

along  the  closed  curve  consisting 

of  a  line  segment  on  the  real  axis 

and  a  circular  arc,  C  ,  as  shown. 

R 

Eq.  (3.30)  can  then  be  written 


I 


—  'R  F(x)e1Xtdx  +  -L  T  F(z)e1Ztdz 
2^  1R  2n  JCR 


(3.33) 


Taking  the  limit  of  I  as  R-+oc  and  applying  The.  1  with  m 
replaced  by  t,  we  see  that 


1  '*00  ,  v  ixt  „  ,  _  _  .  x 

Lim  I  =  —  j  F (x ) e  dx  (3o34) 

*od  2"  _ao 


which  is  the  desired  inversion  integral.  But,  by  the  calculus 
(14) 

o f  residues : 
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2 1 

lim  I  =  i^Res(F(z)e  ;a  ) 

R— co  k  ^ 


(3.35) 


izt  izt 

where  Res(F(z)e  ;a  )  13  the  residue  of  F(z)e  at  the  point 

K 

of  singularity  and  the  summation  is  taken  over  the  singu¬ 
larities  in  the  entire  upper  half  z  plane.  Equating  the  left 
hand  sides  of  Eqs.  (3.34)  and  (3.35)  we  achieve  the  final 
result  : 


f(t)  =  —  F(x)ela;tdju  =  i  £  Res  (F  (z )  e1Zt ;  a  )  (3.36) 

2rr  L«  k  k 


for  t  >  0  and  for  all  a  such  that  Im(a  )  >  0 

k  K 


t  <  0 ;  Now  let  the  path  of  integration  be  as  shown  at  left 

and  let  t  =  |t|.  Then 


1  rR  /  x  -ixt  r  ,  .  -izt , 

I  =  -  F (x ) e  dx  +  F(z)e  dz 


-R 


Again,  taking  the  limit  as  R  oc, 
and  applying  The.  2,  we  have 


1  00 

Lim  I  - - f  F(u))e~1JJtdaj  (3.37) 

r—  217  loo 

Employing  the  same  reasoning  as  for  t  >  0  we  conclude  that 

(3.38) 


f(t)  =  —  7  F(i.)e1Xtdx  = 

2tt  " 

—  CD 


-  i  S  Res(F(z)e 
k 


izt 


;  V 


for  t  <  0  and  for  all  a  such  that  Im(a  )  <0. 

k  k 

When  f(0V  f(0~)  the  value  of  [f(t)]  is 

t=0 
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(3.39) 


[f(t)3  =  7  (f(0+)  +  f CO”) 

t=0  2 


as  has  been  stated  in  section  3.2. 

In  our  applications,  F(z)  will  be  a  rational  function;  so 

that  all  its  singularities  will  be  poles0  For  this  case,  the 

residues  are  most  simply  calculated.  If  F(z)  has  a  pole  of 

(14) 


order  m  at  z=a,  then 


izt  1 

Res(F(z)e  ;a)  = 


m-1 

d  r  .  . m  izt, 

-  — ~r  {  (z-a)  F  z  e  }  (3.40) 

(m-1)!  dzm~1 


z-a 


izt 


since  e  has  no  singularities  in  the  finite  z-plane. 

To  illustrate  the  above  methods,  we  present  the  following 
example . 


Example  3.1 

Consider  a  simple  oscillator  for  which  the  displacement, 

:  (t) ,  is  given  by 

*2  Q(t) 

§  +  2 cn|  +  n  5  =  (3.4i) 

m 

where  Q(t)  is  a  shot  noise  with  a  constant  intensity  function. 
Calculate  the  correlation  function  and  mean  square  of  ?(t). 

Let  the  constant  intensity  function  be  I;  then  from  (2.58) 

VW  =  1  MW  (3,42) 

Thus  Q  is  a  stationary  process.  With  this  correlation  function 
we  have  from  Eq .  (3.23) 
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00 


aQ(x1<  jJ2)  =  I  JJ  6  )exp(-i  (x  t  -i_t0)  )dt,dt 


2  1 


11  2 


1  2 


=  2n  i  6  (x  -X  ) 
1  2 


(3.43) 


by  Eqs.  (3.11)  and  (3.15). 

Therefore,  from  (3.29): 


(&! ,  x^)  -  2tt  I  6  (’JJ1-a'2)H(Ui1)H 


(3.44) 


where  H(x)  is  the  frequency  response  function  associated  with 

A 

(3.41).  Substituting  the  expression  for  ffi  (u>  ,x  )  into  (3.24) 

§  1  2 

we  have 

GO 

CS(tl't2>  =  Jj  Mxi-x2)H(xi)H*(x2) 


or 


exp(i(x  t  -x  t  )  )dx  dx 
1  1  2  2  12 


,  „  I  “  ,  ,  .  2  i'x  (ti  -t  ) 

C5  (fcl '  t2)  =  2^  1  H  (x  )  |  e  1  2  dx 


(3.45) 


by  (3.11).  Since  the  combination  (t  -t  )  appears  in  the 

1  2 

exponent  in  Eq .  (3.45) #  we  see  that  C-  is  a  function  of 

(t  -t^)  only.  Hence,  the  response  is  weakly  stationary  and 
we  can  write : 


00  -) 

T  ^  |  1XT 

C„ (t)  =  H(x) |  e  dx 


(3.46) 


where  T  =  t  — t 
1  2 
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We  can  now  apply  the  methods  of  this  section  to  the 
evaluation  of  the  integral  in  (3.46). 

From  (3.20),  the  frequency  response  function  is 

i  2  2  —1 

i_.(Q  -x  +  2i£Qx)  and  so: 
in 


|h(jl) 


1/m' 


2  2  2  2  2  2 
(Q  -X'  )  +  4C  Q  U) 


This  can  be  written  as 


|h(x)|  = 


1/m 


(x-a^)  U-a2)  (x+ap  (x+a^ 


where  a  =  iC^  +  Q  ,  a  =  i£Q  -  Q  and  Q  =  fK/l-£- 
1  <3  2  d  <3 


(3.47) 


Substituting  (3.47)  into  (3.46),  we  obtain 

®  i'JUT 


c^(-)  = 


J 


dx 


0_2  „  (x-a  ) (x-a  ) (x+a  ) (x+a  ) 

2  m  -“1212 


(3.48) 


Let  us  first  consider  the  case  t  >  0.  We  first  replace 
x  by  z  =  x  +  iy  and  adopt  the  contour  of  integration  shown  in 
Fig.  3(a).  Employing  (3.36)  and  (3.40),  in  the  expression 
for  C_ (t ) ,  we  have 


-CHt  . 

Cp(T)  =  — e'3~  —  [cos  Q,t  +  - ^ 


—  sin  n  t] 

d  JTF  d 


4£Q  m 

for  t  >  0. 

For  the  case  "  ^  0,  we  adopt  the  contour  of  Fig.  3(b) 
Setting  t  =  |t|  and  employing  (3.38)  and  (3.40): 
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(T  ) 


I  e 


-COT 


4CC1  m 


_  [cos  Q  ,T 
3  2  d 


+  — ;  £  —  sin  fl  3 

d 


Hence,  for  all  7,  the  correlation  function  of  5(t)  is 


given  by 

I  e"CalTl  r 

C  (t)  =  - — —  (cos  O  r  +  -T-  ■ . -  ■-  sin  H  |  t  |  (3.49) 

1  4QCi3m  d  J  1-C  d 


Since  the  mean  of  §(t)  is  zero,  its  mean  square  is 


K 


cr 


E[§  (t)]  = 


4C03m2 


(3.50) 


3  * 6  Approximate  Inversion  of  the  Fourier  Transform 

In  most  of  our  considerations,  we  shall  be  dealing  with 
functions  possessing  Fourier  transforms  which  are  rational 
functions.  Hence,  it  would  appear  that  the  evaluation  of  the 
inverse  Fourier  transforms  in  our  applications  is  perfectly 
simple.  However,  when  the  Fourier  transform  to  be  inverted 
has  many  singularities  in  the  complex  plane,  the  resulting 
expression  for  the  inverse  transform  may  be  quite  complicated. 
In  such  cases,  it  is  more  desirable  to  obtain  useful  approxi¬ 
mate  forms  of  the  transform  such  that  its  inverse  is  alge¬ 
braically  simplified.  Here,  we  present  one  such  approximation, 
valid  for  many  of  the  cases  to  be  considered.  Essentially, 
the  approximations  involve  ignoring  the  presence  of  singulari¬ 
ties  for  which  the  residues  of  the  transform  are  small.  We 
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shall  consider  only  transforms  which  are  rational  functions 

and  which  possess  no  singularities  in  the  lower  half  z -plane. 

As  can  be  seen  from  Eq .  (3.38),  the  latter  restriction  implies 

that  the  inverse  transforms  are  zero  for  negative  t.  Such 

(13) 

functions  are  called  causual  functions  of  time;  and  our 

later  discussions  will  invariably  involve  them. 

With  the  above  considerations  in  mind,  consider  the 
evaluation  of  the  integral 


Kt)  = 


,  .  iuut 

e _  d  x 


(t  >  0) 


2"  f(x)G(x) 

where  F(x)  and  G(x)  are  polynomials  and  J(z)  is  analytic  in 

the  upper  ahlf  z  plane  and  possesses  no  zeros  in  common  with 

F(z)  and  G(z).  We  shall  assume  the  following: 

a.  F(z)  has  zeros  of  order  one  at  the  distinct  points 

a  ,  a  / - . . / a  and 
12  n 


(3.51) 


lai  -  V  -  6i 


(3.52) 

b.  G(z)  has  zeros  of  order  one  at  the  distinct  points 


b  , . . o ,  b  and 
2  m 


lb-  -  b-;  I  >  e2 

i  3  z 


c.  Finally: 


\  -  'ai“b:  '  i  V2 


where  v  ,  \  >  0 

1  ^ 


(3.53) 


(3.54) 


With  these  assumptions,  it  can  be  shown  that 


See  Appendix  I. 
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I(t)  =  { 


2"G(a)  F (ul' ) 


e^di;]  (1+T.)(1+T) 


where  n-l  m 

m  6  Y  _  MAX[ | J(b  ) | ] 
OtT  }  =  O  { - ~-T - - -  ) 


m-1  n 
'2  1 


r,  €  _  "y','  MIN[  |  J(a  )  |  ] 


(3.55) 


(3.56) 


0{T  }  =  O  [  MAX  - 1]  } 

2  g  (a+A  . ) 

3 


and  where 

a  =  a  +  A  #  ,  A  .  <  5 


( 3  “  1#  .  o  .  #  n) 


(3.57) 


(3.58) 


that  is,  a  is  some  point  in  the  vicinity  of  the  a.'s. 

If  the  Fourier  transform  can  be  factored  into  the  form 
of  the  integrand  in  (3.51)  in  such  a  way  that  the  quantities 
given  by  Eqs .  (3.56)  and  (3.57)  are  small,  then 


I(t)  =  S’  J 


J(U) 


iuut , 
e  da; 


(3.59) 


G(a)F(uu) 


with  a  relative  error  of  (T^  +  +  T^T^)* 

A  situation  for  which  (3.59) 
is  valid  occurs  when  the 
zeros  of  F(z)  are  clustered 
within  a  small  region  of  the 
z  plane  (Fig.  6)  far  removed 
from  the  zeros  of  G(z)  which 
Fig.  6  are  themselves  widely  scat¬ 

tered.  As  an  illustration  of  this  situation,  we  present  the 
following  example. 
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Example  3.2: 

a)  Consider  the  integral: 


K (t )  =  — —  r  six) 

2tt 


ixt  , 
Be  dx 


2  2 
(o+ix)  +  p 


(3.60) 


where  a  and  g  are  real  and  positive,  and  where  ffi(x)  = 

and  the  n  poles  of 


-yd 

iiy 

1— 

1 

2-1 

1 

X 

Fig .  7 


we  may  write 


F(ju) 


are 


K  (t )  =  -  fc(x) 

2~  Icd 


F(x) 

confined  between  the  real 

axis  and  a  circular  arc  of 

radius  r,/2  with  center  at  the 

origin  (as  shown  in  Fig.  7). 

2  2 

Factoring  the  term  (a+ix)  +  B  , 


^  ixt  . 
-B  e  dau 


(x-ai+B)  (x-ai-B) 


Thus  the  poles  of 


2  2 
(a+ix)  +  B 


occur  at  x  =  ai  +  3 •  We 


shall  assume  that  |ai  +_  p|  »  r  so  that  we  may  apply  Eqs.  (3.55) 
through  (3.57)  . 

In  this  case: 

m  =  2,  6  =  ^  ,  e2  =  2B 


~  Ja2  +  =2  -2.  ,  v2  ~  va2+  o2  +  2 
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MAX  [ | J(bk) | 3 

Furthermore,  since  J(x)  =  1,  -  =  1.  Hie 

MIN  [  | J(a  )  |] 
k 


point  z  =  a  may  be  conveniently  chosen  as  the  origin,  i.e., 
a  =  0.  Hence,  from  (3.55),  (3.56)  and  (3.57)  we  have 


K  (t )  =  [~  — —  J*  ®(x)  eUtdi]  (1+T  )  (1+T  ) 

2"  a  +p  12 


where 


n-1  2  2  1  7 

O  [T  )  =  O  {-  - -  (a  +  P  )  } 

1  n  p 


since 


1  2  77  t)  _  r 2  72 

v  a  +  P  +  -  ~  v  a 


~  2 


0  {T  }  =  0  { 


2r, 


J  a  +  p2 


+  p  ,  and 

} 


(3.61) 


(3.62) 


(3.63) 


b)  Now,  let  ^(t)  «— ♦  ffi(x).  Then  by  the  convolution 
theorem  and  the  fact  that  e~atSIN  gt  U(t)  - ~ -  : 

(a+ix;  +  g2 

-a(t-T) 

K  ( t )  =  j  cp(t)  e  SIN  0(t-T)dT 

0 

Hence,  from  (3.61): 

^(t-T)  e  '  SIN  g"  d"  =  y(t)  {I+T  )  (1+T  )  (3.64) 

"0  a2+B  1  2 

Using  Eqs.  (3.62)  through  (3.64)  we  can  estimate  the  order 

of  magnitude  of  K(t).  Consider  the  case  n=l .  We  have  for 

small  a; 

/  2  T 

o  It  !  =  o  I  -+p  }  *  l 

l  3 
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«  1 


o  {T  }  =  o  t^11} 

2  0 


Hence,  if  O  {cpu}  A  ,  then 


O  {  '  cp(t-T)  e~a  SIN  s~  dr}  <  o  {— } 

"0  "  0 


(3.65) 


Now  when  n  >  2  ,  the  order  of  T  is  always  «  1  from 
~  1 

(3*75).  Hence  for  n  >  2: 


r  “OLT  _  /*  7V  t 

O  (  i  cp(t-~)  e  SIN  gr  dr]  «  o  {£} 
n  3 


(3.66) 


Since  3  is  large,  we  can  write  the  general  result: 


o  {  J  Cp(t-T)  e  1T  SIN  0-  dr}  =  O  {-}  (3.67) 

‘  0  P 

In  particular,  for  a  =  0; 

O  (  %t  v  (t-T )  SIN  0T  dr  }  =  0  {£}  (3.68) 

0  P 

3.7  Reimann*s  Lemma 

The  result  (3,68)  in  the  example  of  the  last  section  is 
a  special  case  of  a  result  known  as  Reimann 1 s  lemma o  This  can 
be  stated  as  follows: 

If  ( c* )  is  of  bounded  variation  in  the  range  a  to  b,  and 
has  an  upper  bound  of  order  one,  then 


O  [  v(")  cos  pr  d-}  =0  {— } 

‘a  P 

b  , 

0  (  ;  (r )  sin  0t  dr}  =  O  (-] 

"a  3 

when  p  is  large. 


(3.69) 


(3.70) 
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All  functions  under  consideration  will  be  assumed  of 


bounded  variation.  Simply  stated,  this  means  that  they  can 
be  represented  by  a  curve  of  finite  length  in  any  finite 
interval  of  time. 


(13) 


An  estimate  of  the  variation  of  cp(t)  is  given 


by 


I  v  ( t  )  -  c?(t  )  I  <  M  It  -t  I 
2  112  1 


(3.71) 


OD 


where  is  called  the  first  moment  of  the  transform,  ®(ud). 
From  (3.83)  we  see  that  cp(t)  is  of  bounded  variation  when 
©(z)  has  no  singularities  at  infinity.  Ihis  is  a  necessary 
condition  for  the  validity  of  (3.80),  so  that  Reimann's  lemma 
is  seen  to  be  closely  related  to  the  results  of  section  (3.6). 

To  conclude,  we  now  have,  both  in  the  time  and  frequency 
domains,  techniques  for  the  approximate  analysis  of  system 
response  to  random  excitation.  In  the  following  chapters, 
we  examine  the  second  moment  response  of  various  systems  and 
apply  the  techniques  of  this  chapter  to  their  approximate 
evaluation . 


48 


CHAPTER  IV 


FORMULATION  OF  SECOND  MOMENT  RESPONSE 

4.1  I ntroduct ion 

In  this  chapter  we  consider  the  response  of  structural 
systems  to  random  excitations  as  measured  by  the  matrix  of 
second  order  moments  of  the  normal  displacements.  In  section 

4.2  the  Fourier  transform  of  the  second  moments  is  derived  while 
in  later  sections  the  set  of  variances  is  given  a  special 
treatment.  The  chapter  concludes  with  a  consideration  of 
response  to  shot  noise  excitations. 

4 . 2  Second  Moment  Response 

The  responses  of  the  j ^  and  k*”*1  normal  coordinates  of 
a  structure  for  which  the  equation  of  motion  is  (2.17)  are 
given  by 

t 

:  (t)  =  ;  h  (t-T  ) Q  (T  )dT  (4.1) 

*  “o  k  1  k  1  1 

t 

§ .  (t)  =  h  (t-'  )Q.  (t  )d-  (4.2) 

D  *  o  3  2  j  2  2 

where  we  assume  that  {q}  is  a  random  vector  process  the  com- 
ponents  of  which  are  zero  for  negative  time. 

From  definition  d.  of  section  2.6,  the  matrix  of  second 
order  central  moments  of  the  normal  coordinates  is  given  by: 
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(4.3) 


[K.  •  (t)  J  =  E[{§  -  }  1_§  .  -  Uc  j  ] 

ID  k  D  ?  • 

K  D 

=  E[{=  }  L?J  3  -  (uc  1  LU?  J 
K  x  Sj 

From  (4.3)  and  Eq.  (2.50)  we  have 


i 

0  0 


cvs]  “  n  ^\(t-,i)'3[cokj(Ti-T2nc'hj(t-r2>'Jd-1 


T  dT 


2 

(4.4) 


where  [c  , t 2)]  is  the  correlation  matrix  of  the  excitation, 

Qkj  1 

IQ(t)},  and  the  impulse  response  functions  are  given  by 
Eq.  (2.22). 

Now,  we  consider  only  the  component  K  (t)  of  the  matrix 
of  second  moments  and  take  the  Fourier  transform  of  it.  Denot- 

A 

ing  the  transform  of  K  (t)  by  K_  .  (x),  we  have  from  (4.4); 

kj  KJ 


A  ^ oo  t  t  —  ixt 

*  (x)  =  |’  tj  '  h  (t-T  )h  (t-T  )C  (t  ,t JdT  dT  ] e  dt 

k3  i-  Vo  k  1  j  2  Qkj  1  2  1  2 


(4.5) 


C  (t  )  is  zero  for  t  or  *■_  <  0.  Also,  h  (t-T,)  and 

Qk3  1  2  1  2  k  1 

h  (t-~  )  are  zero  for  t  >  t  and  r  >  t  respectively.  Hence 
-j  2  12 

the  limits  of  integration  in  the  double  integral  within  braces 
in  (4.5)  may  be  replaced  by  -oc  and  +  or,  so  that 


k,.(x)  =  ■"  ( ; ;  h  (t-T  ih.  (t-T  )cn  (t.t  )dT  d-  } 


-ixt 


kj 


1  j'  2  Scj  1#  2  1  2 


dt 


Assuming  the  order  of  integration  may  be  reversed,  we  have 


-l'JQt 


-  u  l  1>hj't-2»e'"''dt>cQk3«-i--2>d-1dT2 


(4.6) 
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Thus,  we  need  only  take  the  transform  of  the  product  h  (t-T  ) 

k  1 

h  (t-T  )  .  From  (3 « 5)  : 

j  2 

-iuuT 

h^tt-T^)  li^(x)e  ,  h  .  (t-T  )  «-♦  H  .  ( x )  e  11  2 


and  from  the  frequency  convolution  theorem: 


V'-VV'-V  *-  ^  j.'Vly) 


-iyT,  -i'2 (x-y) 

e  H  (i-y)e  dy 


Substitution  of  Eq .  (4.7)  into  (4.6)  yields 

\i(x)  =  2^  ^  [  jp  Vy)Hi(*_y)e  1  dy} 

-J  —  CD  —  CD 


(4.7) 


\'VT2ldVT2 


(4.8) 


Again  reversing  the  order  of  integration,  we  have 


‘*3u)  i3(y,Hi(x-y)  vw 


-i(-  y-T  (y-x)) 

i  ZdT  d~  } dy 

12 


(4.9) 


But,  from  (2«33),  the  term  in  braces  is  nothing  but  the 


th  th 

generalized  cross-spectral  density  of  the  k  and  j  excita¬ 
tion  components,  ffi  ,  expressed  in  terms  of  y  and  y-x .  That 


IS 


®  -i(-  y--  (y-x)) 

CQ  <VVe  d"id"2  =  (y,y~x)  (4-10) 

-®  1  z  1  z  uk] 


Hence : 


kj 


U)  =  —  H  (y)H  (x-y)  ffi  (y,y-x)dy 


2~ 


k  j 


(4.11) 


kj 
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or  in  matrix  notation 


^  00 

[K  (x)3  =  —  j*  I>H  (y)-0[i  (y,  y-x)  11>H .  (x-y)<]dy  (4.12) 

KD  AT1  _co  k  q  j 


^  A 

With  [ffi  (x  ,  x  )3  given,  we  can  calculate  [K  (x)3  and 
Q  1  2  kj 

apply  the  approximation  procedure  of  section  3.6  to  obtain  a 
simplified  representation  of  [ K  . (t)l.  However,  in  the  next 
section,  we  derive  a  special  result  for  the  variances  (i.e. 
the  diagonal  components  of  [K  (t)])  for  which  approximations 

kj 

may  be  obtained  directly,  without  an  integration  of  the  type 
considered  in  (4.12). 


4 . 3  Special  Treatment  of  the  Variances 

The  equation  of  motion  for  the  k^  normal  coordinate  is, 
from  Eq .  (2.17): 

2  i 

+  +  a  =  —  0,  (4.13) 

k  k  k  k  k^n^k 


where  r  (0)  =  5  (0)  =  0. 

k  k 


Taking  ensemble  averages  of  both  sides  of  (4.13)  gives 

(4.14) 


•*  2  i 

u  +  2C  x  p  +  X  u  =  —  e[o  3 

k  kkk  kkrru.  K 


% 

Subtracting  (4.14)  from  (4.13)  and  letting  (§  )  =  y 

k  k 

and  (Q  -  E[Q^3)  =  x,  we  have 


4  +  2C  x  +  x2y  =  -A  x 

dt  K  k  dt  k  m 

k 


This  equation  may  also  be  written  in  the  form 
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X 


(4.15) 


dyi  dy2  ,  2  i 

-  =  y  ,  -  =  -2C  my  -my  +  — 

dt  2  dt  k  k  2  ^  1 

dy 

where  y  =  y  and  y  =  . 

1  2  dt 

If  the  equations  of  (4.15)  are  multiplied  by  y  and  y 

1  2 

in  turn,  we  obtain: 

2 

1  dyi 

- =  y  y 

2  dt  1  2 

-  2 

1  y2  2  2  i 

- =  -2Q  x  y  -  x  y  y  +  -  xy 

2  dt  k  k  2  k  1  2  m  2 

k 

dYl  2 

y  -  =  y 

2  dt  2 


dy2  2  2  1 
y,  — —  =  -2C  <*  y  y  -  x  y  +  —  xy 


i  dt 


k  k  1  2  k  1 


m 


k  i 

Adding  the  last  two  equations  and  taking  ensemble  averages 
of  the  resulting  equations: 
dS, 


11“  =  s 

2  dt  12 


(4.16) 


1  dS22 
2  dt 

dS12 


=  ~  2'  x  S._  -  x  S  +  E[ x  y  ] 


k  k  22 


*  12  \ 


(4.17) 


J  -2CxS  -ns  +—  E[x  y  ]  (4.18) 

dt  22  "kk  12  11  *i 


where : 


S..  =  E[y  ] 
11  1 


S  =  E[y  y  ] 
12  *1*2 
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2 

s  =  e[y  ] 
22  2 


Substitution  of  S  as  given  by  (4.16)  into  (4.17)  and 

12 


(4.18)  gives 


1  dS22 
2  dt 


2 

id  dS 

2  c  X  s  - 
k  k  22  2  dt 


11 


+  e[x  y23 


1 

2 


2 

d  S 


dt 


H  r 

2 —  =  s_  -  C,  u>. 


dS. 


22 


11  2  lr  i 

- 'X  S  +  —  E[xyJ 

k  k  dt  k  11  %  1 


(4.19) 


(4.20) 


From  (4.19)  and  (4.20)  we  may  solve  for  S  and  substitute 
this  into  the  remaining  equation  to  obtain  a  single  equation 


in  S  .  Recognizing  that  S  =  K  (t)  we  may  write  this  as: 
11  11  kk 

d  K  d  K  t  2  2  d  K, 

1  kk  3,  kk  ,  *■ 

+  —  l.  x.  - —  +  (x  +  2Q  x  ) 


dt' 


2  k  k 


dt" 


k  k  dt 


kk  1 

+  2C  X  K 


k  k  kk 


=  itr  E[xy^3  +  2CkXkE[xyil  +  Efxy^l  3 


(4.21) 


which  can  be  factored  as  follows 

'2  H  2 

2  +  S.xk  IT  +  11 
dt 


(—  +  2C  U) )  (7  — ,  +  C  x.  +  x  )  K 
dt  k  k  4  .2  k  k  dt  k  kk 


“  ^  +  2c.  at  )  E[xy]  +  ~  E[xy3 

\  2  dt  k  k  \ 


(4.22) 


Now,  from  (2.24)  we  have 


Y  =  (W  =  J  \«t-)x(r)d- 


j  h  (t-")x(~)d- 
* 


(4.23) 


and 
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/  =  J  [—  h  (t-T  )  lx(T  )dT  =  J  [ -JL  h  (t-T  )  ]x(T)dT 


L.  St  k 


(4.24) 


Multiplying  Eqs.  (4.23)  and  (4.24)  by  x(t),  taking 

ensemble  averages;  and  noting  that  E[x(t)x(r)]  is  Cn  (t,T): 

Tck 


E[xy]  =  [  h  (t-T)Cn  (t,-)dT 

:»  k  °kk 


(4.25) 


E[xy]  =  J  h  (t-T)C  ( t , T ) d T 


(4.26) 


(4.27) 
h,  (t-')Cn  (t,T 


Hence,  Eq .  (4.22)  becomes 

(—  +  2W  (7  A-  +  s,  *v  —  +  *2)  •* 

dt  kk  4  ^^.2  k  k  at  kk 

1  1  d  ->c° 

=  —  (-  —  +  2"  x  )  J  K.  (t-T  )  c  (t,r)dT  - 

\  2  dt  k  k  la\  Qkk  \ 

We  shall  now  take  the  Fourier  Transform  of  both  sides  of 
Eq.  (4.27 ) . 

From  (3.6),  we  have 

,2  ,  2 

1  Kvv(t) 

(4.28) 


(—  +  2C  x,  )  (i  SI-  +  C.  x  —  +  x.  )  K  (t) 
dt  k  k  4  ^  k  dt  ^  kk 


\  2  2  * 

(ix  +  2;kxk)  (-  —  x  +  iC^x  x  +  x^)  (x) 


Since  h^(t-T)  H^(x)e 


—  IT  X 


Erxy]  -*  -i.  ;  l  H^_  (y ) 


*  *e-^ye-it(^-Y)  (t,.)dtdTl  dy 


_  CO  —CO 


~Q. 


kk 


or 


)dT 
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1 

E[xy]  —  j  (y)  ffi  (y,y-x)dy  (4.29) 

by  (4.25)  and  the  frequency  convolution  theorem. 

Employing  (3.6),  we  have: 

E[xy]  *-*  — -  ix  r°°H-  (y)  »  (y,y-x)dy  (4.30) 

2n  k  Qkk 

and 

E[xy]  *-*  -r -  f  iy  H.  (y)  06  (y.y-a)dy  (4.31) 

2tt  -»  *  ukk 

Therefore,  from  (4.28)  through  (4.31)  the  Fourier  trans¬ 
form  of  K  (t)  is 

1  a  ijj  ** 

-  TF  v_(^  ~2^vU)vi+y^Hk:  (y»y-u;)dy 

Kkk(i)  =  - - - - - -  (4.32) 

2  2 

m,  (x-2£  x,  i)  (-  —  +  i£  x  x  +  uu 

k  k  k  4  k  k  k 

4 . 4  Summary  of  the  Response  Formulation 

In  general,  the  second  moment  response  to  random  excita¬ 
tions  is  given  by  Eq .  (4.12) 

[*kj(x)3  =  2^”  0*  £  Hj^  (y)  ][ffiQ(y,y-x) j[  (x-y)  ]dy  (4.12) 

a 

and  an  alternative  form  for  the  diagonal  elements  of  [K  .  ] 

kg 

(i.e.  the  variances  of  the  normal  coordinates)  is  given  by 
Eg.  (4.32) 
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2r 


(4.32) 


Th esc  two  equations  summarize  the  response  formulation  of 
our  work  and  will  provide  the  basis  of  all  subsequent  considera¬ 
tions  . 

A  cautionary  note  should  be  given  on  the  evaluation  of 
the  integrals  appearing  in  (4.12)  and  (4032).  Let  us  denote 
either  of  the  integrands  by  J ( x , y )  and  the  integrals  themselves 
by  I  (x)  .  Then  we  may  write 


I('A)  =  [—  J°°  J(x,y)elytdy] 


2  TT  t=0 


which,  by  (3.3)  becomes 


(4.33) 


From  the  discussions  of  section  3.5  it  is  clear  that  the  first 
integral  within  the  brackets  is  to  be  evaluated  from  the 
residues  of  the  integrand  in  the  upper  half  z-plane.  Similarly/ 
the  second  integral  is  to  be  evaluated  from  the  residues  in 
the  lower  half  plane.  Hence/  I(x)  is  one  half  the  sum  of  the 
results  obtained  when  the  integral  is  evaluated  from  the 
residues  in  the  upper  and  the  lower  halves  of  the  z  plane.  This 
we  represent  symbolically  by 
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I(UJ)  =  i-(I  +  I  ) 

2  u  i 

A  consideration  of  Eq.  (4.12)  for  the  case  j=k  shows  that 
it  does  not  reduce  to  the  corresponding  result  in  Eq .  (4.32) 

in  any  obvious  way.  Consequently,  in  Appendix  II,  we  present 
a  demonstration  that  the  two  expressions  for  the  variances 
are,  in  fact,  equivalent. 

In  the  next  sections,  we  consider  the  cases  for  which 
the  excitation  is  (i)  a  stationary  vector  process  and 

(ii)  a  vector  shot  noise. 


4 . 5  The  Stationary  Case 

Analogous  to  Eq.  (2.52)  and  (2.53)  written  for  the  case 
of  a  single  random  variable,  we  define  a  stationary  vector 
process,  [X^(t)}  as  one  for  which 


[c  (t  ,t  )]  =  E[[x-  a  }  LX.-  u  J  1  =  [C 
\3  1  2  *  Xk  3  Xj 

where  (t^-t^)  =  T*  matr^x  generalized 

densities  is  given  by 


„  (-)j 

spectral 


[®  (x  ,x_)]  =  2- 
xk3  1  2 


6  (x  -X  )[S 
1  2  *kj 


(iu  )3 


(4.34) 


(4.35) 


where 

OD  *“i  JJ-  T 

fSv  (0  3  =  J  tc  <')1e  dT  (4.36) 

1  I-  *kj 

If  [Q  (t)}  is  such  a  process  we  conclude  that,  quite 
apart  from  any  considerations  involving  Eqs.  (4.12)  and  (4.32) 
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li^  *  iunT 

[ce  (T  )  3  =  —  j  !>H  (X ) — 3  [  S  (x)3!>H  (x)^Se  div  (4.37) 

'kj  2n  —  k  akj  j 

from  Eq.  (3.29)  and  an  application  of  (3.15).  Hence: 


IX.]  =  ([Cr  (")])  =  —  r  ISh.  (xjSlts^  (x)]NH*(uu)S3d 


t  =  0  2n  I 


JU 


Now,  let  us  examine  Eq.  (4.12).  From  (4.35): 

[a  (y,y-x)]  =  2tt  6  (x)[sQ  (y)] 

°kj  k  j 

Substituting  this  into  (4.12)  gives 


(4.38) 


(4.39) 


[^kj  =  j  [  ^(y)  J[Sq^_]  6(x)[  H_.  (x-y)  Jdy 


kj 


so  that 


[Kkn]  =  —  j  {  j  N(y)N[s  (y)K'-H  (x-y)S 

2"  —  CO  _  CD  ^  j  3 


-CD  _  CD 


ioit  , 

*  6  (x)  e  dx  j  dy 


Applying  (3.15)  and  noting  that  H  (-y)  =  H  (y) ,  we  have 

3  j 

X  ^  =  dr  u  l>H,  (yK![s0  (y)]I>H  (y)X!dy 

k3  2  *  Scj  j 

which  is  in  complete  agreement  with  (4.38) 

Substituting  (4.35)  into  (4.32): 


.  Mx):m(|  -2'^i  t  y)Hk(y)SQkk(y)dy 

K  (x)  -  ~ 2  t 

m  (x-2"  x.  i)  (-  —  +  iC  x.  x  +  x 
k  "k  k  4  k  k  k 


Taking  the  inverse  transform  and  employing  (3.15),  this 
becomes 


5S 


1C  (t)  =  - 

kk  2nmk  -»  2Cku;ki 


»  2ck’^i-y 

J  -  -  -  Hk(y)Sn  (y)dy 


(4.40) 


kk 


If  both  numerator  and  denominator  of  the  integrand  is 
2  2 

multiplied  by  (ju  -  y  -  2iC  'X  y)  and  using  Eq.  (3.20)  for 
k  k  k 


H  (y)  : 
k 


+  2)  (\-  y2~  2ick\y) 

od  uo  2^  i  k  *  * 

=  —  f  — - - -  s  (y)  dy 

OT*  "  y 


kk  2^ 


-°°  2  2  2  2  222  kk 

m  (  (x  -y  )  +  4C  tu  y  ) 

k  k  k  k 


2222  222  2 
Noting  that  m  ( (ui  -y  )  +  4£  uo  y  )  is  |H  (y)  |  and  expanding 

k  k  k  k  k 

the  numerator  of  the  integrand,  we  have 


k  =  -i  ;  [i  +  -ix 

kk  2*t 


2iC  y 


2£  x  x  2C  it 

k  k  k  k  k 


-  ■  -iy  =•]  |h  (y)  |  (y)dy 


3  1  k 


kk 


(4.41) 


Finally ,  since  Sn  (y)|H  (y ) |  is  an  even  function  of  y, 

wkk  k 

3 

the  terms  involving  y  and  y  in  (4041)  vanish. 


Hence,  from  (4.41): 

***  *  ^  c,  lHx(y)l\kiy)dy 


which  agrees  with  the  results  corresponding  to  the  diagonal 
terms  of  both  sides  of  Eq .  (4.38). 

In  later  sections,  we  shall  repeatedly  use  the  result 
(4.38)  to  check  results  obtained  for  non-stationary  excita¬ 
tions  via  (4ol2)  or  (4.32). 


60 


4.6  Shot  Noise  Excitations 


If  [Q^ (t) }  is  a  vector  shot  noise,  then  according  to 
(2.59)  : 

[c  1  =  [I  (t  )]  6(t  -t  ) 

Qkj  jk  1  2  1 

So  that,  from  (3.23): 


=  -  J'-Iik(tl)nj6  (t2_t1)eXp(~i(UJltl~X2t2))dtldt; 


jk 


;  ci  k(tn  a 

_  03 


-it (x  -x  ) 
1  2 


dt 


(4.42) 


Or 


[a.  (y.y-x)l  =  [i_,  (x)3 


(4.43) 


where  I  (x)  is  the  Fourier  transform  of  the  intensity  function, 

i  v(t). 

lk 

Substitution  of  this  into  (4.32)  gives 


[K 


kj 


(x)j  =  ^7  (y)^j[hk  (x)  ][\H^  (x-y)-<]dy  (4.44) 


Without  further  information  about  the  intensity  functions 
we  cannot  carry  out  the  integration  of  (4.44).  Thus,  (4.44) 
represents  the  most  general  result  that  can  be  obtained  via 
(4.12). 

However,  we  can  carry  the  determination  of  the  variances 
of  the  response  further  than  (4.44)  by  employing  (4.32). 

Substitution  of  (4.43)  into  (4.32)  gives 
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(4.45) 


57  '  2CkV+y)Vy>dy 


\(x'2Ck\i)(-  r +  iCkV+V 

The  integral  on  the  left  hand  side  is  easily  evaluated 
when  we  realize  that 

f  Hk(y)dy  =  2r  \(y)e+;LYtdy)t=0 


J  yH  (y)dy  =  -2"i  {-f  J  It  (y)elytdy]  } 

loo  k  dt  2rr  ^  K  t=0 

from  the  definitions  of  section  3.2.  Employing  (2.22),  (3.20) 

and  the  above  formulae,  we  have 

1  ~^VUJkt 

H  (y ) dy  =  2tt  -  e  Sin  ju  t  =0  (4.46) 

1®  k  ‘"kd1")?  kd  t-0+ 


O  ■,  -CV«L  -C  x  t 

yH  (y)dy  =  -2"i  {  -  [  -  e  Sin  iu  t 

k  2  id  m  kd 

kd  k 


1 

+  “ —  e  cos  x  ,t]  +  ] 

kd  t_0 


m 


(4-47) 


lm 

k 

where  we  have  noted  that  t=0  is  a  discontinuity  point  of  the 
derivative  of  the  impulse  response  function .  Finally,  sub¬ 
stitution  of  (4.46)  and  (4.47)  into  (4.45)  yields 
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(4.48) 


K  (x) 
kk 


i  W”' 


2  2  2 
m  (x-2C  x  i)  -iC  x,  x-x  ) 
k  k  k  4  k  k  x 


W"  can  now  take  the  inverse  Fourier  transforms  of 
both  sides  of  (4.48)  to  obtain  the  variance  of  the  kt^1  normal 


coordinate  in  terms  of  a  convolution  integral  between  I  (t) 

K  K 


and  the  inverse  of 


integral 

1  p° 

2^ 


2 - ~2 - 2 

2m  (x-2C  x-  i)  (^x_  -  i£  x  x-x  ) 
k  k  k  4  k  k  k 


ixt 

i  e  dx 


.  Hence  the 


2m2(x-2C  x  i) (i—  -  iC  x  x-x2) 
k  k  k  4  k  k  k 


must  be  evaluated.  Upon  replacing  x  by  z  =  x^-iy  and  employing 

Eqsc  (3o36)  and  (3.40),  the  result  is  found  to  be 
~2£  x  t 


k  k 


2m2  x2 
k  kd 


(1  -  cos  2x  t)U(t) 
kd 


2  7"2  T 

2m  (x-2'  x  i )  (—  -2;  xl  x-x  ) 

k  k  k  4  k  K  k 


(4.49) 


where  U(t)  is  the  unit  step  function  and  x  =  x  N 1-' 

kd  k  “k 


Thus,  using  the  time  convolution  theorem  and  assuming 

that  I  (t)  is  zero  for  negative  time,  we  have  the  result 
kk 


1  ,.t 

*kk(t)  "  "22  .  Ikk(t"')  6 

2m  x  0  *k 

k  kd 


-2'  x  ' 

K  K  (1-cos  2x  t )d" 
kd 


(4.50) 


which  is  the  general  expression  for  the  variances  of  the 
response  to  shot  noise  excitations.  We  conclude  this  chapter 
with  an  illustration  of  the  application  of  Eq .  (4.50). 
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Example  4,1: 


Consider  the  variance  of  the  response  of  a  one  dimensional 
oscillator  to  a  shot  noise  excitation,  for  which  the  intensity 
function  is  a  unit  step.  If  its  equation  of  motion  in  terms 
of  the  displacement,  §,  is 

••  2  l 

§  +  2£Q|  +  o  §  =  -  Q 

m 


then  the  mean  of  §  is  zero,  and  from  (4.50): 


_  2  i  ^t  -2en(t-T) 

K(t)  =  E[§  (t)]  =  — — 2  J  I  (t ) e  (1-cos  2w  (t-T))dT 

2m  Q  0 
d 


i  2 

where  0  =  0  J  l-£  .  Since 

d 


I(t)  = 


I  T  >  0 

_  0  T  <  0 

where  I  is  some  constant,  we  have 


(t)  = 


-2;nt  ,,t  2CflT 


2  2 
2rn  Q 


(1-  COS  2f)  (t-T))dT 
d 


which,  after  a  direct  integration  by  elementary  methods, 
becomes : 


K(t)  = 


2  T 

4m  Cfl 


[1  -  e  t(l+  Sin  2Q^t  +  2  Sin  fit)] 


This  result  is  in  complete  agreement  with  that  obtained 
in  Ref.  (6).  When  t  tends  to  infinity,  K(t)  tends  to  the 
value  obtained  in  the  Example  of  section  3.5: 
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Lim  K(t)  =  — - - ^ 

t-*«°  4m  CH 


so  that  for  large  t,  K (t)  approaches  its  stationary  value. 
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CHAPTER  V 


APPROXIMATE  FORMULATION  OF  RESPONSE 

5 . 1  Introduction 

In  this  chapter,  we  use  the  results  of  sections  3.6  and 
3.7  to  obtain  approximations  to  the  second  moment  response  of 
a  structural  system  to  random  excitations.  This  is  accomplished 
by  a  consideration  of  Eqs.  (4.12)  and  (4.32)  for  a  particular 
class  of  input  processes,  described  in  the  next  section.  A 
special  treatment  of  the  variances  is  given  in  section  5.3 
since  Eq.  (4.12)  is  amenable  to  direct  approximation.  In  the 
remaining  sections  we  treat  cases  involving  both  shot  noise, 
and  filtered  shot  noise  excitations. 


5 . 2  Slowly  Varying  Excitation  Processes 

Recall  the  expression  for  the  generalized  spectral  density 

matrix  of  [q]  in  terms  of  the  correlation  matrix,  [ Cq  ]. 

ij 


Replacing  by  y  and  x^  by  y-ut): 


[*0  (y.y-x)]  =  jJ[Cq  (t  ,t  )lexp(-i(t  y-(y-uu)t  ))dt  dt 

wij  -oo  ij  1  z  1  2  12 


(5.1) 


Now,  consider  the  transformation 

t  =  t  -t  and  t*  =  t 
12  2 


(5.2) 
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Since  the  Jacobian  of  the  transformation  is  one,  Eq . 


(5.1)  becomes 

00  .  .  ★ 

[ $  (Y/Y-'A)]  =  JJ  [cQ  (T,t*)]  e  iy  e  1Xt  dt*dT  (5.3) 

Uij  -«  ij 

We  can  also  write  this  as 


[*Q  (y.y-u)]  =  j  tr_  (t,iu)] 

i  j 


-lyr  , 

e  dT 


13 


(5.4) 


where 


[r,  .(T,i)]  =  j  [CQ  (T,t*)]  e-ixt  dt* 
13  _oo  ij 


(5.5) 


From  Eq.  (3.83),  we  see  that  an  estimate  of  the  variation 


of  C  (~,t  )  with  respect  to  t  is  given  by 

Qij 


:Q  (T*t2^ 

ij 


c-  (T,t*) |  <  M. ■  |t*  -  t*| 


'Q.  . 
13 


13 


(5.6) 


MiD  =  y:  J  M  \t  (t,®) |d® 

Hence,  the  variation  of  the  correlation  functions  with 

* 

respect  to  t  depends  on  the  behavior  of  the  generalized 
spectral  densities  with  respect  to  x. 

We  designate  as  slowly  varying  those  non-stationary 

processes,  [Q.  (t)},  whose  correlation  functions  vary  slowly 

★ 

(in  a  sense  to  be  specified  later)  with  respect  to  t  when 
the  transformation  (5.2)  is  effected.  This  restriction  implies 
that  the  M  (i,j  =  1 , N)  of  (5.6)  are  small  for  all  t.  It  can 

ij 

be  seen  that  the  M  are  small  when  most  of  the  singularities 

ij 

of  T  (^/Z)  are  near  the  origin  of  the  complex  plane.  For 

ij 
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convenience  in  later  discussions  it  will  be  assumed  that  all 

the  singularities  of  T  (t,z)  lie  within  some  small  distance 

i  3 

to  the  origin. 

When  the  damping  coefficients  of  the  system  are  small, 

_2 

(C  =  0(10  ))  then  the  methods  of  sections  3.6  and  3.7  may  be 

used  to  approximate  the  second  moments  of  the  response.  In 
general,  this  can  be  done  through  the  use  of  Eq .  (4.12)  only 

when  the  integration  in  y  is  first  performed.  However, 
approximations  to  the  variances  can  be  obtained  directly  from 
Eq.  (4.32)  and  in  a  general  forme  Thus,  we  give  the  set  of 
variances  a  special  treatment  in  the  next  sections. 

5 . 3  Approximate  Evaluation  of  the  Variances 

In  the  following  discussion  we  assume  that 

a.  The  damping  coefficients  of  the  structure,  C^(k=l,N) 

_2 

are  small  (of  the  opder  of  10  ) 

b.  The  natural  frequencies,  w  ,  are  large  (>  10^). 

k 

c.  The  excitation,  [^(t)},  are  slowly  varying  as 

described  in  the  last  section  in  the  sense  that  all  the 

singularities  of  7  ('#z)  li-e  with- 

kk 

in  the  shaded  region  shown  in 
Fig .  8 . 

d.  The  functions  T.  .  (t,z) 

kk 

(k=l , . . . , N)  can  be  expressed  in 
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the  form 


kk 


(t,z)  = 


Jk(T,z) 

Fk(T'z> 


where  J  is  analytic  in  the  finite  z  plane  and  has  p  poles 

*  Fk 

at  a  ,  .  a  in  the  upper  half  z  plane. 

1  P 

Having  made  these  assumptions,  we  now  consider  Eq .  (4.32) 

.  .  uu 

Recognizing  the  term  (—  -  2b  jj  i  +  y)  in  the  numerator  is 

2  1c  k 

i-(x  -  2Q  x  i)  +  (y  -  £  w  i),  (4.32)  becomes 

2  k  k  k  k 


*  2^~  [  J  Hk(y) 


*kk  (y/y-uy) 

2D 


dy  +  J  (y_ckUjki) 


$  (y.y-o;) 

H  (y)  -  dy  1 

k  D 


(5.7) 


where  D,  =  (-  —  +  iC  x,  x  +  x  )  ,  D  =  (x  -  2C,  U),  i)D, 

1  4  kk  k  2  kkl 


And,  employing  (5.4)  and  (5.5),  we  have 

(t,uu) 


*kk(*>  -  i^r  ;  I  t;"vy)e  dy]  kk 

k  _00  —00 


2d  (uj) 
1 


+  [  (y-s,  X  i )  H  (y)  e  dy  j 

k  k  k 


T,  ,  (t  ,  oj) 
-lyT^  ^  kk  x 


D  (ou) 
2 


}  dT  (5.8) 


1  ^0  1  -iyT  rkk(T'^ 

‘  ^  d  Hk(y)e  dy> 


\  1 


2D  (uj) 


,  oo  •  r  k  (T,uu) 

+  [—  .  (iy  +  “kxk)Hk(y)e  Y  dyj  -  3  dT 


2r  i« 


i  D  (x) 
2 


(5.9) 
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since  (y)  is  the  transform  of  a  causual  function. 

We  shall  now  estimate  the  relative  order  of  magnitude  of 
the  two  terms  in  (5.9). 

00 

First,  notice  that  the  integral  — i-  J  H  (y)e  dy  (t  <  0) 

2tt  -co  k 

is  h  (-t),  so  that  by  Eq .  (3.6),  —  (y)iy  e  1“^Tdy  is 

k  2rr  ^  * 

-  —  (h  (-t)).  Hence,  from  (2.22)  and  assumptions  a.  and  b.: 

dT  k 

1  -iyr  . 

0  [  (iy  -  C,  w  )  H  (y)  e  dy  } 

2TT  ^.00  K  k  k 


=  \  0  1^7  J  Hk(y)e  lyTdy  } 


(5.10) 


And  since  h  ( —  t )  and  -  —  (h  (-t))  are  both  exponentially 
k  dr  k 

decaying  sinusoids  with  the  same  frequency  and  exponential 
factor  : 

0  [2^  ^_(iy  "  CkXk)Hk(y)rkk{T,y)e"lyTdydT^ 


—  00  —  00 


=  xk  0  ^  1°  I  Hk(y)rkk(T'Y)e  ly^dydT;i 


(5.11) 


Taking  the  inverse  Fourier  transform  of  (5.9),  we  see 


that  K  (t)  involves  the  two  integrals 
KK 


1 

2tt 


rkk(T,x)  ioit  ^  ^  i 

-  e  dr  and 


2d  (ai) 


2tt 


r»rkk(T,U))  i*t 

-  e  doj 


Using  (3.36)  and  (3.40)  and  noting  assumptions  c.  and  d. 
we  have 


n  ,1  r°°  ikk(',Uj)  lout  , 

0  \  —  - - -  e  dx  ] 

2r  1™  2D.  (tu) 


r  n  S1  r°°  rkk(',uj)  iuL't  , 

Ck  0  ijF  J  -  e  dx] 


(5.12) 
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Hence,  combining  Eqs „  (5.11)  and  (5.12),  we  conclude 

that  the  first  term  in  (5.9)  gives  a  time  function  that  is  of 

the  order  of  C  /&  multiplied  by  the  order  of  the  time  function 
k  k 

resulting  from  the  second  term. 

With  this  result,  from  (5.9),  we  have 


kk 


1  00  0  °°  _  -j  ^  ^ ) 

(t)  =  —  J  j  J  (y-iCv^v)H  (y)e  Y  — - -  dy  dT  ] 


_  co  2irm  _oo  —  c 

k 


k  k  k 


D  (id  ) 
2 


iuit 

e  djj  (1+T^) 


2tt,ti  •)  J 


Tmk  -c 


(y-iCkui  (y) 


-iyT  ri 


rkk(T'”»  ia)t  ,  ,  , 

-  e  dju  j  dy  dT 

d2(.) 


(1  +  t3) 

ck 

where  0  (T_]  =  0  ( —  ) 

3  \ 

_  ......  .  1  j-"  rkk(T,U))  ittjt. 

Consider  the  integral  —  J  _  e  dju. 

2rr  _ qp 


be  written 


D2(iu) 


(5.13) 


This  can 


I  = 


-4  (t  ,  x) e^xtduu 


2 n  -cc  F  (UJ-2C ,  x,  i)  (x2-4iC  ID,  UU-4UJ2) 

K  K  K  K  k  k 


(5.14) 


from  assumption  d.  and  (5.7). 


71 


The  positions  of  the  poles 
of  the  integrand  are  shown 
in  Fig.  9.  Since,  by  assump¬ 
tion  a.,  the  poles  of  the  term 

— _ - i -  are  far 

(x  -4iC  X  X-4x  ) 
k  k  k 

distant  from  all  the  other 


Fig .  9 


poles,  the  treatment  given  in 
section  3.6  is  certainly 

applicable,.  Referring  to  section  3.6#  we  can  set 

2  2 
(x  -4iC  x  x  -  4x  )  =  G(x) 
k  k  k 

Fk(",x)  (x  -  2Ckxki)  =  F(x) 
n  =  1  +  p 
m  =  2 


6.  =  4C  x 

1  k  k 

e  =  4iu 

2  kD 


1 

a  =  0 


2x  «  y  ^ 
kD  2 


b  =  2Q  a),  i  +  2m  ,  b  =  2C  x  _ 

1  k  k  kD  2  k  k  kD 


i  -  2x, 


From  (3.69)  to  (3.71),  it  can  be  seen  that 
li”  rkk(T,x)elAtdlu 

i  =  r-|  — - 1  (1+t  )(i+t  ) 

\  2 "  -»  (x-2C  0)  i)  1  2 

k  k  k 


(5.15) 


where 
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o  [T  ) 


=  O  [ 


2(4CkUJk)P(2UJkD)2MAX  [|jk(bj)|] 

(1+p)  (4uJkD)  (2JJkD)P+1MIN[  |  Jk  (a  )  |  ] 


2P  MAX  [ | J  (b  )  |  ] 

O  {  - k  - } 

(1+p)  MIN  [|j (a.)|] 

*  D 


(5.16) 


O  {T  )  <  O  {  |  (- 


G(a-6  /2) 


G(a) 


")  G  (a )  |  } 


o  £C2} 

k 


(5.17) 


Finally  a  substitution  of  (5.15)  into  (5.13)  results  in 
the  approximate  Fourier  integral 


^  (i) 

kk 


»  (y-iC  \)H  (y)  A 

-  j  - -k--k--k -  5,.,.  (y.y-ui)dy 


2rrm  m  -a  ou  -  2C  U)  i 
K  k  k  k 


kk 


(5.18) 


which,  upon  inversion  gives  K^(t)  correct  to  within 
O  +  T^}  (neglecting  higher  order  terms  in  the  expansion 

of  (1+T^) (l+T^) (1+T^)),  where 


2PyP  MAX  [ | J  (b  . )  | ] 

0  { T  )  =  o  {  - k.  - 5-J - 

1  1+P  min  [ I J  (a  ) I ] 

*  j 


(5.19) 


o  (t  ]  <  o  {c  1 

2  k 


o  (t  }  <  o  {—3 

3  U) 

J  k 


(5.20) 

(5.21) 


It  should  be  noted  that  for  shot  noise  excitations,  the 
first  term  of  (5.7)  is 
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1  GO  I-.  («)) 

J  Hv(y)  - dy  = 


ikk(“,) 


2*rm  J  „  k 
k 


2D^ (m) 


4nm  D  (x) 
k  1 


J  Hk(y)  dy 


and  vanishes  identically.  Hence  the  error  to  which  this  term 
gives  rise,  is  zero. 

When  the  singularities  of  T  _  (t  ,  x)  lie  within  a  distance 

kk 

€  from  the  origin  such  that  e  >  2Q  x  ,  it  is  obvious  from 
k  k  x  k 

the  preceding  arguments  that  if  €  «  2 jlj  then  formulae  (5.18) 

k  k 

Gk 

to  (5.21)  are  still  valid  provided  Q  is  replaced  by  -  . 

k  2  m 

Thus,  it  is  now  possible  to  clarify  the  sense  in  which 

the  excitations  considered  here  are  “slowly  varying".  The 

excitations  are  slowly  varying  when  the  function  7  (t,x)  (of 

xx 

Eq.  5.5)  takes  significant  values  only  for  |  uu  |  <  €  where 

e  «  2uu  . 

kD 

To  illustrate  the  application  of  the  above  ideas  we  shall 
consider  the  following  example. 


Example  5.1 

Let  us  consider  once  more  the  situation  of  Example  4.1. 

For  this  case,  the  spectral  density  is 

My#y-a)  =  i(x)  •  i  (5.22) 

A 

where  l(x)  is  the  Fourier  transform  of  the  unit  step  function, 

U  ( t )  .  Since  U(t)  (x)  +  — —  ,  (5.22)  becomes 

ix 

«(y.y-x)  =  I(-6(x)  +  -X)  (5.23) 

ix 
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Substituting  the  first  term  of  (5.23)  into  (5.18)  gives 


„®  (y-iCkxk)Hk (y) 
-*  “-^kv 


I  n6(x) 


dy  e 


ixt  , 
dx 


=  -  1  1 

2  2-mQ2  (2QQ) 


j  (iy  +  C^)H(y)dy 

*00 


1 

2  2 

8rn  f<  £ 


(5.24) 


by  (3.11)  and  (4.47).  It  can  be  seen  that  the  exact  expression, 
Eq.  (4.48)  yields  the  same  result.  Thus,  for  the  first  term 
of  (5.23),  (5.24)  involves  no  approximation. 

Similarly,  in  general  the  expression  (5.18)  becomes 


Ml)  =  \  -*M-  (5.25) 

2m  Q  i  x-2C^i 

Substituting  the  second  term  of  (5.23)  for  I(x)  above  and 
taking  the  inverse  Fourier  transform,  we  have 


I  >“  i  e1X  dx  _  ii  r  1  _  e  2^nt- 

2- (2m2:.2)  -®  ii(x-2Cni)  2m2f.2  4CCU  2;ni 

=  -  [1  -  2e'2Cnt]  (5 

8m2cn 

/^(t)  is  just  the  sum  of  expressions  (5.24)  and  (5.26)  with 
suitable  correction  terms: 


<(t)  = 


4m2C03 


e  2"nt)  (1+T  +T  )) 
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and  from  (5.19)  and  (5.20): 

O  (T  )  =  O  [Cl 

0(t2)  =  o  £c2} 

since  J(x)  =  1,  and  P=1 . 

Thus  : 

K(t)  =  1 -  [1  +  1  -  e"2COt(l  +  T)]  (5.27) 

4m2C03  2 

where  T  is  of  order  Q .  This  result  can  be  directly  verified 

by  comparing  it  with  the  exact  solution,  (4.51).  In  Fig.  10, 

the  two  expressions  are  shown  graphically.  Here,  the  quantity 
2  3 

—  ~  is  plotted  as  a  function  of  ftt .  Now,  when  the  unit 

1  -iuut  ^ 

0  A 

step  commences  at  t=t^,  by  Eq.  (3.5),  I  (t-t^)  «-♦  e  I('j)). 

And  from  (5.25)  we  see  that  K (t)  is  merely  shifted  in  time  by 
tg.  Hence,  for  a  white  noise  excitation  modulated  by  a  step 
function  beginning  at  t^: 


*(t) 


r 


i 


(1  -  e-2^Q(t-to)  )  (l+T)  ,  t  >  t 


4m  QQ 
0 


(5.28) 


t  <  t 


where  O  [t]  =  O  [Q  +  Q  }  . 


Example  5.1  shows  the  validity  of  Eq.  (5.18),  when 

C  (T,t*)  has  points  of  discontinuity  but  is  otherwise  slowly 
Qij 

varying  in  t*.  For  suppose  that  C  (T,t*)  has  a  finite  jump 

Qij 

of  magnitude  C  at  t*  =  t*.  'I'hen  we  may  write 
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ZKlt)n?Cl 


O 

* 

O 

LL 


t= 

© 


in 


£ 


o 


in 


O 
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ill  3TT 


(5.29) 


(c) 

cn  (T,t*)  =  C  (T,t*)  +  CU(t*-t  ) 

Qij  o 

<  =  > 

where  (r,t*)  is  continuous  in  t*.  The  generalized 

ij 

spectral  density  of  Cq  is 

ij 

*  A (c)  1  — imt 

(Y/Y-x)  =  (y ,  y-x)  +  C(r6  (jj)+  0  (5.30) 

Qij  Qij  ^ 

(c)  (c) 

where  $  corresponds  to  C  in  (5.29). 

°i3  ^ 

But  the  second  term  of  (5.30)  is  the  spectral  density 
of  Example  5.1  with  I  =  C.  Then  from  (5.18)  and  (5.28): 


"kk't) 


1/2 T  (y-iCk'J;k>Hk(''>  1(C),  _  i«tj 

=  — L -  J  -  I  (y,  y-o))dy  e  dtu 

kk 


2~W  ""  ,J)  “2CkU,ki 


+  ~  (1  ~  exp(-2Ck«M  (t-tQ)))  (1+T) 


.  2,  3 

4m  ck\ 


(5.31) 


5 . 4  Approximate  Variance  Formulation  in  the  Time  Domain 

In  some  cases,  it  may  be  advantageous  to  work  with  the 
slowly  varying  case  in  the  time  domain;  that  is  to  express  the 
inverse  transform  of  Eq.  (5.18)  in  a  general  form.  We  can 
write  (5.18)  in  the  form 


\k(i)  = 


A  (id) 


Vk(*_2CklV)i 


^  CD  A 

A(x)  =  —  J  (iy+"kXk)Hk(y)  $kk(y'y_x)dy 


(5.32) 
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From  the  frequency  convolution  theorem  and  the  definition  of 
the  generalized  power  spectral  density: 


j  [C  OJ  h  (t-T)  +  f~  h  (t-T)]c_  (t.T)dT  A  (uu)  (5.33) 

-kkk  ^4-  V  o_. 


0 


dt  k 


'Q. 


kk 


Substituting  (5.33)  into  (5.32)  and  employing  the  time  con¬ 
volution  theorem,  we  have 

.  t  t  -Cx  (t-t  ) 

Kkk(t)  =  - 2  I J  e 


m  uu 

k  k 


k  k  k  '  1  2 


+  k1  h<Vt2n  CQkk<tl' t2)dtldt2 
"W1  _  1 


(5.34) 


where  the  fact  that  U(t)  e 
From  Eq.  (2.22)  : 

[C.  X  h  (t  -t  )  +-rf— h(t  -t  )] 

kkk  1  2  dt^  1  2 


was  used. 


l‘-2!iVu 


-C  x  (t  -t  ) 

=  - e  k  2  Sin  J'vn(ti  _t2^ 

®kD”k  k  k  kD  1  2 

•ww, .  „  „  . 

-  -k\  e  sin  WW 

-C  x  (t-t) 

+  x  e  cos  x  (t  -t  ) 3 

kD  kD  1  2 

i  -C.x  (t  -t  ) 

=  —  e  k  cos  uu  (t  -t  ) 

m^  kD  l  2 


Substitution  of  this  result  into  Eq.  (5.34)  gives 

3  ti  -2C  x  (t-t  )  -r  x  (t  -t  ) 

1  npi  k  k  i  kkl2 


‘kk(t)  2*2  J  J  6 

"HA  °  0 


COS  X  (t  -t  ) 
kD  1  2 


'  C°kk<tl't2)dtl'dt: 


(5.35) 
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The  region  of  integration  is 


Fig.  11 

(5.35)  is  symmetric  with 
integral  of  (5.35)  is 
entire  shaded  region  of  Fig. 


shown  as  the  cross-hatched 

area  in  Fig.  11.  But  from 

the  definition  of  CA  (t  ,  t  ) 

Qkk  i  2 

it  follows  that  it  is  sym¬ 
metric  in  t  and  t  .  Hence# 

1  2 

the  entire  integrand  in  Eq. 


respect  ^ ^  auu  «^ . 
one  half  the  integral  taken  over  the 
11.  'That  is: 


K 

kk 


(t) 


e~2^kxkt  t  t  C  \(t  +t  ) 

— 2~ 2 -  j  J  e  1  cos  (t  .t_J 

2mkUJk  °  ° 


kD  1  2  Qkk  1  2 


•  dt 

1  2 


(5.36) 


where  from  (5.19)  to  (5.21),  this  is  correct  to  within  terms 

of  order  Q  .  This  is  the  general  approximation  to  the  variances 
k 

of  the  response,  and  proves  to  be  more  convenient  than  (5.18  ) 
when  the  generalized  spectral  densities  cannot  be  easily  cal¬ 
culated.  The  following  example  illustrates  this  advantage. 


Example  5.2 

Suppose  that  the  simple  oscillator  of  Example  4.1  is 
excited  by  a  random  force,  Q(t),  for  which  the  correlation 
function  is  an  exponential  cosine  in  T  enveloped  by  a  slowly 
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varying  deterministic  function,  g(t).  That  is: 

-a  |  t  | 

^(t^t  )  =  CQ  e  cos  aQT  g(ti)g(t2)  (5.37) 

where  a  is  small  and  t  =  t^-t^  .  And,  the  equation  of  motion 


is 

5  +  2C  +  n2§  =  2-^-)  (3.4i) 

m 

It  can  be  seen  that  the  generalized  spectral  density  of 
Q  corresponding  to  (5.37)  is  quite  complicated.  Therefore, 
we  use  (5.36)  in  calculating  the  variance  of  the  response.  We 


shall  consider  only  the  case  for  which  ft  is  in  the  vicinity 

0 

of  0. 


Substitution  of  (5.37)  into  (5.36)  yields 

-2Cnt  t  t  cft(t,+t„) 


\k(t)  =  co 


J  I 


2m  n  0  0 


e  12  n  (t  -t  ) 

D  1  2 


-a  |  t  | 

e  cos  0  t  g  (t  )g(t  )  dt  dt 

0  ^1^2  ^  2 


-2£0t 


t  t  Cn(tx+t2)  -a | t2-t1 | 


O  2o2 

2m  n 


J  J 


0  0 


g(tl)9(t2) 


-  [cos  (nD+nQ)  (t2-t1)  +  cos  (nD-£V  ■'dtldt2 


But,  by  Reimann '  s  lemma,  when  I  I  <<;  (Gp+fl-)  ~ 


then  to  order 


iya0i 

2n 
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C„  e  *  t  t  £fl(t  +t  )  -a|t  -t 


Kkk<fc>  " 


-rr-iS  e 

4m  Q  0  0 


1  2  e  2  1  cos(0  -0  )  (t  -t  ) 

DO  2  i 


9  (t  )g  (t  )dt  dt 
1  2  12 


(5.38) 


Or,  if  we  let 

u  ■  1ST  •  €  '  c  tf  '  Ti  *  eVi  ■  t2 
0  0 

nQ 

r\  s  efljjt  and  X  =  —  =  1  +  eg 

“D 


efl  t 
D  2 


then  Eq.  (3.42)  becomes 

-2cnt 


h  n  (t.+t  )  -dx|t2-t  | 

^TTT  J  J  «  2  «  9(T.)g(T.) 


Kkk(t)  '  °202n2  2  J  J  e 

4m  QflG  0  0 

D 


.  cos  p  (t2”Ti )dT1dT2  (5.39) 


to  order  - .  This  is  in  complete  agreement  with  the 

2% 

result  obtained  in  reference  (7). 

This  concludes  our  special  treatment  of  the  variances.  In 
the  remaining  sections  of  this  chapter,  we  show  how,  for  two 
important  cases,  methods  analogous  to  those  in  preceeding 
pages  may  be  applied  to  the  approximations  of  the  entire  matrix 
of  second  moments. 


5 . 5  Applications  for  Shot  Noise  Excitations 

Suppose  that  a  structural  system  governed  by  Eqs.  (2.17) 
is  driven  by  a  vector  shot  noise  excitation.  Then  according 
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to  (4.43)  of  section  (4.6)  the  matrix  of  generalized  spectral 


densities  of  the  excitation  is  given  by 

U  (y,y-u j)]  =  [i  (uu)]  (4.43) 

Q-m  3  k 


where  I  (x)  is  the  Fourier  transform  of  the  intensity  function 

I  (t) .  Hence  the  components  of  the  matrix  [K  .  (x)3  are  given 
jk  kj 

by  (4.12)  in  the  form 

A 

/S  •  (U')  CD 

K  (a)  =  — 1 -  J  h  (y)  H  (uc-y)dy  (5.40) 

k3  2tt  k  ] 


By  Eq.  (3.20)  we  can  write  (5.40)  in  the  form 


I.  .  (uu) 

K  (a)  =  j 

K;|  2tti 


dy 


nyn^  (y-a^) (y-a2) (x-y-a^) (x-y-a4) 


where 


ai  =  i£k%  + 


a^  =  iC.  x,  -x 
2  k  k  kD 


a^  =  iQ.u.  +  x  *  a  =  i£.x.-x, 

3  3  3  3D  4  3  3  jD 


J 


r  (5-41) 


1  IV 

»  “HC  < 

X 

2 

Fig.  12 


When  y  is  replaced  by  z  the 
integrand  has  simple  poles 
at  the  points  shown  in  Fig.  12 
Integrating  over  the  upper 
and  lower  half  planes  along 
the  contours  of  Figs.  3(a) 
and  3(b)  we  have 
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(K,_J  = 


iTki(ui) 


[ 


u  2m,m.'X,  _  (a-i6  . -u  . )  (a;— i  6  . -v  .) 


'WkD 


kj  kj'  ~~kj  'k  j ‘ 
1 


(uu-i 6  +v  .  )  (uu-i 6  +u  . ) 
kj  kj  kj  k] 


3 


(5.42) 


(K  )  « 


il  .  (’Ju) 

Kj 


[ 


1  2m  m.uu.  (uu— i 6,  u,  )(uu-i6.  .+  v,  .) 

k  ]  ]D  kj  kj  kj  kj 


( jj— i 5  v  . )  (uj— i 6  .+  u  .) 

k3  k}  kj  k] 


1 


(5.43) 


where 


6n  .  =  C  u)  +  C  .uu . 

*3  k  k  3  j 


u  a  (uu  +  uu  )  y 

kj  kD  jD 


v,  .  a  (u>,  —  id  .  )  , 

kj  kD  jD  / 


(5.44) 


and  from  (4 . 34) ,  K  =  ±( (K  )  +(*).). 

k3  2  kj  u  k]  t 

Now,  we  assume  that  the  intensity  functions,  I  (t),  are 

3k 

slowly  varying  in  the  sense  that  the  singularities  of  I^(x) 

lie  within  a  distance,  e,  from  the  origin  where  c  ~  6^<<  ^kD 

or  x  .  We  also  assume  that  I  (x)  can  be  expressed  in  the 
jD  3* 


form 


v*>  - 


J.  .  (uj) 


(5.44) 


F  .  (uu) 
k3 


where  J  (z)  is  analytic  in  the  finite  z  plane  and  F,  .  (z)  has 
k  j  k  J 

p  poles  at  a  ,  ...,a  in  the  upper  half  z  plane. 

1  P 
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We  shall  first  treat  of  the  case  for  which  v  = 

kj 

(x  -  x  )  is  of  order  e.  Consider  the  first  term  within  the 
kD  jD 

brackets  of  Eq.  (5.42);  its  inverse  transform  is 


1  r 

J 


~  ,  .  ixt 

I,  .  (x)  e  d-ju 

kD 


(x-i6,  u  . )  ( uu-i 6  v  ) 
*3  kD  kj  kj 


(5.45) 


Referring  to  section  3.6,  we  can  set 

(x-i6  a  .)  =  G(uo) 
kj  kD 

a  =  i6  +  v 

kj  kj 


so  that  G(a) 
Also : 


m  = 


n  = 


6  = 


€  ^  = 


-2x  . 
JD 

1 

P+1 

€ 

1 


u. 


kj 


2  J 


(5.46) 


(5.47) 


And  substitution  of  (5.46)  and  (5.47)  into  Eqs .  (3.69)  to 

(3.70)  gives 


I(t)  =  -  T- 

2"  2jj 


I  .  (x)  e"*"1  fcdx 
kD 


-,n  x-i&  . — v 

kj  kD 


(1+T  +T^ ) 


(5.48) 
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1 


3 


1  e  p  lJ<i6k-i+  vki)  I 
o  It]  =o(4-  (-^-)  - ^ - — — 

1  P  1  Ukj  MIN  ( |j(a  ) |] 

(5.49) 

0  (T  )  <  O  (—) 

2 

We  can  get  similar  results  for  the  other  terms  of  the  inverse 

A  A 

transform  of  (K  (x))  and  (K.  (x))  .  The  results  are 

ij  u  ID  l 


-i/2rr 


(K  •  (t)) 


kj  u  4m  m  x  x 

k  j  kD  jD 


“  ixt  ~ 

J  e  I .  .  (uu ) 

_»  ID 


x— i6  v.  .  x-i6,  .+  v,  . 

kD  kl)  kD  kD 


-]  dx  (1+T  +T  ) 
1  2 


(K  .  (t)) 

kj  X 


-  J"  <><=  i  (uj) 


4m  m  x  x  -00 
k  j  kD  jD 


ID 


x-i6,  .  +  v  .  x-i6,  v.  . 

kD  kD  kD  kD 


]  dx  (1+T  +T  ) 
12 


Hence,  the  approximate  Fourier  transform  of  K^(t)  is 


\j(X)  = 


-il  (x) 
Dk 


[ 


3  (5.50) 


4m  m  ll  u:  x-i6  . +  v  .  'X-i6_  .  -  v, 

k  j  kD  jD  kD  kD  k3 


fhis  gives,  upon  inversion,  the  second  moments  correct  to  with- 

i  p  lJ<i6k-i+  I 

in  terms  of  order  ( -  (— €  )  - ^ ^ c  }  . 

P+1  Pkj  MIN  (|j(ak)|3  Pkj 

We  shall  now  express  (5.50)  as  a  convolution  integral. 


Since 
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2  e 


-6,  .t 

kD 


cos  v  .t  ♦— *  -  i  [ -  +  - 

uj-i6  .+  v  .  u)-i6  -  v 


kj  kD 


kj  kj 


by  the  time  convolution  theorem,  we  have 

i  -t 


*kj  2m  m  x,  uj  .  J  J'jk 

J  k  j  kD  jD  0 


-6..t 

J  I  e  -*  cos  <3t} 


kj 


where 


(1+T  +T  ) 
1  2 


p  |  J(i6,  >  v,  J  | 


JSJ _ *2_  } 


O  {T  ]  =  0  [  (~— ) 

1  P+1  pkj  MIN  UjC^)  I) 


O  {tJ  =  O  {-£-} 


> 


Hkj 


(5.51) 


(5.52) 


and  for  small  v  .  • 
kD 

Now  suppose  that  0(v,  )  ~  0(x,  )  or  0(x  )  and  hence  that 

KJ  K  D  jD 

0(v  )  »  e.  Consider  the  first  term  of  the  inverse  transform 

kj 

of  Eq •  (5,42);  given  by  Eq.  (5.45) 


Kt)  =  -±  J 


I  .  (x)e1Uitdx 
kD 


2"  _cd  (x-i6  uj  (x-i6, 
k  j  k_ 


(5.45) 


kj  kj 


In  this  case,  we  can  set 

°m  - ukj,(*-itkj-  v 

a  =  0 

so  that 

G  (a )  =  ( i  6  +  u_  .  )  ( i  6  .+  v  .)  (5.53) 

kj  kD  kD  kD 
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Also  as  in  the  previous  case: 


m  =  2 


n  = 


6  = 


€  = 


P 

e 


b  ,  b 
1  2 


la  v  |  =  2x  .  >2  MIN  ( x  ,  uu  ) 

k}  kj  1  3D  -  kD  jD 


kj 

Ukj 

i6,  .+  u  ,  b  ,b  =  i6  +  v  . 
k}  k}  3  4  kj  kj 


>  (5.54) 


And  substitution  of  (5.53)  and  (5.54)  into  Eqs.  (3.69)  through 
(3.70)  gives 

1  ~°°  A  i  fllf- 

(5.55) 


 1 1 
2n  G(a) 


I(t)  =  -i  — ±-  J  I  .  (x )  e  dx  (1+T  +T  ) 

—  cr>  ^3  1  ^ 


o  {t  ]  <  o  {. 


eP  1  MAX  [  |  J0\)  I  3 


p  MIN(xkD'ju;jD)  vP.  MIN  [|j(bk)|] 


O  {T  }  <  0  {^} 

VkD 


(5.56) 


Performing  the  same  analysis  for  the  other  terms,  we  have: 

21  (t)  6  . 

K,  ,  (t)  =  - ^-—2 - J1  (1+T_  +T_+T.,  )  (5.57) 


kj 


m  m  u  v 
k  j  kj  kj 


1  2  3 


where  the  orders  of  magnitude  of  T  and  are  given  by  (5.56) 


and 


O  (T  }  =  O  {(-^)  ) 

3  ^kj 


(5.58) 
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1 


for  the  terms 


T  arises  from  the  approximations  made  of 
3  G(a) 

in  (5.42)  and  (5.43). 

Equations  (5051)  and  (5.57)  form  the  complete  approximate 
response  formulation  to  vector  shot  noise  excitations.  We 
consider  one  application  in  the  following  example. 


Example  5.3 

Suppose  that  we  have  obtained  the  second  moments  of  the 

response  of  a  system  by  experimental  measurements,  and  that 

we  know  the  excitation  is  a  vector  shot  noise  but  the  specific 

intensity  functions  are  unknown.  Eqs.  (5.51)  and  (5.57)  make 

it  possible  in  certain  cases  to  estimate  quite  easily  these 

intensity  functions  from  a  knowledge  of  K  .  (t). 

*3 

Let  us  first  consider  the  case  v  o{e),  so  that  we 

kj 

employ  Eq.  (5.51): 

1  t  ~6k’T 

K  (t)  ~  - -  J  I  k(t-T)  e  3  cos  v  T  dT 

^Vji  o  3k  k;l 


VjWiD  'o 


t  ,  ,  '6kj(t_T) 

I  (t)  e  J  dT 


(5.59) 


.o  order 


11 

^3 


Assuming  that  I  (t=0  )  =  0  and  by  repeatedly  integrating 
3k 

(5.59)  by  parts,  we  obtain  the  series 

V  (t)  ~  r~~ — r —  '_I-ik{t)  ~  r  [tk(t)  '  e~  tbk(0):i+***:i 

kD  2r\mj'xkDXDD5Dk  3  6  3  3k 
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Now,  that  I,  ■  ( t )  is  slowly  varying  implies  that  its 

V 

derivative  is  small.  With  a  =  O  {— — }  where  a  is  small,  we 

Ij^ 

have 

I  (t)  ~  (2m.  m  m  uu.  6  K  (t))  (5.60) 

lk  k  ]  kD  jD  I)k  jk 

correct  to  order  a. 

A  consideration  of  (5.57)  shows  that  for  the  case  v  »  e, 

kj 

we  have  2  2 

ro.  m  #  u,  .  v,  . 

I  (t )  ~  - j — ^  K  (t)  (5.61) 

3k  2  6 

kj 

Thus,  with  K  (t)  =  0  at  t  =  0,  the  intensities  may  be 

kj 

estimated  from  a  direct  examination  of  the  second  moments. 

5 . 6  Applications  for  Filtered  Shot  Noise  Excitations 

Here  we  treat  a  case  of  the  greatest  importance  to  the 
estimation  of  Titan  III-C  flight  loads. 

Suppose  that  a  shot  noise  with  intensity  function  I(t) 
is  the  excitation,  Q (t ) ,  to  the  system  of  (3.41): 

5  +  2C  r4  +  n2§  =  (3.4i) 

m 

The  displacement,  is  called  a  filtered  shot  noise.  Now 

let  each  normal  coordinate  of  an  N  degree  of  freedom  system 
be  driven  by  the  displacement  of  (3.41).  That  is: 

1  +  2C  x.  L  +  X2;.  =  ^1  (k=l,  N)  (5.62) 

k  k  k  k  k  k  rt^ 
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We  shall  apply  Eq.  (4.12)  and  the  methods  of  the  previous 

sections  to  the  approximate  calculation  of  the  second  moments 

of  the  §  . 

k 

From  (3.29),  the  generalized  spectral  density  of  §  is 

A  A 

(w.  ,u>-)  =  H  (x  )  $  ('JO,  ,uu  )H*(uj  )  (5.63) 

1  2  ls±2  2 

and  by  (4.43)  the  spectral  density  of  the  shot  noise, 

A 

i  ( il  ,  x  )  is  given  by 
s  1  2 

VW  =  ^xi-UJ2)  (5.64) 

Thus : 

i_(x.,ii^)  =  H(X-  )I  (U)  -U)  )H*(tD  ) 

1  2  1  1  2  2 

or 


where 


My»y-n  =  H (y)  I  (x)H*  (y-x) 


1  2 

H  (x)  =  -  (f/ 
m 


2  -1  ^ 

-x  +  2iCQx) 


I(x) 


I(t) 


-ixt 
e  dt 


► 

J 


Thus,  for  this  case  (4.12)  becomes 


(5.65) 


(5.66) 


<k  •  (■L)  =  I(a)u  H(y)H(y)H(x-y)H  (x-y)dy 

—GO  3 


since  H  (y-x)  =  H(x-y). 

From  Eq.  (3.20)  this  can  be  written  as 


K 

r 


*3 

(x)  = 


=  r (x)i (x) 

1  >“ 

—  H  (y)H(y)H(x-y)H. (x-y)dy 
2"  k  ] 


(5.67) 


(5.68) 

(5.69) 
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In  what  follows,  we  shall  attempt  to  evaluate  approxi¬ 
mately  the  inverse  transform  of  T (x)  (denoted  by  Y(t))  so 
that 

t 

Kk  (t)  =  j*  I(t-T)y(T)dT  (5.70) 


As  in  the  last  chapter, 
are  small  (of  the  order  of  10 
varying  function. 

From  (3.20),  F(x)  may  be 

r  (j )  =  L-  _ • 

2~  1®  (y_a  )  (y-a2) 


we  assume  that  £ ,  £  ,  and  C 


-2 


)  and  that  I(t)  is  a  slowly 


expressed  as 

1/mj 

(x~y-a3) (x-y-a4) 


l/m 

(y“a5 ) (y-a  ) 


dy/m 


(x-y-a  )  (x-y-a  ) 
5  o 


(5.71) 


where 


i£,  x 

+  X 

,  a  =  i  £  x 

-  X 

k  k 

kD 

2  ~k  k 

kD 

iC  .x  . 

+  X 

,  a  =  i  £  .  x  . 

-  x  , 

3  3 

ID 

4  3  3 

3  D 

i  £  n  + 

n 

,  a^  =  i£Q  - 

n 

D 

6 

D 

> 


y 


(5.72) 


We  first  assume  that  lx,  -x  I,  lx,  -F  |  ,  | x .  -H I  «  x,  , 

1  kD  JD1  1  kD  3  D  1  kD 

x  ,  or  0  o  Upon  replacing  y  by  z  and  employing  the  calculus 
jD  D 

of  residues  the  integral  of  (5.71)  can  be  evaluated.  Subsequently, 
the  inverse  transform  of  F(x)  can  be  calculated  approximately, 

* 

according  to  the  methods  of  preceding  sections.  The  result  is 

*  See  Appendix  III. 
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V  (t) 


[Re  (- 


H  (a  )  e 


i  (a  -fa  )  t 
1  4 


m  m  m  8  jd,  uj  .  0 
k  J  kD  ]D  D 


a4  -  a6 


i (a  +a  )  t 
H  (a  )  e  H  (a  )  e 

+  Re  ( - — )  -  Re  ( - - — 


i  (a1+a6)t 


a  -  a  a  -  a 

2  6  .  ,  .  4  6 
i (a  +a  )t 

H  (a  )  e 

-  Re  ( - )  ] 

a^  -  a 
2  6 


1  Hk<a5>e 

+  — - — 7  [  Re  (---5- 
8x  Ci 


i (a  +a  )t 
4  5 


-)  -  Re  ( 


Va5> 


i  (a  +a  )  t 
.  5  6 


JD 


a  .  -  a 
4  , 


a  —  a 
4  6 


1  Hi(a5)e 

+  — - - -  [  Re  (—3 - - - 


i (a  +a  )t 
2  5 


8"kDn 


H.  (a  )e 

)  -  Re(J - 2_ 


i  (a  +a  )  t 
v  5  6' 


)3 


-)  3 


a  -  a 
2  6 


a_  -  a 
2  6 


(5.73) 


correct  to  terms  which  we  denote  by  T  +T  • 

12 


Eq.  (m  .15)  : 


According  to 


0  [1^3  < 

2[MAX  ( |  a}  )|3P  lMAX(|a(  +  )-a(’)  |  j  MAX[  |j(a(  +  ))  |] 

0  ( - £ - ^ - - - 

(p+2) [MIN (la ^-a  ^ ) ][MIN( | a^  +  ^-a  | ) JP+ZMIN[ |  J(a  | ] 

k  j  k  j  K 


0  [T  }  <  [ 
2 


2  MAX  ( | a^  } | ) 
MIN  ( |a^+) | ) 


(HI. 15) 


where  these  are  the  orders  of  magnitude  relative  to  y  (t)  .  This 
can  be  approximated  still  further  when  we  note  that 
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H(V  -  un2-  \  +  2Vk  -  2CnW  -  2inD(cn  ‘W3’1 


-  C(t!„  -  *2  +  cV  +  c2'.!2  -  Jcnc.n  )  -  2in  (Cti-C.  o),  )] 

DkDkk  kk  D  k  k 


-1 


[(nD  -  “kD1  -  2inD(C°  -  Vk)fl 


to  order  C  -  Now  we  expand  this  in  a  Taylor  series  in  terms 


of  the  quantity  (QCl  -  C,UU,  ): 


H(a )  _  ~z  *2  +  ~  ? 1  2  £ 

-L  n  n  . 


2inD<r-n-Ck\» 

-  J  •  •  • 


o  -a).  o~-x,  (n^  -  uu?  ) 

D  kD  D  kD  D  KD 


Since  the  £’s  are  small,  this  series  converges  rapidly  and  we 
can  say  that 


H<V  *  r^-T- 

Q  -  uj, 

D  kD 


(5.75) 


to  order  less  than  MAX(C«  v)  since  the  magnitude  of  the  term 

C 

2iVn-  r  V 

- - -  is  of  order  one. 

q2  ~ 

D 


Also  if  we  expand 


a  -a  k  (C  w  -  £fl)  +  (f2  “  0)  ) 

46  J  3  D  jD 


in  terms  of  the  quantity  Q  ju  -  :'Q  ,  we  obtain 

j  j 


L- 


i(ClXj’Ca) 


(a  -a  )  fi  -X  Q  -tu  Q  -a) 

4  b  D  JD  D  jD  D  3D 


-1  +  ... 
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i(C.u)  -CO) 

Since  the  magnitude  of  - 3 — 3 -  is  of  order  MAX(C,C.): 

fl  -uj  .  j 

D 


1 


a  -a 
4  6 


n  -uj  . _ 

D  JD 


(5.76) 


to  order  MAX  (C/C.)» 

3 

Combining  (5.75)  and  (5.76)  we  have 


»<*!>  i 

vae' 


(5.77) 


to  order  less  than  MAX(C,C  /C  )» 

3  k 

By  a  similar  procedure,  we  can  obtain 


H(a3) 

a  -a 
2  6 


VS5> 

a  -a 
4  6 


W 

a  -a 
2  6 


(n2-uj2  )  (nT-uT  ) 

D  jD  D  kD 


(uj2  -n2)  (n  -uu .  ) 

kD  D  D  3D 


U  -D-nD>  <  VW 


(5.78) 


(5.79) 


(5.80) 


all  to  order  less  than  MAX  (C/C./C  )• 

3  k 

Finally,  by  virtue  of  Eqs.  (5.72),  Eq .  (5.73)  and  Eqs. 

(5.77)- (5.80) ,  we  have: 


95 


-  (Ck^k+G -j,JJ  ) t 

Y.  .  (t)  ~  A[  B,  e  J  3  cos  (id  -uu  .  )t 

ki  1  kD  jD 


"\ 


-(■'■oj  +cn)t 

-B  e  cos(uj,  -n  )t 

2  kD  D 


B3  G 


-  (G  .uu  .  +Gft)t 
3  3 


cos (x  -Q  ) t  +  B  e 
jD  D'  4 


-2GClt 


A  = 


2  2  2  2 

8i  •  )  (A  -UJ,  ) 

kD  ] D  D  D  jD  D  kD 


(5.81) 


B  =  20  +  uu,  +  UJ 

1  D  kD  jD 

UJ  .  _ 


B2  -  n  (^n  +  uu,  )  +  (ft  +  U)  •  n) 

A  ftD  D  kD  D  jD 


UJ. 


kD 


B3  nn  'D  jD 


(n-+  UJ-)  +  (nD  +  \d) 


ui  .  U) 

B  —  — ( ’JD  +  Q  )  +  — — —  ( fJU  +  Q  ) 

4  Q  V  kD  DJ  Qn  jD  V 
D  D  J 


J 


which  yields  a  result  correct  to  (T  +T  +T  )  where  O  [T  }  and 

12  3  1 

O  {t^}  are  9^VGn  ky  (III. 15)  and  arises  from  the  approxi¬ 
mations  of  the  frequency  response  functions  so  that  O  {T^}  < 

O  {MAX(£, £,£.)}• 
k  D 

It  is  instructive  to  note  the  results  for  the  following 
special  cases. 

A.  j=k  : 

From  (5.81)  : 

\*(t)  "  8.2  ~  *2  2  2  t2'V  *W  2  ^  + 

8x  w_(U  -uu  ) 
kD  D  d  kD 


(equation  continued ) 
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(Q  +  u)  )  e 
D  kD 


-(C  u)  +CO)t 
k  k 


cos  (uu,  -  Q  )t 
kD  D 


*  ^  +  0  ,  e-2Cnt] 

n  kD  d 

D 


(5 


B.  id  .  =  n ,  Q  .  =  c 
-D -  — ] - 

Letting  x  =  Q  -A ,  where  A  is  small,  Eq.  (5.81) 
JD  D 


becomes 


8\dV2VA<V”W 


_ _ _ _  [  (3Q  +u)  )e 

1  D  kD 


-(CkVcn)t 


COS(\D"VA)t  "  PnD  +  \D}  S 


- (c  »  +cn)t 

k  k 


cos  ( x  -n  )t  -  (-^  20  +  G  +  ttt  )  e  2^ntcos  At 
kD  D  nD  D  D  kD 


+  2?i  +  n  +  U)  )  e  2^nt  ]  +  T(A) 

Q  D  D  kD 
D 


where  T(A)-*0asA-*0*  Hence: 

3Q  +  ua 

vv.  (fc)  ~  - 

A  -0 


LIM  v  (t)  ~  - 3-  -■  r~ 

K3  i6x  cr(n2-u;  ) 


,,n  ~(CV«J  +C0)t 
D  kD  e  k  k 


kD  D  D  kD 


•  LIM  {i  (cos (xkD~nD+A)t  -  cos (xkD~nD)t) } 


A— 0 


3x  +n 

kD  D 


-2cnt 


3  2  2 

16x.  w  (Ci-X  ) 
kD  D'  D  kD 


LIM  (— (1-cos  At)} 

A—0  A 


-(Vk+cn)t 


16®kDD(°D  -  XkD> 


tsin(GD  -  U^Jt  (f 


.82) 


.83) 


97 


by  L 'Hospital  1 s  rule. 

Hence,  when  the  central  frequency  of  the  filter  coincides 

th 

with  the  damped  natural  frequency  of  the  j  mode : 


Y.  (t) 


- (Cvw  +cn)t 

(3VW  «  * 

1S“*D°£  <nD  -  \D> 


t  sin(H  -x  )t 
D  kD 


(5.84) 


From  Eq .  (5.83);  with  yj  =  Q  -A 

kD  D 


V  .  (t)  =  — ^-r t  e  sin  At  +  T(A) 

8  a4  a 

OD 


where  LIM  T('-)  =  0  „  Hence, 

A— 0 

-2£nt 


LIM  Y,  .  (t)  = 


t  e 


Vn 


k j  '  *  „  ^4 


8  a 


OD 


t2  e-2Cnt 


LIM  (SinAt) 

A— 0  A 


8  a 


OD 


Thus,  for  the  case  au  .  =  0,  C.  =  id.  =  £.  =  £  : 

3  3  k  k 


Y,  .  (t) 


t2  e'2Cnt 


8  n 


OD 


(5.85) 
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CHAPTER  VI 


ESTIMATION  OF  PAYLOAD  STRUCTURE 
FLIGHT  LOADS 


6 .1  Introduction 

In  Chapter  V  the  basic  method  of  approximate  calculation 
of  the  second  moment  response  was  presented.  A  variety  of 
simple  problems  were  considered  in  order  to  illustrate  the 
basic  application  of  the  approximation  technique.  In  this 
chapter  the  more  realistic  problem  of  the  estimation  of  flight 
loads  in  a  payload  is  treated.  In  particular,  we  return  to 
a  consideration  of  the  problem  that  was  introduced  in  Chapter  I. 

We  take  note  of  the  following  basic  facts: 

1.  The  dynamic  inputs  to  the  payload  system  are  assumed 
to  be  the  booster/payload  interface  accelerations  of  the  six 
rigid  body  degrees  of  freedom. 

2.  These  interface  accelerations  were  obtained  by  a 
numerical  calculation  of  the  response  of  an  analytical  model 
(consisting  of  assemblage  of  discrete  elements)  of  the  booster 
structure.  These  booster  models  were  driven  by  loads  derived 
from  chamber  pressure  measurements  of  the  rocket  engines  during 
27  test  firings. 

Thus  the  subensemble  of  excitations  to  the  payload  system 
consists  of  the  27  sets  of  six  acceleration  time  histories. 
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We  assume  that  each  set  is  a  sample  of  a  random  vector  process 
whose  vector  components  are  the  six  rigid  acceleration  time 
histories  at  the  interface. 

In  Section  6.2  we  present  the  basic  deterministic  equations 
for  the  calculation  of  loads  in  the  payload  structure.  In  6.3 
and  6.4  we  discuss  the  assumptions  made  about  the  interface 
accelerations  and  express  these  assumptions  in  analytical  form. 

In  sections  6.5  and  6.6  we  show  how  the  statistical  parameters 
of  the  acceleration  process  can  be  estimated  under  the  assump¬ 
tions  made.  Finally  in  section  6.7  is  given  the  approximate 
formulation  of  the  second  moment  response  of  a  payload  to  the 
booster/payload  interface  accelerations.  This  chapter  concludes 
with  the  application  of  the  approximate  formulation  to  a  simple 
two  degree  of  freedom  payload  structure.  This  example  exhibits 
all  the  main  characteristics  of  the  calculation  of  response  of 
a  real  payload  structure. 

6.2  General  Formulation :  Equations  of  Motion 

We  assume  that  the  modal  frequencies  and  corresponding 
mode  shapes  of  the  payload  structure  are  known.  Then  the  problem 
may  be  formulated  entirely  in  terms  of  the  normal  coordinates 
of  the  payload.  According  to  assumption  B,  section  1.1,  the 
rigid  body  displacements  are  prescribed  at  the  booster/payload 
interface.  This  case  is  exactly  that  considered  in  Example  2.2, 
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section  2.3.  Thus,  from  Eq.  (2.19)  the  equations  of  motion 
of  the  payload  system  are 

[hi  +  2!M.^  Ch}  +  C *n)  =  [C]  [a]  (6.1) 

ii  i 


where 


[Til  ^  normal  coordinates  of  the  payload 

=  [cp]  [m3[u},  where  {u}  is  the  payload  relative 

displacement  vector  and  [m]  is  the  mass  matrix. 

tr  the  modal  frequencies 

i 

C  =  the  modal  damping  coefficients 
i 

{a}  =  the  booster/payload  interface  accelerations 

=  (x,  y,  z,  6  ,  9  ,  6  } 

x  y  z 

[C*  is  a  matrix  of  constants 

The  member  stresses,  { S} ,  are  related  to  the  normal  coordinates 
by 

[s)  =  [G][h}  (6.2) 


where  [g^  is  some  matrix  of  constants. 

In  what  follows,  we  shall  assume  that  the  modal  damping 

-2 

coefficients  are  of  order  10 

Now,  the  interface  accelerations  are  idealized  as  a 
random  vector  process.  We  wish  to  find  the  matrix  of  the 
second  order  centra]  moments  of  the  response,  {,n},  from  a 
knowledge  of  the  second-order  statistics  of  the  interface 
accelerations.  Once  this  is  accomplished,  the  second  order 
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(6.3) 


moments  of  the  member  stresses  are  given  by 


where 


[s.  .3 

ID 


[g3[k  •][g]T 

k  l 


[s.  .  ]  ^  the  second  order  central  moments  of  the  member 
ID 

stresses 


[K  j  the  second  order  central  moments  of  the  response 
kl 

£r»3 


Thus,  the  major  problem  is  the  determination  of  the 
statistics  of  the  normal  coordinate  response.  ihe  first  step 
in  this  task  is  the  description  of  the  interface  accelerations 
considered  as  a  random  process.  In  this  regard,  we  establish 
certain  basic  assumptions  in  the  following  section. 


6 . 3  Idealization  of  the  Booster/Payload  Interface  Accelerations 

Ihe  basic  assumptions  about  the  nature  of  the  booster/ 
payload  interface  accelerations  were  presented  in  section  1.3. 

We  restate  them  here  in  the  following  form. 

A.  (i)  Ihe  interface  accelerations  are  the  response  of 
a  lumped  parameter,  second  order,  linear  system 
to  random  excitation.  We  assume  thau  chis  linear 
system  is  the  structure  of  booster  itself. 

(ii)  Ihe  second  time  derivatives  of  the  random  excita¬ 
tion  to  the  booster  can  be  idealized  as  a  non¬ 
stationary  vector  shot  noise  process. 
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B.  The  intensity  functions  of  the  components  of  the 
vector  shot  noise  are  the  same  to  within  multi¬ 
plicative  constants . 

C.  The  excitations  to  different  normal  modes  of  the 
booster  structure  are  uncorrelated. 


These  assumptions  must  be  expressed  in  mathematical  form. 
Consider  assumption  A(i),  and  suppose  that  the  booster  structure 
has  N  degrees  of  freedom.  According  to  this  assumption,  the 
response  of  the  booster  to  the  initial  random  excitation  may 
be  expressed  by  an  equation  of  the  form  of  Eq.  (3.29).  That 
is, 

[i-  (x  ,x  )]  =  O-HL  (uj  J'vld  (x  ,x  (x  (6.4) 

'b  1  2  *>k  1  Q  1  2  bk  2 

(NXN)  (NXN)  (NXN)  (NXN) 


where 


[$_  (x  ,  x  )  J  =  The  matrix  of  generalized  spectral  densities 

;b  1  2 

of  the  normal  coordinate  response  of  the  booster  structure. 

th 

H^^(x)  The  frequency  response  function  of  the  k  booster 
normal  coordinate,  in  the  form  given  by  Eq .  (3.20). 


“^Q^X1'X2^  H  'Ihe  matr^-x  generalized  spectral  densities 
of  the  initial  random  excitation,  assumed  to  have  p  non-zero 
normal  coordinate  components. 

Employing  the  time  differentiation  theorem  for  Fourier 
integrals,  we  see  that  the  generalized  spectral  density  matrix 


103 


of  the  acceleration  response  of  the  booster,  [$  (x  ,  uj  )]  is 

ab  1  2 

given  by 

[$  (x  #x  )]  =  t^H.  (x  )0[x  x2  $  (x  ,x  )](Sh 

ab  1  2  \  1  1  2  Q  1  2  X  2 

(NX)  (NXN)  (NX)  (NX) 


(6.5) 


2  2 


where  xs  t*le  9enera^^-ZG<^  spectral  density 

matrix  of  the  second  time  derivatives  of  the  random  excitation 
components.  Now  the  interface  accelerations  are  related  to 
the  booster  normal  coordinate  accelerations  by  a  matrix  trans¬ 
formation.  Denoting  the  interface  accelerations  by  [a^] ,  and 


the  booster  accelerations  by  [a  },  we  have 

b 

Uj)  =  [C]  {a^} 

(6x1)  (6xN)  (Nxl ) 


(6.6) 


where  [c]  is  a  constant  matrix,  and  the  acceleration  components 
are  designated  as  follows : 


ui}  = 


/ •• 

X 

qi 

y 

q2 

z 

i 

q3 

’6x 

q4 

0 

y 

q5 

6 

zj 

Iq6i 

Then  from  Eq.  (6.6)  and  (6.5)  and  the  definition  of  the 
generalized  spectral  density,  we  have 
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/N  - 

whore  [$  (uu  ,  x  )  J  is  the  generalized  spectral  density  matrix 

a-r  1  2 


of  the  interface  acceler at ions <>  Eq .  (6.7)  now  embodies 

assumption  A(i). 

According  to  the  definition  of  a  general  vector  shot 

noise,  Eq.  (2.59);  assumption  A(ii)  states  that  the  covariance 

matrix  of  the  second  derivatives  of  the  random  excitation  has 

the  form  [i  (t)]6(t  -t  )  where  the  I  (t)  are  the  general 
jk  2  1  jk 

shot  noise  intensity  functions.  Hence,  according  to  Eq. 

(4.43)  of  section  4.6: 


(6.8) 


Eqso  (6.7)  and  (6.8)  now  yield  the  analytical  statement 
of  assumptions  A: 


(6XN)  (Nx  N)  (NXN)  (Nx  N)  (NX  6) 


(6.9) 


Now,  assumption  B  states  that 


I3k('Ll_J'2)  “  1  (xl_'x2)rKjk  ' 


(6.10) 


where  I (x  -x  )  is  the  transform  of  the  common  intensity 
1  2 
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function  and  [K  is  a  matrix  of  constants.  Thus,  if  assump¬ 
tions  A  and  B  are  valid,  we  have 

[3  (x  ,x  )}  =  [c]|>H,  (X  )<][K  11>H*  (x-  )<j[c]T-I(x  -x  ) 

aI  1  2  Jk  \  2  12 

(6.11) 

from  Eqs.  (6.9)  and  (6.10). 

Finally,  according  to  assumption  C.  the  cross-covariances 

between  different  normal  coordinate  components  of  the  initial 

random  excitation  are  all  zero.  Thus,  from  Eq.  (6.8)  and  the 

definition  of  the  cross-spectral  density, 

[I  (X  -x  )]  =  !>I  (X  -X  JO 
3k  1  2  kk  1  2 


or. 


(6.12) 


Substitution  of  this  result  into  Eq.  (6.11)  results  in 


[3  (x  ,x  )■)  =  [C][^H.  (xXlt-'K  .  (x  )'<l[c]T»i(x  -X  ) 

ai  1  2  \  1  31  X  2  1  2 

But  since  equals  -/k j  j  "'■s C''  v'K we  have 


[$a  (x1'x2)^  "  CcK^v1T‘;j'JC,‘Hb^(iu1)x3C,*H^(*^)^3l>^^ 


11 


[Cj  I(i1-x2) 


(6.13) 


Similarly,  we  can  factor  out  the  quantities  -  in  the  H,  (uu) 

\  \ 

(see  Eq.  (3.20)).  In  summary,  Eq.  (6.13)  may  be  written  in 
the  form 


106 


(6.14) 


*  J  l(*.  ) 

[$  (u)  fut  )]  =  [m]Oh  (uu  )^['R  - 

ai  1  2  H  1  \  2  1(0) 

(6x6)  (6xN)  (NxN)  (NXN)  (6xn) 


where 

[m]  =  Vi(0)[c3[-ViT<l  [  —  1 

1  mk 

Vx)  =  [n* '  "2+  2iCokV° 


oo  —it  (it  -U)  ) 

I(jji-^2)  =  J*  e  2  I(t)  dt 


(6.15) 

(6.16) 

(6.17) 


and  where  Q  and  Q  are  the  natural  frequency  and  damping  ratio 

iC  K. 

of  the  k^  booster  normal  mode. 

Equations  (6.14)  through  (6.17)  completely  specify  our 
general  assumptions  concerning  the  nature  of  the  booster  inter¬ 
face  accelerations.  In  the  following  section  we  present  some 
of  the  data  relevant  to  the  interface  accelerations  and  discuss 
the  validity  of  assumptions  A  through  C. 


6 . 4  Comparison  of  the  Idealization  With  Relevant  Data 

At  our  disposal  we  have  two  sets  of  data  concerning  the 
interface  accelerations. 

The  first  of  these  sets  is  composed  of  the  acceleration 

time  histories  themselves.  As  stated  in  section  6.1  these 

*  This  differs  from  the  form  of  Eq .  (3.20)  in  that  the 

reciprocal  of  the  component  of  the  modal  mass  matrix  has 

been  factored  out. 
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time  histories  were  obtained  from  the  calculated  responses 
of  analytical  models  of  the  booster  structures  to  excitations 
derived  from  booster  engine  chamber  pressure  measurements. 

Ihese  measurements  were  obtained  for  27  rocket  engine  test 
firings  during  the  booster  engine  cut-off.  Thus  we  have  27 
sets  of  acceleration  data,  each  set  consisting  of  the  6  rigid 
body  acceleration  components  at  the  booster  payload  interface. 
Ihe  coordinates  for  these  components  are  as  illustrated  in 
Fig.  1  of  section  1.2. 

In  Figures  13  through  18  we  present  these  acceleration 
time  histories.  In  each  figure  is  shown  a  superposition  of 
each  of  27  cases  of  one  acceleration  component.  It  can  be 
seen  that  the  mean  of  the  x  components  clearly  indicates  the 
thrust  decay  of  the  booster  engines.  Ihe  remaining  components 
have  approximately  zero  mean  values.  In  the  sections  that 
follow,  we  shall  be  mainly  interested  in  these  zero  mean 
acceleration  components.  Ihis  will  imply  no  loss  of  generality 
in  the  remaining  discussions  since  the  mean  of  the  x  component 
may  be  easily  calculated  and  subtracted  from  the  original  time 
histories  to  yield  a  set  of  time  histories  with  a  zero  mean 
value . 

Ihe  second  kind  of  data  obtained  consists  of  the  power 
spectral  density  estimates  of  the  acceleration  components. 

We  define  the  power  spectral  density  of  the  and  k^ 
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Fig.  16 
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ACCELERAT IONCG  S) 
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Fig.  17 
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ACCELERAT IGNCG  S) 


T I  ME  C  SECONDS ) 


Fig.  18 
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acceleration  components,  $  ,  (x)  as 


$  (x)  =  E[A.  (uc)A  (uu)  ] 

3  k 


(6.18) 


.th 


where  A  (x)  is  the  Fourier  transform  of  the  j  acceleration 
3 

component.  From  elementary  probability  theory,  an  estimate 
of  the  quantity  is  ^ 


27  * 

-q-  A  ( x ) A  (uu)  =  §.  .  (uu) 

27  k=l  3  ^ 


(6.19) 


where  $  .  (x)  is  the  estimated  value  of  the  quantity  defined 

~k  3 

by  equation  (6.18).  The  Fourier  transforms  of  all  of  the  6x27 
acceleration  time  histories  were  obtained  for  frequencies  in 
the  0  to  50  Hz.  range  and  the  six  power  spectral  densities 
$^(x)  (k=l,..0,6)  were  estimated  according  to  Eq.  (6.19). 

These  functions  are  presented  in  Figures  19  through  24.  It 
can  be  seen  that  the  non-zero  mean  value  of  the  x  component 
is  reflected  in  the  large  peak  at  x=0  on  the  spectral  density 
plot.  However,  the  x  component  power  spectral  density  still 
provides  useful  spectral  information  at  frequencies  higher 
than  ~  10  Hertz . 

There  is  a  simple  relationship  between  the  power  spectral 
densities  and  the  generalized  spectral  densities.  For  the 
acceleration  components  with  zero  mean;  the  diagonal  components 
of  [i  (x  , x  ) j  evaluated  at  x  =x  =x  are  just  the  power 


aT  1  2 


1  2 


spectral  densities  defined  by  Eq.  (6.18).  Hence  the  estimates 
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Fig.  2.0 
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ENERGY  DENS  MEANSC  THETA-XD 
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Fig, 
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of  the  power  spectral  densities  are  estimates  of  the  diagonal 


elements  of 


V’ 


These  estimates  are  only 


first  order  statistics  since  we  have  provided  no . calculation 
of  their  confidence  intervals.  A  detailed  calculation  of  the 
higher  order  statistics  of  the  spectral  density  measurements 
is  beyond  the  scope  of  this  work.  So,  for  the  present  dis¬ 
cussions  we  shall  assume  that  the  estimated  power  spectral 
densities  are  valid  representations  of  the  exact  power  spectra 
as  defined  by  Eq.  (6.18)* 

Having  presented  the  data  discussed  above,  we  shall  now 
discuss  the  validity  of  the  assumptions  stated  in  the  last 
section.  We  shall  examine  these  assumptions  one  by  one  in  the 
light  of  the  observed  characteristics  of  the  time  histories 
and  of  the  power  spectral  density  measurements. 


Assumption  A(i): 

The  acceleration  time  histories  clearly  exhibit  a  random 
nature.  Hence  if  these  accelerations  are  to  be  regarded  as 
the  response  of  a  linear  system  to  some  initial  excitations, 
we  must  conclude  that  the  excitations  initially  applied  to 
the  linear  system  are  themselves  random. 

Examining  the  acceleration  time  histories  more  closely 
we  see  that,  for  the  components  with  zero  mean,  the  number 
of  zero  crossings  is  roughly  equal  to  the  number  of  peaks. 
This  observation  corresponds  well  to  a  description  of  a 
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superposition  of  narrow  band  processes 


(2) 


Such  processes 


are  typical  of  the  response  of  lightly  damped  mechanical 
systems  to  broad-band  random  excitation. 

Hence  the  characteristics  of  the  time  histories  themselves 
suggest  the  validity  of  assumption  A(i).  This  assumption  is 
further  strengthened  by  the  fact  mentioned  in  section  6.1; 
that  the  individual  accelerations  were  obtained  by  calculating 
the  responses  of  an  analytical  model  of  the  booster  structure 
to  some  initially  applied  excitation.  We  can  then  conclude 
that  the  linear  system  mentioned  above  is  the  booster  structure 
itself.  Of  course,  in  the  27  cases,  different  booster  models 
were  used.  We  hypothesize  a  booster  structure  which  has,  in 
some  sense,  the  “average"  properties  of  all  the  booster 
structures  used  in  the  calculation  of  the  interface  accelera¬ 
tions.  In  the  following  discussions  we  shall  use  the  term 
"booster  structure"  to  denote  this  hypothesized  booster 
structure . 

From  the  foregoing  observations  we  conclude  that  the 
interface  accelerations  are  components  of  response  of  the 
booster  structure  (in  the  sense  described  above)  to  an 
initially  applied  random  excitation.  Thus,  equation  (6.7) 
is  indeed  appropriate. 
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Assumption  A(ii) : 

Equation  (6.7)  predicts  the  following  form  for  the  power 
spectral  densities  of  the  interface  accelerations: 

=  [C]^  (uu)^3[uu4i  (uj,ul)^H^(u;)  ^[C]T  (6.20) 

J  V  Q  V 


Or,  expanding  the  matrix  product: 
N  N 


’  m=ln=l  CjA  <“’)<“'4  ?0  (")Ckn 

j  in  ±n  -l  j  m  mn  n 


y|  A 

=  ^  (“>  *C  (w,ti))Q 

£=1  3  1  Jli  kl 

N  4 

+  E  E  C.  H.  ( uu )  (id  3n  (uj,uj))H*  (uu)C. 

_ -i  /  im  t>  v  t>  kn 

m— l  r\j£m  J  m  mn  n 


(6.21) 


where  (x  ,uu  )  is  the  generalized  cross-spectral  density 

mn  1  2 

of  the  mth  and  n components  of  [Q]  .  In  particular  for  j=k, 
we  have  the  spectral  densities  shown  in  Figs.  19  through  24: 


N 


(*)  =  E  c2  |r  (ug)  |  (u)4$  (oj , uu) ) 

kk  1=1  k l  '  QlJL 


N  N  *  4. 

+  m?1  ScmScn  Re(Hb  (Uj)Hb  $Q  (a,'W,) 

m-J-  n=l  m  un  mn 

n^m  (6.22) 

Now,  suppose  that  the  (x)  is  the  force  input,  displacement 

l 

th 

output  frequency  response  function  of  the  i  booster  normal 
mode.  That  is,  assume 

V  x)  in  the  form 


2  2  "1 
H,  (JU)  =  [n  -  «  +  2i  C  Q  uu] 

bA  k  Ok  k 


(6.23) 


124 


as  in  Eq.  (6.16).  With  C  «  1,  |h  ( jj )  |  peaks  sharply 

°k 

near  tx  =  Q  and  drops  off  to  negligible  values  for  |uu-Q  | 
k  k 

equal  to  several  multiples  of  the  bandwidth,  2  Q  Cl  .  Thus, 

Ok  k 

when  looking  at  the  power  spectral  density  measurements  for 

til 

frequencies  in  the  vicinity  of  the  l  modal  frequency  of  the 


booster,  we  should  see  a  peak  of  the  form 

V 

i 

Q, 


W*1  ~  CktlHb,(“,)|2(“’4  *Q,,(U‘'“')) 

*j  jL  Jo 


N 

E 


_  ScjlScn  Re(Hb  (i)Hb  $Q  (“»“*))  (6*24) 

n-1  U  0  VP.n 

n/m 


'l  ~n  £n 

th 


where  we  have  assumed  that  the  k  modal  frequency  is  several 

bandwidths  distant  from  any  other  modal  frequency. 

Turning  to  the  data  we  note  that  the  power  spectral 

density  estimates  exhibit  sharp,  distinct  peaks,  each  of  which 

closely  resembles  a  function  of  the  form  of  |h^  (id)  |  or  of 

the  form  Re  (H^  ('x)H^  (x))«  From  (6.24)  and  this  observation, 

°  i  n 

4  ~ 

we  conclude  that  the  quantities  (x  $  (x,x))  must  be  nearly 

mn 

constant  near  modal  frequencies  of  the  booster.  For  mathe- 

4  * 

matical  simplicity  we  assume  that  the  (x  (x,x))  (m,n  = 

ymn 

1 ,  .  o  . ,  N)  are  constants  for  all  x.  Hence  for  ^  : 


2  2  - 

(x-  X,  4  (X _  . ou_)  )  =  K  F(x  -  x  )/F(0) 

1  2  Qmn  1  2  mn  12 


(6.25) 


(m,  n  =  1 ,  .  .  . ,  N) 
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where  the  K  are  constants  and  F(x)  is  some  function  con- 
mn 

tinuous  at  0)  =  0.  But,  it  is  obvious  that  the  right  hand 
side  of  Eq.  (6.25)  is  just  the  cross-spectral  density  of  two 
components  of  a  general  vector  shot  noise  process. 

Thus,  we  conclude  that  with  the  (x)  given  by  (6.23) 

i 

the  properties  of  the  spectral  density  measurements  require 
that  the  random  excitation  to  the  idealized  structure  is  a 
broad  band  nonstationary  process.  For  convenience,  we  can 
assume,  to  a  good  degree  of  approximation,  that  the  random 
excitation  is  a  vector  shot  noise  process. 

Assumption  B 

It  must  be  noted  that  the  magnitudes  of  the  acceleration 
time  histories  exhibit  a  general  time  trend  which  is  slowly 
varying  compared  to  the  rapid  periodic  fluctuations.  The 
accelerations  with  zero  mean  rise  to  maximum  values  in  less 
than  one  second  and  rapidly  decay  to  zero  after  one  second 
(not  shown  in  Figs.  13  through  18).  This  general  time  trend 
we  associate  with  the  variation  of  the  intensity  functions  of 
the  shot  noise  process  which  excites  the  booster  structure. 

Moreover,  as  the  plots  of  the  time  histories  show,  the 
general  time  trends,  as  measured  by  functions  which  envelope 
the  acceleration  maxima,  are  of  the  same  general  form  and 
duration.  It  can  then  be  concluded  that  the  intensity  functions 
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of  the  initial  random  excitation  to  the  booster  are  approxi¬ 
mately  the  same  except  for  multiplicative  constants.  'Thus 
Eq.  (6.11)  is  approximately  valid* 

Assumption  C 

According  to  assumption  B  the  general  time  trends  of  the 
interface  accelerations  are  the  result  of  the  time  variation 
of  a  nonstationary  shot  noise  intensity  function.  If  this  is 
true  it  can  be  inferred  from  the  time  histories  that  this 
intensity  function  begins  at  zero  at  the  time  origin,  increases 
to  its  maximum  value  after  approximately  one  second  and  then 
drops  back  to  zero  a  few  seconds  later.  Thus  it  is  apparent 
that  the  magnitude  of  the  Fourier  transform  of  the  intensity 
function  assumes  a  maximum  value  near  a  frequency  of  the  order 
of  one  Hertz.  This  frequency  is  much  less  than  the  center 
frequencies  of  any  of  the  peaks  observed  in  the  power  spectra 
measurements  presented  above.  Assuming  the  magnitude  of  the 
transform  of  the  intensity  function  drops  off  rapidly  enough 
from  its  maximum  value,  we  conclude  that  the  generalized 
spectral  densities  of  the  interface  accelerations,  as  given 
by  Eq .  (6.11),  take  on  significant  values  only  when  (uu^-uu^l 

is  small. 

In  the  light  of  this  assertion,  examine  Eq .  (6.11). 

Expanding  the  matrix  product  we  have 
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$ 


N 

(u),  ,UU  )  =  E  C  C  K  H  <UL)  )H*  (uu  ) 
aT  1  2  D  &  b£  1  bx  2 

*kj  Z_1 


N 


+  E  E  C  C  K  H  (uu  )H"  (uu  ) 
m=l  n/m  jm  kn  mn  cm  1  bn  2 


I(U)1-UJ2) 


(6.26) 


Now,  from  the  spectral  density  measurements  we  see  that, 

in  almost  all  cases,  the  center  frequencies  of  the  peaks  are 

separated  by  several  band-widths.  Thus  the  value  of  |h  (uu  )| 

om  i 

is  quite  small  at  the  center  frequency  of  a  different  booster 

normal  mode.  Therefore,  with  luu  -uu  I  small  the  quantities 

1  2 

★ 

H  (x  )H  (x  )  (n/m)  are  much  smaller  than  the  quantities 

bm  1  bn  2 

H  (x  )H*  (x  ).  Now,  the  quantities  C  C  are  certainly 
hi  1  hi  2  j t  k l 

comparable  in  magnitude  to  the  quantities  C  C  (m/n) .  We 

jm  kn 

conclude  that  if  the  K  (m^n)  are  less  than  or  at  most  com- 

mn 

parable  to  the  then  the  second  term  within  the  braces  in 

Eq.  (6.26)  is  negligible  compared  to  the  first  term. 

Now  suppose  that  we  calculate  the  second  order  central 
moment  response  of  a  payload  structure  using  Eq .  (4.12),  with 

A 

r$  (x  ,x  )]  as  the  excitation.  The  second  moments  can  be 
al  1  2 

expressed  as  the  sum  of  two  terms.  One  of  the  terms  arises 
from  the  first  term  of  Eq .  (6.26)  and  the  other  term  arises 

from  the  second  term  of  Eq.  (6.26).  From  preceding  discussions 
it  can  be  concluded  that  the  second  group  of  terms  in  the 
second  moment  responses  can  be  neglected  in  comparison  with 
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the  first  group  of  terms.  That  is,  the  second  moments  of  the 


response  of  a  payload  to  excitations  given  by  Eq .  (6.26)  can 

be  approximated  by  the  response  to  excitations  characterized 
by 

N 

$  (’jj 

'Thus  if  the  diagonal  elements  of  [KmnJ  are  not  much  less  than 
the  off  diagonal  terms.  Then  approximately  correct  results 
may  be  obtained  with  the  use  of  the  approximate  spectral 
densities  given  by  Eq.  (6.27). 

Writing  Eq .  (6.27)  in  the  matrix  form 

(6.28) 

we  see  that  the  off  diagonal  elements  of  [ K ^ m 3  are  effectively 
zero.  This  assertion  is  equivalent  to  the  statement  of  assump¬ 
tion  C.  We  have  thus  demonstrated  in  some  measure  the  applica¬ 
bility  of  this  assumption. 

Having  established  our  basic  assumptions  about  the  nature 

of  the  interface  accelerations  we  devote  the  next  two  sections 

to  a  more  detailed  estimation  of  the  statistical  parameters. 

As  Eq.  (6.14)  shows,  we  have  yet  to  determine 

(a)  the  constants  Q  and  C  ,  describing  each  booster 

k  Ok 

normal  mode 


,U)  )  ~  £  C.  C  K  „H  ( jj  ) H*  (w  ) 

1  2  £=1  ]  AjI  W  b l  1  bf  2 


(6.27) 
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(b)  the  elements  of  the  constant  matrix  [m],  and 

(c)  the  intensity  function,  l(t). 

In  section  6.5  we  discuss  the  determination  of  the  power 
spectral  density  parameters,  (a)  and  (b),  and  in  section  6.6 
we  estimate  the  intensity  function,  (c) . 


6 . 5  Estimation  of  the  Spectral  Density  Parameters 

According  to  Eq.  (6.14)  the  power  spectra  presented  in 
Figs.  19  through  24  have  the  form 

N  2  .  .2 


* vk  (W)  =  £  M  |h  (ui)  I 

kj i  oi 


where  H^tw)  =  -  uj  +  2i£QK 


-1 


(6.29) 

(6.30) 


From  Eq.  (6.29)  and  (6.30)  we  can  show  that  from  the  shapes 
of  the  peaks  in  the  power  spectra  the  quantities  0^  and  £ 
(k=l,.o.,N)  can  be  determined,  and  from  the  magnitudes  of  the 
peaks  the  absolute  values  of  the  elements  of  [m]  can  be  found. 
But  first  we  must  identify  and  label  the  booster  structure 
normal  modes* 

In  Figure  25  is  shown  a  plot  of  the  number  of  the  more 
significant  spectral  density  peaks  centered  at  any  given 
frequency.  We  see  that  the  positions  of  the  peaks  cluster  at 
certain  locations.  We  assume  that  the  positions  of  the  peaks 
within  each  grouping  of  peaks  correspond  to  an  estimate  of 
one  modal  frequency  of  the  booster  structure.  These  groupings. 
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(Hertz ) 


6  in  number,  are  thus  associated  with  particular  booster 
normal  modes.  We  have  indicated  this  correspondence  by 
numbering  each  grouping  as  shown.  The  peaks  within  the  grouping 
labeled  "mode  5",  for  example  are,  assumed  to  represent  the 
response  of  the  5^  normal  mode  of  the  booster  to  the  initial 
shot  noise  excitation. 

With  the  booster  mode  identified  as  indicated  above,  we 
see  that  each  significant  peak  in  the  estimated  spectral 
densities  corresponds  to  one  of  the  frequency  response  functions 

of  a  booster  normal  coordinate.  For  a  peak  corresponding  to 

th  th 

the  k  booster  normal  mode  in  the  j  acceleration  component 

spectral  density  the  spectral  density  is  given  approximately 


by 


“  iHbku)| 


2  2 
M 


M 


jk 


jk  (n2-*2)2  +  4G2  nV 

k  Ok  k 


(6.31) 


This  is  valid  as  long  as  the  peak  in  question  is  not  within 

a  few  band  widths  of  another  peak.  Assuming  that  C,  is  small 

Ok 

-2v 

(of  order  10  )  we  may  write 


uu  =  ft.  _ 
MAX  kD 


(6.32) 


Ok 


2  D 


(half  height  width  of  i  (x)  )  (6.33) 

3  j 


MAX 


|M.,  ~  2  u  *  .  .  (t) 

3k  Ok  MAX  3  3  MAX 


(6.34) 


to  order  Q  ,  where  i  .  .  (id)  is  the  maximum  value  of 

Ok  3  3  MAX 
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§  .  .  (x)  and  n  =  fl  7l-C2  . 

33  kD  K  Ok 


*.  (w) 

DJ 


MAX 


$  .  .  ('JU 

DD  MAX 


) 


The  quantities  given  by  Eqs.  (6.32)  to  (6.34)  were  cal¬ 
culated  for  each  peak  in  each  of  the  spectral  density  estimates. 
The  resulting  numerical  values  are  presented  in  Table  I.  The 
calculations  for  the  x  or  q^  component  spectral  density  were 
also  included  since  the  large  peak  at  jj  =  0  has  little  effect 
on  the  size  and  shape  of  the  20  Hertz  peak.  However  a  cal¬ 
culation  of  the  M  associated  with  the  x  component  is  not 

jk 

included  in  Table  I.  As  an  example  of  the  details  of  the  cal¬ 
culation  consider  the  deter¬ 
mination  of  the  parameters 
associated  with  the  peak  near 
20  Hertz  on  the  power  spectral 
density  of  the  z  acceleration 
component.  The  position  and 
magnitude  of  the  maximum  and 
of  the  half  power  points  are 
indicated  by  Fig.  26.  With  the  values  of  these  as  given,  we 
proceed  as  follows: 


Fig.  26 


(a)  Q  =  JJ  =  20.33  Hz 
v  '  kD  MAX 

=  127.76 


radians 


sec . 


from  Eq .  (6.32). 
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TABLE  I 


Spectral  Density  Parameters 


Acc  ELE  RATION 
Component 

i 

Booster 

ttoDE 

K 

n„. 

(  E«.  (6.31)) 
(radians  / s tc.) 

(E$.  (6.32)) 

1  Mjk  1 

(Ea. (6.33)) 

1 

5 

1 29.848 

•  o43 

2 

1 

48.167 

.  1  05 

1  K„  I 

=  10.79 

n 

2 

80.634 

.  07  0 

|K«I 

=  5.11 

»! 

3 

93-201 

.049 

|Mm| 

=  7.35 

It 

5 

1 24.617 

.  03b 

|M„| 

=  10.58 

II 

6 

153- 93B 

.055 

1  M„| 

=  44.23 

3 

2 

79.583 

.076 

|M„| 

=  35.62 

ti 

5 

127.756 

.u43 

1  *'■  JS  1 

=  87.71 

4 

3 

96.340 

.  1  ^4 

1  *'-  t.}  1 

=  8.394 

it 

5 

126.711 

.046 

|K,SI 

=  0.512 

5 

2 

ol .682 

.  Oo9 

1  M«| 

=  0.190 

it 

1^3.673 

.061 

IKs-rl 

=  0.377 

»i 

5 

1 29.852 

.039 

lM«l 

=  0.201 

0 

1 

4b. 171 

.105 

|M..  1 

=  0.068 

it 

2 

80.634 

.071 

1  km  1 

=  0  •  07  4 

it 

5 

1 30.859 

.046 

1  Kts  | 

=  0.271 

it 

6 

1 51 .844 

.050 

1  1 

=  u.ol6 
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(b) 


(c) 


The  half  power  points  are  located  at  19.34  and  21.10 
Hertz.  Hence  the  half  height  width  is  (21.10- 
19.34)  Hertz  =  1.76  Hertz.  Therefore,  from  Eq. 
(6.33)  : 


1.76 


Ok 


2(20.33) 


=  0.043 


Finally  from  Eq.  (6.34) 


2  /  -  2 

| M35  |  ~  2  (0.043)  (127.76)  V0. 387x10 


=  87.71 


In  a  similar  way,  all  the  other  values  listed  in  Table  I 
were  obtained. 

As  would  be  expected,  the  center  frequencies  and  damping 
coefficients  of  any  one  booster  mode  fluctuate  randomly  from 
one  acceleration  component  spectral  to  another.  In  Table  II, 
the  average  values  of  these  center  frequencies  and  the  maximum 
values  of  the  damping  ratios  are  presented.  Table  III  presents 
the  matrix  [|m  |]  with  the  booster  modes  designated  as 

explained  above.  We  shall  use  the  contents  of  Tables  II  and 
III  for  the  example  at  the  end  of  the  chapter. 

It  should  be  noted  that  we  have  obtained  only  the  absolute 
values  of  the  components  of  [m].  The  signs  of  the  components 
can  be  determined  once  the  off  diagonal  terms  of  [  $  .  ('Jj)I  have 

been  estimated.  However,  this  involves  the  calculation  of  15 
quantities  similar  to  those  given  by  (6.18).  In  some  cases, 
this  great  mass  of  calculation  can  be  avoided 0  In  the  example. 
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TABLE  II 


Booster  Normal  Modes 


Booster 

Mode 

k 

Average  Modal 
Frequencies 

O-KP 

(  RADI  ANS  /SECOND  ) 

Max  1  MUM 
Damping 

Coefp l Cl  E  NTS 

(  Yok 

1 

46.169 

.105 

2 

80.633 

.  07  6 

3 

94.770 

.104 

4 

103.673 

.061 

5 

128.380 

.046 

6 

152.891 

.055 
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TABLE  El 

Elements  of  the 
Matrix  [IfiJ] 


MJ]  = 


Mode  No. 


i 

2 

3 

4 

5 

6 

T 

T 

10.79 

5.11 

7-35 

0.00 

1  o .  58 

44.23 

0.00 

35.65 

0.00 

C.00 

87.71 

0.00 

0.00 

0.00 

0.39 

0.00 

U.51 

0.00 

% 

0.00 

0.19 

0.00 

0.33 

u.20 

0.00 

5. 

0 . 68 

0.07 

0.00 

0.00 

0.27 

0.62 

137 


we  show  how  the  signs  of  the  components  of  [m]  can  be  obtained 
through  physical  considerations  alone. 

In  the  next  section  we  consider  the  estimation  of  the 
intensity  function,  I(t). 

6.6  Estimation  of  the  Intensity  Function 

We  here  assume  that  the  intensity  function,  I(t),  is 
slowly  varying  in  the  sense  described  in  section  5.6.  Now, 
we  have  idealized  the  interface  accelerations  themselves  as 
the  response  of  an  elastic  structure  to  uncorrelated  compon¬ 
ents  of  a  vector  shot  noise.  It  follows  that  the  response  of 
the  booster  to  one  uncorrelated  excitation  component  is  un¬ 
correlated  with  the  response  to  another  excitation  component. 

(i)  th 

Thus  if  K  is  the  variance  of  the  j  acceleration  due  to 

aj 

the  i  modal  component  of  the  vector  shot  noise,  then  the 
total  variance  is  just  the  sura 

N  (i) 

K  .  (t)  =  2  tcv  (t) 

al  i=l  aj 

However,  if  the  shot  noise  is  slowly  varying  then 

(i) 

Example  6  of  section  5.5  shows  that  K  (t)  is  approximately 

a  j 

proportional  to  the  intensity  function  of  the  shot  noise  com¬ 
ponent.  Since  we  have  assumed  these  intensity  functions  to 
be  the  same  within  multiplicative  constants  the  total  variance 
of  any  one  acceleration  component  will  be  proportional  to  the 
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intensity  function,  I(t).  That  is 


*  At)  «  i  (t ) 

*3 


(6.35) 


Now,  for  the  acceleration  components  with  a  zero  mean 

value,  the  Tchebychev  inequality  states  that  curves  along 

which  v'K  .  (t)  is  constant  are  those  for  which  the  probability 
a  1 

of  reaching  the  magnitude  of  acceleration  given  by  the  curve 
is  a  constant.  Hence,  if  we  construct  a  function  which  closely 
envelopes  all  the  ensemble  members  of  each  acceleration  com¬ 
ponent,  we  can  conclude  that  ic  (for  all  j)  is  proportional 

aj 

to  the  function  represented  by  the  curve.  A  suitable  function 
for  this  purpose  is 


C(t)  = 


A(l-cos  Gt) 
H 

0 


0  <  t  < 
0  >  t  > 


2tt 

G 

2n 

G 


(6.36) 


G  was  evaluated  by  noting  the  time  position  of  the  maxima 
attained  by  each  time  history,  averaging  these  values  and 
dividing  2~  by  the  resulting  average.  The  approximate  value 
obtained  in  this  way  is 

e  =  4.833  radians  (6.37) 

sec . 

With  this  value  for  e,  the  curve  given  by  (6.36)  was  found 
to  be  a  good  estimate  of  the  enveloping  function. 

Since  *Jk  06  C(t)  (all  j),  relation  (6.35)  gives 

aj 
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I(t)  cc  C  (t) 


or 


I(t)  = 


C  (1-cos  et) 
0 


0  <  t  <  — 


0  >  t  >  -^TT 


(6.38) 


where  C  is  an  undetermined  constant. 

We  shall  need  the  transform  of  I(t)  for  the  response  for¬ 
mulation.  A  straightforward  evaluation  of  the  Fourier  integral 
by  elementary  methods  yields 


-IX 


2tt 


i  (e 


e  -1)  6  e4 


I  (uu)  =  C 


2  2  2  2 
x(x  -e  )  (ol3  -4e  ) 


and  since 


we  have 


2tt/ e 


I ( 0 )  =  [  C (1-  cos  et)  dt  =  —  C 

"o  G 


2rr 

-ix —  ^  c 

•  /  e  .  2  5 

; ,  .  i (e  -1)  ~  e 

I  (uu)  _  1 

^  “  - 2 - 2 - 5 - 7 

1(0)  x(x  -  e  )  (x  -  4e^) 


(6.39) 


From  (6.39),  it  is  apparent  that  I(t)  is  indeed  slowly 
varying  since  the  major  components  of  its  Fourier  spectrum 
are  centered  at  ~  0.77  Hertz  while  the  booster  frequency 
response  functions  are  centered  in  the  region  7  to  35  Hertz. 


Having  obtained  a  crude  but  serviceable  model  of  the  inter¬ 
face  accelerations,  we  now  turn  to  the  problem  of  calculating 
the  response  of  a  payload  to  these  excitations. 
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6.7 


Response  Formulation 


'The  two  fundamental  relations  are 

2 

{V,}  +  2[>Q  +  Ox.^^}  =  [Cl  [a]  (6.1) 

11  1 

(MXl)  (MXM)  (MXM)  (MX  6)  (6x1) 

and 

,  *  T  hlD.-X  ) 

C*a(u)1#a»2)3  =  [M]|^Hbi{au1Klt-Hbj  (a)2fl[M  ]  — - - (6.14) 

I  (a ) 

According  to  Eq «  (4.12),  the  Fourier  transforms  of  the  second 

order  central  moments  of  [p}  are  given  by 

J  (y,y-x)H_.  (x-y)dy  (6.40) 

In  this  case: 


[$«  (y,y-x) }  =  [c]  [M]  t>H,.(y)H  (y-x)^[M]  [c] 

y^-j  bi  bi 


T  I  (x) 

1(0) 


(MX 6) (6xN) 


(NXN) 


(NX  6 )  (6xM) 


so  that 

t  \  N  6  6 

%  (y,y-x)  =  E  E  E  c  m 


1(0)  i=l  m=l  n=l  ml 


M  C  (H (y)H  (y-x)) 
i  ni  kn  bi  bi 


(6.41) 


Substitution  of  (6.41)  into  (6.40)  yields 


N 


K  (x)  =  S  E  E  C  .M  .M  .  c. 
k3  ,  ,  m]  mi  ni  Jen 


i=l  m=l  n=l 


m3  mi  ni 


(6.42) 


1  2  2  -1 
where  H  ( x )  =  —  (X.  -  x  +  2iC,  x,  x) 

k  m^  k  k  k 
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Tbe  integral  on  the  right  in  Equation  (6.42)  is  nothing 
but  the  integral  of  Eq .  (5.69)  of  section  5.6.  Hence  the 

inversion  of  the  Fourier  integral  gives 


N  6  6  t  (i) 

K  (t)  =  S  Tj  S  C  M  M  C  [ -1  ^  y  ("  )dT  ] 

i=l  m=l  n=l  mi  ni  kn  0  i(0)  *3 

(6.43) 

(i) 


The  quantity  y  (t)  is  given  by  Eq.  (5.81)  with  0  and  Q 
replaced  by  and  respectively: 

(i) 

Vkj  (t)  =  A[Bi  6  3  3  COs(\D'J'jD)t 

-<Ck'\+£oini,t  ,  „ 

-Be  cos  u;  -0  )t 

2  kD  jD 

-(C  x.+£  .n.  )t  -2:  .n.t 

-Be  3  ^  Oi  i  cos  -Q  )t  +  B  e  O3-  i  ^ 

3  jD  D  4 


2  2  2  2  -1 

A  =  [8  r  a  0 ,  (Q  -x  ) (ft  -x  ) j 
kD  jD  iD  iD  jD  iD  kD 


(6.44) 


B  =  2Q .  _  +  x„  +  x . 

1  iD  kD  jD 

x  .  _ 


B  =  pp —  {Cl .  +  xVn)  +  (0  +  x  ) 

2  Cl  .  iD  *D  iD  j D 

lD 


X. 


B  = 


kD 


-  (n.  +  w .  )  +  (0.  +  o),  ) 

3  Qi  D  lD  lD  kD 


X  .  X 

B  =  (x  +  0.  )  +  (x.  +  Cl  ) 

4  aiD  kD  id  niD  dd  d 


y 


where 


w  =  ft../l-C . 

iD  1  Oi 
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X  =  X  V-L-C 
kD  k  k 


=  X  Jl-L 

\ 

!  2 

x  =  x  Vl-C 
jD  j  j 


(6.44)  yields  a  result  for  the  convolution  integral, 

^  (T)  3t)  which  is  correct  except  for  the  term 

0  1(0) 

„t  ,  /+.  (i) 

(T1+T2+T3)  -4-T'~ — L-  .  (T)  d~)  whore,  from  (III  .15)  and 

section  5.6: 
o  {T  }  < 


O  l 


2 [ MAX ( |  *  |  )  JP+1[MAX(  *  \  )  ^MAXp  J  (a^+)  )  |  ] 

(P+2 )  [  MIN  ( |  a^  +  ha  f +  )  |  )  ][MIN(|a(  +  )-a(_)  |  )  ]P+2 }  MIN[  |  J  (a  (_) )  |  ] 
k  3  k  3  k 

(-) 


} 


2MAX(|cO  1 1  ) 

o  {T  }  <  O  { - — -  ) 

2  MIN ( | a  +  |  ) 


O  [T^]  <  O  {  MAX  (Coi«  Ck»  C.)l 


(6.45) 


It  is  understood  that  the  locations  of  the  poles  of 
frequency  response  functions  are  denoted  by  a  •  In  formulae 

K 

( + ) 

(6.45),  a  denotes  all  combinations  of  two  a  which  involve 

k  j 

(-) 

sums  of  two  of  the  frequencies  x,  ,  x.  and  Q.  while  a 

kD  JD  ID  k 

denotes  all  such  combinations  involving  differences  between 
two  of  these  frequencies. 
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Equations  (6.43),  (6.44)  and  (6.45)  summarize  the 

approximate  formulation  of  the  second  moment  response  of  a 
payload  structure  to  the  interface  accelerations  as  idealized 
by  Eq «  (6 .14) . 

Once  the  approximate  second  moments  have  been  determined, 
the  member  stresses  can  be  found  from  a  simple  matrix  trans¬ 
formation  (e.g.  Eq«  (6.3)). 

To  illustrate  the  application  of  Eqs.  (6.43)  through 
(6.45)  we  present  the  following  example. 


Example  6.1 

A  two  degree  of  freedom  structure  is  driven  by  a  combina¬ 
tion  of  the  z  and  Gy  components  of  the  booster/payload  inter¬ 
face  accelerations.  Suppose  that  the  normal  equations  of 


motion  are 


/ ..  > 

- 

r  e  N 

r  2  i 

X 

-  1 

■  +  2; 

03  0 

1 

4 

X1 

►  + 

x  0 

1 

2 

j 

K1 

k  x 

0  a 

X 

0  x 

1 

2 

l  2  J 

V  2- 

L  2  J 

2, 

C  03  C  03 

11  2  1 

2  2 
C  03  CO) 

L  1  2  2  2 


r  •• 

z 

0 

>1 

(6.46) 


where 


x  ,  x 
1  2 


=  payload  normal  coordinates 


0  0  025 

16.5  Hz, 

20.5  Hz, 

1.0 

217  .4 


(6.47) 
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and  were  chosen  so  that  the  magnitude  of  the  peaks 
centered  near  16.5  Hertz  on  the  9^  spectral  density  equals 
the  magnitude  of  the  peak  centered  near  20  Hertz  on  the  z 
spectral  density.  A  comparison  of  Eq.  (6.46)  with  (6.1)  shows 


that 


[c]  = 


2  2 

0  0  C  x  OCX  0 


1  1 


2  1 


2  2 

0  0  C  'X_  0  C  x  0 


1  2 


2  2 


And,  from  Table  III: 
[C][M] 


0  M  C  X  +M  _C  X 
32  1  1  52  2  1 

2  2 
0  M_ _C  x  +M_  C  _x 
32  i  2  52  2  2 


2  2  2 

0  C  x  M  C  X  M_  +C  x  Mr_ 

2  1  54  1  1  35  2  1  55 

2  2  2 

0  C  x  M_.  C  xM+CX  M 
2  2  54  1  2  35  2  2  55 


Let 


2  r ^  .  (y)H,.  (y-x)-<l 

3  i(0)  ° 


Then  Eq .  (6.13)  becomes 


[$  (y,y-x)l  =  [ m 3  I>h.O  [m] 
a  D 


so  that 


t$Q(y.y-x  )3  -  [ cj r m! [m j t[ c j t 


o 

o 

(6.48) 


Or, 
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[$Q(y»y--u)  3 


2  2  2 
[H  (M„C  +M  C  )  +  H  (C  M  )  +  H  (C  M  +C  Mcc)  ] 
2  32  1  52  2  4  2  54  5  1  35  2  55' 


(6.49) 


4 

2  2 

■ju 

UU  JU 

1 

1  2 

2  2 

4 

U) 

UU 

1  2 

2  J 

To  proceed  further,  we  must  ascertain  the  signs  of  the 

M's  in  Eqe  (6.49).  Let  us  focus  our  attention  on  the  rela- 

2 

tive  signs  of  M^  and  since  the  term  H^(C^M35  +  C2M55^ 

will  most  affect  the  results. 

First,  we  notice,  after  a  consideration  of  Eq.  (6.14) 

••  •• 

that  the  Fourier  components  of  z  and  6^  in  the  vicinity  of 
booster  mode  5  are  proportional  to  M  and  respectively. 

It  can  be  seen  from  Fig.  27 
that  the  accelerations  along 
the  z  and  9  directions  in 

y 

the  normal  mode  centered  at 

—  20  Hz.  can  be  considered  a 

mode  of  vibration  of  a  beam 

in  the  x-z  plane.  If  we 

assume  this  to  be  the  case, 

then  it  can  be  concluded  that 

o 

z  and  0  at  20  Hz.  are  180 

y 

out  of  phase.  Hence  their 

Fourier  components  are  opposite  in  sign;  i.e.,  M  and  M^r 

j  j  o  J 
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have  opposite  signs. 


Assuming  the  above  argument  to  be  correct  and  referring 
to  Table  III,  we  can  set 


M 

=  + 

35.65 

32 

M 

=  + 

0.19 

52 

=  + 

0.38 

54 

M 

=  + 

87.71 

35 

M55 

_  T 

0.20 

Then,  with  and  given  by  (6.47): 


(M_  _ C  +M  C  )  =  5922.84 

32  l  52  2 

2 

(CM)  =  6824.41 
2  54 

2 

(M  C  +M  C  )  =  1956.29 

35  l  55  2 


(6.50) 


Finally,  the  result  corresponding  to  Eq.  (6.43)  is 


K 

kj 


(t) 


2  2 
x,  x  . 

k  3 


t 

[5922.84  J  I(t-T) 

“o 


("■)<3t 


t 

+  6824.41  J  I(t-") 
'  0 


(4) 

Y  (t  )  dT 

kj 


(6.51) 


t 

+  1956.29  ■  I(t-~) 

v  0 


(  T  )  d  T  3 


1 

1(0) 


I(t)  =  _e  (1 

1(0) 


2  ^ 

cos  e  t) 

► 


e  =  4.833 


(6.52) 
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Fig.  26 
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T I  ME  C  SECONDS  ) 


where  the  are  given  by  Eqs.  (6.44). 

kj 

With  the  numerical  values  given  by  Table  II  and  Eq.  (6.47) , 
Eq.  (6.51)  may  be  integrated  numerically  to  give  the  desired 
second  order  moments  of  the  system  response. 

The  integrals  of  (6.51)  were  evaluated  by  the  trapazoi- 
dal  rule.  Fig.  28  shows  the  resulting  time  histories  of  K  , 


*12  and  K22‘  In  Fi9s*  29(a)  and  (b)  3^*^  and  3V*22  (i.e., 

3o  and  3a  )  were  plotted  along  with  a  superposition  of  the 
1  ^ 

system  responses  to  each  of  the  27  sets  of  acceleration.  These 
were  obtained  by  a  numerical  integration  of  the  deterministic 


equations  of  motion.  It  can  be  seen  that  3V*  and  3*Jk 

11 

bound  most  of  the  deterministic  solutions  to  the  response. 

This  is  to  be  expected  since  by  the  Tchebycheff  inequality, 

Eq.  (2.28),  there  is  a  probability  of  approximately  0.89  that 
at  any  particular  time,  the  responses  are  less  than  three 
times  the  corresponding  standard  deviations. 

It  is  instructive  to  calculate  the  relative  errors  induced 
by  the  use  of  Eq.  (6.51) .  That  is  we  wish  to  calculate  the 
orders  of  magnitude  of  the  error  terms  associated  with  each 
of  the  three  terms  of  (6.51)  by  Eq.  (6.45).  As  an  example, 
we  shall  calculate  the  error  terms  of  K^^(t). 

Since  the  first  mode  of  the  system  is  tuned  to  the  center 
frequency  of  the  4th  booster  mode  the  second  term  in  Eq. 

(6.51)  gives  the  largest  contribution  to  *^ (t) .  Hence,  the 
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relative  errors 


associated  with  the  quantity  J 

0  1(0) 


(4) 


y  (')d^  will  be  calculated.  Ihe  positions  of  the  poles 


of  H^(x)  and 


are 


a  = 
1 

i 

Gu^  + 

tl) 

Id 

a> 

NJ 

II 

i 

Qx  — 

IX 

ID 

a  = 
3 

i 

G  n 

04  4 

+  n 

4D 

a  = 
4 

i 

G  n 

04  4 

-  0 

4D 

And,  we  then  have 

a^  +  *  =  ai+a-3  =  1  (C'x  +G  n  )  +  (su  +n  ) 
1  1  3  1  04  4  ID  4D 


(  +  ) 

a  =  2Q 
2 


4D 


=  -20 

4  4D 

(  +  ) 

a5  =  a2+a4  =  i(C'Xl+C04n4)-(XlD+r*4D) 
and  a^_)  =  i(C^  +  Q4V  +  (»1D“  ^4D> 

(-) 

a  =  i('x  +  0  )  +  (0  -x  ) 

2  1  04  4'  4D  ID 


(6.53) 


Hence : 


MAX(|a^  }|)  ~  Ct»1+C04n4  MIN(|a][+)-  a_j  +  )  |  ) 

MAX ( | a^+) -a j  }|)  -  MIN(|a^+)-a^-) | )  ~  2f>4D 


MAX(!a^  } !)  ~  Cx1  +  C04n 


(equation  continued ) 
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(6.54) 


.  <  +  ), 

MIN(|ak  |)  ~  2fl4D 


Now,  consider  the  Fourier  transform  of 


I(t) 

1(0) 


/v  -i'JU-^5- 

I  (J  )  _  i  (e _ -1) 


2  .5 


U) 


(uu2-e2)  (uu2-4e2) 


1(0) 

Since  the  first  factor  has  no  singularities  in  the  complex  uu 
plane,  we  set 


J  (x )  - 


•  2  TT 

i(e'1U)~  -1) 


(6.55) 


The  remaining  factor  has  4  poles  so  that 

P  =  4 

And  from  (6.53)  and  (6.55),  we  have 


(6.56) 


r  I  /  (  +  )  \  I  1 

MAX  [  |j(ak  )  |  ] 


i(e  1  c  -1) 


T+T 


MIN  L | J (a^  J ) | ] 


e(;'xl+C04n4)2e 
(  =  h  +  *04R4> 


(->2n 

i  (e  1  e  -1) 


TT 


e(^Xl+:04  ‘4} 


2n 
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Hence 


MAX  [ 1 J(a^+)) H 
MIN  [  |  J(a][~)  )  |  ] 


(6.57) 


Collecting  the  results  (6.54),  (6.56)  and  (6.57)  and  substi¬ 

tuting  them  into  (6.45)  we  obtain 

4 

(Cu),+C  OJ 

O  (Tl  <0(1  ..  1  °4  4  ) 

1  3 

**  0  (  I  c4  ) 

3  04 


2(Cx  +C  0  ) 

0  (T  )  <  0  {  - - - 04  4  }  «  O  (2  C  ) 


20 


04 


4D 


o  (t3)  <  o  {  max(G,C04)}  =  O  {C04l 


Or,  the  order  of  magnitude  of  the  total  error  is 


0  ly  y  t3)  <  o  (3C04) 


O  [0ol83] 


(6.58) 


Hence,  with  as  given  by  (6.51)  the  exact  variance  of  the 

first  mode  is  K^(l+T)  where  0{t]  <  0  {0.183}.  Hie  upper 
bound  to  the  variance  is  then 

S.U.P.tK  3  ~  (1*2) 

in  a  similar  way  the  errors  due  to  the  approximate 
formulation  can  be  calculated  for  the  other  moments  and  upper 
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bound  estimates  can  be 
since  the  calculations 
the  application  of  the 


made.  This  will  not  be  done  in  full 
already  made  are  sufficient  to  illustrate 
approximate  techniques  of  Chapter  V. 
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CHAPTER  VII 


SUMMARY  AND  CONCLUSIONS 

Here  we  summarize  the  basic  results  of  this  work.  We 
have  presented  an  approximate  evaluation  of  the  second  central 
moment  response  of  a  lumped  parameter,  second  order,  linear 
system  to  a  class  of  random  excitations.  This  class  of  ran¬ 
dom  processes  we  have  called  "slowly  varying"  nonstationary 
processes,  since  the  nonstationary  variation  of  the  correla¬ 
tion  functions  of  these  processes  is  small  compared  to  the 
variation  of  the  impulse  response  functions  of  the  dynamic 
system  under  consideration.  When  the  excitation  falls  under 
the  above  description  we  have  demonstrated  that  the  resulting 
response  of  a  structural  system  may  be  conveniently  approxi¬ 
mated  by  relatively  simple  analytical  expressions. 

These  analytical  simplications  make  possible  the  esti¬ 
mation  of  response  statistics  from  various  statistics  of  the 
excitations.  We  have  illustrated  the  possibility  by  a  con¬ 
sideration  of  the  problem  of  the  estimation  of  flight  loads 
in  a  payload  during  booster  engine  cutoff. 

The  scheme  of  analysis  outlined  in  Chapter  VI  seems  to 
offer  the  following  advantages: 
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1.  Integration  of  the  equations  of  motion  of  a  pay- 
load  system  for  27  forcing  functions  to  obtain  the  design 
loads  is  avoided.  The  scheme  proposed  in  Chapter  VI,  in  com¬ 
parison,  involves  only  marginal  computation  time. 

2.  Under  most  of  the  present  schemes  the  analysis  is 
highly  dependent  on  a  prescribed  engineering  policy,  i.e.,  the 
designer  must  first  decide  which  of  the  27  test  data  are  to 

be  used  for  the  structural  qualification  of  a  proposed  design. 
This  practice  can  produce  unreasonably  conservative  demands 
on  the  structural  design.  Under  the  proposed  scheme,  however, 
design  policies  can  be  formulated  from  the  analysis  itself. 
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APPENDIX  I 


Here  we  show  how  Eqs .  (3.55)  through  (3.57)  were  obtained, 

From  the  theory  of  the  inversion  of  Fourier  integrals 
I(t)  is  the  sum  of  the  residues  of  the  integrand  in  the  upper 
half  plane  multiplied  by  2ni.  From  Eq.  (3.40)  and  assumptions 
a.  and  b.  of  section  3.6  it  follows  that 


n  i  J(ak>  i\t 

I(t)  =  i[  E  Res  (— ; a  )  -  e 

F  k 


k=l 


G  (a  ) 
k 


m 


+  2  Res  (—  ;b  )  - —  *  ^ 

k=l  G  k 


J  (b  )  lb  t 

k  -i 

e  J 

F  (b  ) 
k 


(1.1) 


Now,  since  a  and  b  are  in  the  upper  half  plane,  their 
k  k 

imaginary  parts  are  positive  and 

ia  t  ib  t 

|e  k  |  <  1  ,  |e  k  |  <  1 

for  all  k. 

There  fore : 


I  v>  1  .  J(bk)  lbkt !  „  & 

j  Res  (—  ;b  )  -  e  |  <  2_> 


m 

2 

k=l 


F(b  ) 
k 


i  J(V 

Res  (—  ;b  )  - — 

k=l  G  k  F(b  ) 


,  ^  MAX[ | J (b  ) | ] 

<  m  *  MAX[|Res(—  ; b  )  | ]  - 

G  k  MIN[  | F (b  )  |  ] 


(1.2) 


But,  from  (3.53)  and  (3„54)  we  see  that 
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MAX  [  |  Re  s  (—  ,  b  )  |  3  < 

G  k  ^m-1 


(1.3) 


MIN  [  |  F(b  )  |  ]  >  y11 
*  1 

Hence, 

m  ,  J  (b  )  ib  t 

|  E  Res  (—  ;b  )  - -  k 

k=l  G  k  F (b  )  e  y 

k  2  1 


(1.4) 


e  |  <  "  n  MAX  [  |j(b  )  |]  (I -5) 


The  expression  on  the  right  is  an  estimate  of  the  upper 
bound  of  that  part  of  the  integral  derived  from  the  residues 
of  -  . 


Similarly : 

n  1 
E  Re  s  (—  ;  a,  ) 


J (a  )  ia  t  n  ,  J  (a  ) 

>  ,-v,  - e  k  |  <  E  |  Res  (—  ; a  )  - *_ 

k=l  F  k  G(a  )  k=l  F  k  G(a  ) 


MIN[ | J  (a  )  l  ] 

n  •  MIN  [ | Res (—  ,a  )|]  - 

F  k  MAX[ | G(a  ) | ] 


(1.6) 


From  (3.52)  and  (3.54)  : 

MIN  [  |  Res  (—  ;a  )  |  ]  >  —V 
F  k 


m 


MAX  [ |G(a  ) |3  < 

k  2 


(1.7) 

(1.8) 


Therefore : 

n 


0  U  E  Res  (—  ja^ 


k=l 


J  (a  )  ia  t 
— S-  e  k  |) 

G  (a  ) 
k 


o  [ — MIN[  I  J(a  )  I  1 } 

n-l%m  k 

1  2  (1.9) 
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With  (1.5)  and  (1.9) ,  we  can  rewrite  (1.1)  as 


1  JK> 

I  (t )  =  [i  S  Res  (~  ; a  )  - —  e  ]  (1+T  ) 

F  k  G(a  )  1 

k 


n 

s 

k=l 


(1.10) 


where  T  is  a  collection  of  terms  such  that 
1 

n-1  m 

m  6.  Y,  MAX[  |  J(b  )  |  ] 

0  1 1 T.  U  <  t  — - - —  1 

n  e  y"  MIN[  |  J(a  )  |  ] 

2  i  k 


(1. 11) 


Furthermore,  suppose  that  a.  =  a+A^,  where  z  =  a  is  some 
point  in  the  vicinity  of  the  a  (k=l ,  .  . . , n)  .  Then,  we  may 

K 

write : 

...  k  ...  =  — L-  -  [— i-  -  - I -  ] 


G (a  . )  G(a)  G(a)  G(a+A  ) 
3  j 


where  we  note  that  1^.1  <  6  .  Hence; 

3  1 


(1  +  T  ) 


(1.12) 


G(a^)  G (a ) 


o  [t  }  <  o  [max  [■ 


G(a) 


G(a +^) 


-  1]  } 


(1.13) 


and. 


n  i  J  ( a,  )  i  a^  t  n  i  a.  t 


v  i\  .=  * - —  ~  k  )  =  — —  (  f.  Res  (—  ;  a  )J(a  )e  k  ) 


i  (  Res  (-*  ;a  ) 
k=l  F  k 


G(ak) 


g(a)  k+1 


(1+T2) 


F  k 


(1.14) 


Substituting  the  right  hand  side  of  Eq.  (1.14)  into 
that  of  (1. 10)  we  have 
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n  -A- Cl  L- 

I  (t)  =  — —  (i  12  Res  (—  ;a  )  J(a  )e  ]  (1+T  )  (1+T  ) 

G(a)  k=l  F  k  k  12 


But  the  term  in  braces  is  nothing  but  the  integral 
1  n°°  J  (uu )  i&t 

—  J  -  e  dx.  Hence  we  may  write  I  ( t )  as 


-oo  F  (uu ) 


Kt)  =  [- 


J 


J(X)  iut 


where 


2nG(a)  F(i) 

n-1  m 


d-X  }  (1+T  )  (1+T  ) 

1  2 


0  [T  }  =  0  [  , 

1  m-1  n 


m  6^  y2  MAX[ |j(bk) I  3 


n  e2  MlN[|j(ak)|] 


0[T  }  <  0  [MAX  - i]  } 

2  G(a+A . ) 

3 

which  was  to  be  shown. 


(1-15) 


(1.16) 


(1.17) 


162 


APPENDIX  II 


In  section  4.4  we  have  presented  two  alternative  expres¬ 
sions  for  the  variances  of  response: 


Kkk(A)  = 


)  =  r:  j  h  (y)  $  (y,y-x)H  (x-y) 
^  -c»  k  ukk  K 


(II. 1) 


i  r  /•*> 


Kkk(x)  = 


1^2  -2Ck\i+y)Hk(y)^Qkk(y,y-x)dy 
- - - 

\<X-2C*V,(-  4  +iCkV  +  V 


(II. 2) 


Here,  we  show  that  these  two  results  are  equivalent. 

First,  consider  (II. 1),  and  make  the  transformation 

t  =  -  t2  ,  t*  =  t  (II. 3) 

Then,  by  (3.23);  and  since  the  Jacobian  is  one: 


~  1  ..  *  -i  (Ty+t  uj)  * 

Kkk(x)  =  T-  .  H.(y)H  (ii’-y)jj  CQ  (T,t  ) e  d-dt  dy 

KK  2  ~  K  k  —a*  vkk 


1  -iTy  ,  *  -it*'u  * 

=  —  :  I  H  (y )  H  (x-y)e  dy}c  (",t  )e  Xdt  d" 

O—  w  **  v  U.  - 


2-  ico  k  k 


kk 


-i"  t 


[J  a  (y)H  (x-y)e  dy]c  (T,t  ) 
2-  loo  ‘0  I.  K  k  Qkk 


*  -it*x, 

1  e  dT 


+  j*  Lc’  a  (y)H  (x-y)e-1  ydy'cQ  ("/t*)e“lt  Xd7]  dt 
-°°  loo  k  k  wkk 


(II. 4) 


where  we  have  reversed  the  order  of  integration  and  have 

expanded  the  integral  over  t .  In  the  second  integral  in  t 

* 

we  let  7  =  so  that 
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,  * 

-it  ii 


dr 


+  j"  I  C  (t  *,  t*) dT*  }dt* 
0  z  vkk 


(II. 5) 


where 


and 


I.  =  j  H  (y)H  (x-y)e  1  ‘  Ydy 

1  _  co  k  k 


I  -  f  H  (y)H  (ai-y)e1  Ydy 
2  1»  k  k 


Now,  the  coefficients  of  the  exponentials  in  both  of  the  above 
integrals  can  be  written 
^  (y)Hk('i-y)  = 


(y-a  ) (y-a2) (x-y-a^ ) (x-y-a2) 


(II- 6) 


with 


ai  lask*k  +  XkD  '  a2  1XkCk  “  \d 


so  that  I.  and  I  are  to  be  calculated  by  the  methods  of 
1  2 

section  3.5.  According  to  Ihe.'s  1  and  2  of  that  section, 

I  and  1^  are  integrated  over  the  upper  and  the  lower  half 
z  planes  respectively.  Performing  the  contour  integration, 
we  have 


'iAkD  r 

I  =  -  L  - 

1  x-a1-a2 


x-2a 


-i- (x-a^) 

£ - —  1 

x-2a 


.  y  iT *an  iT*a0 

TT  3  /X,  1  2 

_kD  r  e  e 


2  x — a  —  a  x  —  2  a 
1  2  1 


x-2a 


(II. 7) 
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Substitution  of  (II. 7) into  (II. 5)  yields 


kk 


1  _  oo  m 

(a')  =  J  (J 


TTl 


2tt 


-iT  (x-a  ) 
1 


0  x  (x-a  -a  )  (t'_?a 

U  kD  12  (Wa1 


-iT  (x-a  ) 

^ - -  ] 

uu-2a~ 


c  (t  / 1  ) e 
ukk 


-it*x 


dT 


TTl 


iT*a1  it*a2 


-  ^ - 1 

°  \D(X_ara2)  U,-2ai  ,JJ-2a2 


CQ  (-*,t*)e  lt*Xd’r*  }  dt* 

\k 


(II. 8) 


But,  by  (1.3): 

~  -iT(x-ai) 

J7  e  Cq  (T,t*)e  d-  dt* 

kk 

00  ia  (t  -t  )  -i  (t  -t  )uo  -it  id 

=  n  e  1  2  e  2  e  2  CQkk<Vt2,dVt2 

00  ia  (t  -t  )  it  ub 

=  rr  e  e  1  R  (t  ,t  )dt  dt 

J  „  x  1  2  1  2 


If  we  now  let  (t  -t  )  =  :  as  before  but  set  t^=t  # 

Cn  (t  ,t  )  is  the  same  as  before  since,  by  its  definition, 
ukk  1  2 

C~  (t  ,t  )  is  symmetrical  in  t  and  t  .  Hence: 

°kk  1  2  12 

“  -iT^"al>  — it  *x 

e  e  C_  (",t*)dTdt* 

O  O 

vkk 


lalT  -ixt* 


c 

—  CD 


cD  (T,t*)d-dt* 
wkk 


(equation  continued) 
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“  “iT(oi-a  ) 

I  I  e  e  C  (t , t*)dTdt* 

a  °kk 

°°  ia?T  _•  ■*-* 

=  J'J  e  e  1U)'*  C  (T,t*)dTdt* 

-oo  Kk 


(II. 9) 


Employing  relations  (1.9)  and  setting  t  =  in  Eq. 
(1.8)  we  have 


00  i/x.  ia  i  iaT 

k  <x)  =  i  ;;  — —  r- - - —  I  c  (T.t*) 

kk  2  x-a, -a_  x-2a,  x-2a  ^kk 


1  2  1 
-ixt  * 


dTdt  * 


(II. 10) 


Now  we  consider  Eq.  (1.2)  which  can  be  written  in  the 


form 


K  (x)  =  - AtLi - 

kk  (x-a  -a  )  (x-2a  )  (x-2a  ) 

12  1  2 


A(x)  =  i-  J  4(-  ^  +  a  +a  -y)H  (y)  £  (y,y-x)dy 

2Tr  L-  2  12k  ukk 


(II. 11) 


where  a^  and  a^  are  given  by  (II. 6).  We  can  express  "(ix)  as 


:.(x)  = 


4 (-  %  +  a  +a  -y) H  (y)$.  (y.y-x)dy 


2tt  la,  2  1  2  J  '  *k  w  #  Q 


kk 


=  2n*  1"  C  4(-  |  +  ai+a2-y)\(y)e_;LTydy  } 


_  OO  —  oo 


—  ixt  * 

C  ("  ,  t * )  e  d"dt* 

°kk 


=  ~  £’  I,’  C_  (T,t*)dT  +  j*  l‘  C„  ("*,t*)dT*  } 


2”  la,  •  o  1  Qkk 


0  2  °kk 


-it  *x 

e  11  dt* 


(11.12) 
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where  ^  and  t*  are  as  in  (II. 3),  t*  =  -t  and 


J1  =  j  4<~  J  +  ai+a2“y)H]c{y)  e_1  y  dy 
•  00 


(11.13) 


»  »,  ix 

I  =  4(-  i  +  a  +a  -y)H  (y)  e  dy 

2  i.  2  12k 


(11.14) 


I  is  zero  since  it  must  be  integrated  over  the  lower  half  z 

plane  and  H  (y)  has  no  singularities  there, 
k 

1 

We  can  write  I  as 

2 


iT  *y 


,x  _  _  ,  __x  e  '  dy 


I  =  4  (~  -  a  -a  +y) 

2  -®  2  12  (y-ap (y-ap 

and  a  contour  integration  over  the  upper  half  z  plane  yields: 


ia  t*  ia_T* 

i’  =  — -  [  (x-2a  )  e  -  (x-2a  )  e  j 

2  \D  2  1 


(11.15) 


Finally#  a  substitution  of  (11.15)  into  (II.  11)  results 


in 


lar 


ia  t 
2 


;  W  =  U  ‘  ^5[w-2a2)e  *  -(+-2Ve  12 


kk(x)  2  2- 

— ~  OD 


(j;-a^-a^)  (x-2a  )  (x-2a2) 


*  C  (T,t*)dTdt* 
kk 


1  •  • 

2«i  J 

_  Ol 


.  /  ia  t 

iakd  [?  1 


X  a-^  a 2  x*“ 2a^ 


ia. 


X“2a , 


:\k(' 


t*)e 


-ixt 


*d~dt* 


(II-16) 


which  is  identical  to  (I. 10). 

Hence  the  expressions  (II. 1)  and  (II. 2)  are  equivalent 
formulations  of  the  variances  of  the  response. 
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APPENDIX  III 


In  the  following,  the  transform  T (x)  (defined  by  (5.69)) 
is  evaluated-  The  inverse  transform  is  subsequently  obtained 
in  an  approximate  form  by  the  methods  of  section  3.6.  Accord¬ 
ing  to  Eqs.  (5.71)  and  (5.72),  F(x)  can  be  expressed  as 

i/mk 


(JJ)  = 


1/m  . 
3 


2t  (y-ai)(y-a2)  (x-y-a3) (x-y-a4) 


(III-D 


1/m 


dy/m 


(y-a5) (y- 

V 

(x-y-a5)  (x- 

al =  1 

<1  x_  +x.  , 

'k  k  kD 

a2 

=  i£,  'JO,  —  Ju. 

k  k  kD 

a  =  i 

C  .  X  .  - Kt)  .  , 

a  „ 

=  iC  .x  .-'JJ 

3 

]  3  3D 

4 

3  3  3D 

II 

rd 

cn  +  nD 

a 

6 

7  ~  Fp 

Fig.  30 


(III. 2) 


According  to  the  remarks  made 
in  section  4.4,  F (x)  is 

F  (x)  is  the  integral  obtained 
u 

by  a  contour  integration  taken 

over  the  contour  of  Fig.  30(a) 

in  the  limit  as  R  oc. 

1 

Similarly,  7^  (x)  is  obtained 
by  an  integration  in  the  lower 
half  plane  along  the  contour 
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W>  H  > 


of  Fig.  30(b)  in  the  limit  as  oc.  Thus,  replacing  y  in 
(III.l)  by  the  complex  number  z  and  employing  Eqs.  (3.36), 
(3.38)  and  (3.40)  we  obtain: 


f  -1 _ 1 _  _ 1 _ 

u  X  “  11  2i  (x-a  -a  )  (x-a  -a  )  H^al'  (x-a  -a  )  (x-a  -a  ) 
kD  13  14  15  16 


'  £' 


2x. 

kD 

(x-a  -a-. )  (x-a  -a  ) 
2  3  2  4 

1 

1 

2nD 

(x-a5-a3) (x-a5-a4) 

1 

1 

2Hd 

(x"a6"a3) {X'a6'a4) 

1  1 

11 "  2x 

jD  (x-a3~a1 )  (x-a3-a 

1 

1 

2x  . 

3D 

(x-a  -a  ) (x-a  -a  ) 
4  1  4  2 

1 

1  .  1 

2Q 

D 

(x-a  -a  )  (x-a^-a^ ) 

1 

1  .  ] 

2nD 

{X'a6_al)(x"a6-a2) 

H  (a  )  - ± - 

z  (x-a^-a^)  (x-a^-a^) 


Hk(a5> 

W 


6  5' 


®  (x-2a  )  (x-a  -a  ) 

6  6  5 


(III. 3) 


H(a  ) 


3  5‘ 

1 


3  6 ; 


H  (a  )  - 

4  (x-a  -a  ) (x-a  -a  ) 
4  5  4  6 

1 


H  (a  ) 
i  -> 


H.  (a  > 


(x-2a^ )  (x-a6-a5) 


(III. 4) 

Now,  we  examine  the  transform  I (x )  —  [P  (x )  +  P  (x)] 

2  u  * 

=  <  .  Since  I(t)  is,  by  assumption,  a  slowly  varying  function, 

*3 

we  may  approximate  the  transform  of  K  .  (t)  by  the  methods  of 

*3 

section  3.6  so  that  upon  inversion  an  approximately  correct 


form  of  K  (t)  is  obtained.  As  was  stated  in  section  5.6, 
k3 

c,  and  £  are  assumed  to  be  small  (of  order  10  ).  We 

3  k 
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A 

assume  that  the  singularities  of  I(x)  are  no  further  than 

approximately  MAX  [lx  -x  ,  lx*  -Ci  |  #  |  -0  |  }  from  the 

kD  j  D  kD  Dl  1  jD  D1 

origin  of  the  complex  plane  and  that  lx  -x  .  I,  |x\  -ft  |, 

^  1  kD  jD  1  kD  D1' 

|tt>  -ClJ  are  much  smaller  than  the  corresponding  quantities 

lx  +x  |,  |x  +Q  L  |x  +Q  I .  In  other  words,  the  damped 
kD  jD  1  kD  D1  jD  D1 

natural  frequencies  of  both  the  k*1^1  and  normal  coordinates 

are  in  the  vicinity  of  the  centre  frequency  of  the  narrow 
band  filter. 

Consider  the  inverse  Fourier  transform  of  the  first  term 
of  I  (x)T  ( i)  given  by  (III. 3): 

j  1  (»)  r  (^)e  dx  = 

-.00  ^ 


iH  (a  )  cd 
-1  3 


■,J>kD 


2^ 


I  (x)- 


ixt . 
e  dx 


(x-a1-a3) (x-a^-a  ) (x-a^a^ (x-a^a^ 


(III- 5) 


where  from  Eqs.  (II. 2) 


al+ 

a  =  i  ( '  x 

3  3  3 

Z  x  ) 
~k  k 

+  ( X,  +  X 

kD 

V 

a4  “  1<Jjxj 

VV 

+  <\D- * 

a  + 
1 

a5  ‘  i(Ck\ 

cn)  + 

nD» 

a  + 
1 

"e  "  i('kXk 

+ 

-n)  + 

‘Sd-  °d> 

(III .6) 


Using  the  notation  of  section  3.6,  we  can  set 

a  ~  a  (C  j“)j+  :k^k)  +  <\D-  »jD> 
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1 


1 


1 


1 


\ 


G  (a ) 

Jlil  = 

F  (x) 


(a-a1~a3) (a-a1~a5) 


2  a)  (C2  +x  )  4  uu  H 

jD  D  jD  jD  D 


U*1 


(x-a  -a  )  (uu-a  -a  ) 
14  16 


£  «  I X  .  ~  Q  | 

2  1  jD  D1 


0) 

-U)  .  1 

,  |  X  - 

nl } 

kD 

JD 

1  1  kD 

D1 

jD'2nD'l 

“j^D1 

,  I  2x. 

1  /  1  2x. 

(III. 7) 


J 


-  jD  D'  1  kD 

n  =  p+2 
m  =  2 

where  I(x)  is  assumed  to  have  p  singularities  and  where  we 
have  neglected  terms  involving  the  Q's.  Hence: 

(a+A^-a1-a3)  (a+l^-a^a^) 


G  (a )  ^ 

G(a+A^) 

(a-a1~a3)  (a-a^a^ 


(a-a^-a  )  (a-a^a^) 

(a-a  -a  ) (a-a  -a  ) 

13  15 

(  j"2xjD)  j-AjD-“D^  "  2xjD^'XjD+“D^ 


2x  (uu  +0  ) 
jD  jD  D 

l  -  2Md  -rVWV 


2jd  .  (UD  ,  +  n  ) 
JD  JD  D 


j<*»jD-V 

21  . 

3 

2XDD('XjD+aD) 

rw 

UD  +0 
jD  D 
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since  ju  .  ~  Cl  as  stated  earlier.  Therefore, 

jD  D 


MAX 


G.(a.). 


G(a+A  ) 

j 


-  1 


26. 


~  MAX 


2A 


XDD+nD 


2  MAX  tl'*kD-aiDl-KD-nD|] 

jj  .  +0 
jD  D 

Hence,  equations  (3.69)  to  (3.71)  give 


—  I  (*Ju)i  (ju)  e  do) 
2tt  0 


\ 


1 

2tt  J 


-1 


I  (  uu) 


ixt 


where 


2  1®  8x  x  w  (x-aj_-a4)  (x-a^-a^) 

kD  jD  D 

*  (1+1^)  (i+t2) 

p+1  2 

2  6  y  MAX  [  |  J(bk)  |  ] 
o  (T  )  =  o  [  - — 5 - 

1  (p+2)|xj  -n  IyJ  MIN[|j(a  )|3 


dx 


(III  .8) 


r  (hi. 9) 


O  [T  )  <  O  {- 


2  max  [I^D-^DM^D-npl) 

x  +  n 

jd  D 


J 


where  b  are  the  poles  of 
k 


-  and  a  are  the  poles  of  - 

G(x>  k  F(x) 


The  result  (III. 9)  can  be  generalized  to  include  all  the  other 

terms  of  (III. 3)  and  (III. 4).  Let  us  denote  the  combinations 

of  two  a  (see  III. 2)  which  involve  sums  of  two  of  the 
j  <  +  > 

frequencies  uj,  »  UJ  ■  and  Q_  by  a  and  those  combinations 
^  kD  jD  D  k 


172 


It  is  seen  that  each  term  of 


involving  differences  by  a 


(-) 


(  +  ) 


(III. 3)  and  (III. 4)  involves  two  of  the  a^  and  two  of  the 

a^  in  the  denominator.  Hence,  in  each  term,  we  can  let  the 
k 

quantity  involving  the  quantities  a^  be  — —  and  the 

<+)  !  F(X) 

quantity  involving  the  a  be  -  . 

k  G  (uj  ) 

For  the  terms  of  (II. 3)  z=a  is  as  given  in  (III. 7)  and 


4uu  0 
jD  D 


G  (a ) 

and  for  the  terms  of  (III. 4); 


(III. 10) 


a  =  i  (C.tu.  +  Cvtu  )  +  (x  -  x  ) 

j  3  k  k  jD 


G«»  *\0% 


(III. 11) 


And  for  both  (III. 3)  and  (III. 4),  we  may  set: 

e  >  MIN  {  |a*  +  h  a  |  } 

2  1  k  j  1 

61  <  MAX  {  |  a^|  }  |) 


and 


Y  >  MIN  (  |  a  ^  +  t  a  ^  *  |  } 
1  k  j 


'  <  MAX  {  |aj  +  h  a  ^  *  |  } 

2  k  3 


n  =  p+2 
m  =  2 


G  (a ) 


G (a+a  ,  ) 


-  1 


.  (->. 

MAX  {a.  } 

MIN  +  M 
k 


J 


(III. 12) 
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Hence,  the  approximate  expressions  for  F  (r )  and  T0{  ju  ) 

u  *• 


are 

r  (x)  =  - i - [- 

u  8x  x  n 


-H(ax) 


H(a2) 


(x-a  -a  ) (x-a  -a  )  (x-a  -a  )  (x-a  -a  ) 
kDjDD  14  16  23  25 


(*-[• 


Hk<V 


8Q  n  t*-a  -a  )  (’*-a  -a  )  (x-a  -a  )  (x-a  -a  ) 
D  jD  45  56  36  56 


r  (x)  = 

i 


[- 


-H(a3) 


H(a4) 


8x  x  u  (x-a  -a  ) (x-a  -a  )  (x-a  -a  ) (x-a  -a  ) 
kD  ]D  D  23  36  14  45 


[- 


"Hj (a5) 


H.  (a  ) 
3  6 


8:,2n  (x-a  -a  )  (x-a  -a  )  (x-a  -a  )  (x-a  -ac ) 
D  kD  25  56  16  56 


(III.13)& (III. 14) 


which  upon  multiplication  by  I (x ) 


and  inversion  gives  the 


second  moment  response  correct  to  within  the  terms  T  *v-(t) 

1  K  J 

+  T  .  (t)  where 
2  *3 

0  (T  }  < 

/  \  p+1  (  +  )  /  \  2  /  i\ 

2[MAX(|a^  |)]  [MAX(|ak  -a*  ;j)3  MAX[]j(a^  ;)|] _ 

(p+2)rMIN(|a^+)-af+) | ) ] [ MIN ( | a*  + *-a f “ ) | ) ]P+2MIN[ | J (a*~ * ) |] 

2MAX(|F_)  I  ) 

0  (T  }  <  0  {  - - -  ]  (in. 15) 

2  MIN(|a^  +  >|) 
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We  shall  now  obtain  the  inverse  transform  of  F  (x) 

u  x  9 

given  by  (III. 13)  and  (III. 14).  By  the  methods  of  section 
3.5  this  can  easily  be  accomplished.  For  example/  the  inverse 


transform  of  the  first  term  of  (III. 13)  becomes 

M/  v  ixt, 

.  m  .  “H  (a  )  e  duu 

1  r  _ i  1 


2n  8xkD«jDnD  (x-a^-a^)  (x-a^) 


H (a  )  i(Va4H  i(al+a6),; 

1  C  5 -  *  S -  ] 


8x  uo  0 

kD  jD  D 


a  -a 
4  6 


(III. 16) 


a  -a 
6  4 


Doing  this  for  all  the  other  terms  of  (III. 13)  and  (III. 14) 
and  noting  the  relations 


H  (a^ )  =  H*  (a^ )  ,  H^a^  =  H*^)  '  H(a4)  =  H*(a3)  ' 


(III. 17) 


H.  (a  )  =  H.  (a  ) 
D  6  3  5 


and 


(a3-  a5)  =  -(a4-  a^J*  ,  (a^a,.)  =  -(a^)*  )  (III. 18) 


we  obtain,  after  much  algebra: 


1  .  i'Xt 

V  (t)  =  J  ru(x)  e  duj 


u 


8x  x  Q 
kD  JD  D 


[2Re( 


Hfa^e 


i (a  +a  )  t 
1  4 


a  -  a 
4  6 


-)  -  2Re(- 


Hfa^e 


i (a  +a  )t 
1  6 


-)  1 


a  -  a 
4  6 


(equation  continued ) 
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sn2  u>. 

D  jD 


[2  Re  ( 


H(a5)e 


i(a.+a  )t 
4  5 


-)  -  2Re  ( 


1 (a  +a  )  t 

D  O 


)  3 


a  -  a 
4  6 


a  -  a 
4  6 


1  p”  iiut 

vi(t)  =  IT  J  e  ** 


8  x  x  n 

kD  jD  D 


[  2  Re  ( 


H(a3)e 


i(a1+a4)t 


)  -  2Re(- 


H(a3)e 


i (a  +a  ) t 
3  6 


-)  3 


a  -  a 
2  6 


a  -  a^. 
2  6 


1  H  (a  )e 

+  — ^ [2  ReC-3 - 


i (a  +a  ) t  i (a  +a  ) t 

2  5  H .  (a  )  e  5  6 

- )  -  2Re  (-J - 


-)  1 


sn  x'  „ 

D  kD 


a  -  a 
2  6 


a  -  a 
2  6 


And  since  T(x)  =  +  !*.);  Y(t)  =  — (Y  +  i  ,  )  so  that  from  the 

2  u  *  2  u  £ 


above  we  have 


(t)  =  t- 


8x  x  Q  m  m  m 
kD  jD  D  k  j 


[Re  ( 


H (ax )  e 


i  (a1+a4)t 


a  -  a 
4  6 


i(a  +a  )t  i(a  +a  )t 

H  (a  )  e  2  3  H  (a  )  e 

+  Re  ( - - - )  -  Re  ( - - - ) 


a  -  a 
2  6 


a  -  a 
4  6 


-  Re  ( 


H(a3)e 


i(a3«6)t 


•)  1 


a  -  a 
2  6 


[Re  ( 


W 


8x  Q 
jD  D 


a  -  a 
4  6 


-)  -  Re( 


w 


i(a5+a6)t 


■)  1 


a  -  a 
4  6 


(equation  continued) 
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TZ 

8x  n 

kD  D 


[Re( 


H j  (a5)e 


i  (a  +a  )  t 
2  5 


-)  -  Re  (■ 


H.  (a5)e 


i(Va6,,: 


-)  1  1 


a  -  a 
2  6 


a  -  a 
2  6 


[1  +  Tx  +  T21 


(III. 19) 


where  the  orders  of  magnitude  of  T  and  T  are  given  by 

1  * 


Eqs .  (III. 15) 
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